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1.0 INTRODUCTION 

TRW Systems has performed a series of Application Studies for NASA 
Goddard Space Flight Center under contract NAS5-21188. These studies have 
focused on three-axis stabilized control system designs employing momentum/ 
reaction wheel and magnetic or thruster torquing. Tasks included in these 
studies have been conceptual design, system tradeoffs, control law develop- 
ment, performance analysis and simulation. The current study investigates 
two distinct but related aspects of the Earth Viewing Applications Labora- 
tory (EVAL) shuttle missions. The first is evaluating the applicability of 
the gimballed Instrument Pointing System (IPS) to EVAL missions by comparing 
the IPS capabilities with the EVAL requirements. The other area of study 
is assessment of means of stabilizing the shuttle orbiter attitude in earth 
viewing orientations for prolonged periods without use of the orbiter gas 
reaction control system (RCS). 

The IPS, currently under development by the European Space Agency 
for Spacelab, has primarily been considered as a gimballed platform for 
stellar or solar viewing experiments which require pointing to an inertially 
fixed target. One objective of this study is to conduct analyses of an 
earth oriented, shuttle-borne IPS and determine which modifications, if 
any, should be made to the current IPS design to make it suitable to earth 
viewing applications. Shuttle orbits from as low as 150 km up to 1000 km 
are to be considered with small and large .(maximum permissible size/weight) 
IPS mounted payloads. Specific tasks included in this phase of the study 
are: 

• Generation and graphic presentation of parametric data relating 
payload weight and dimensions to the slewed inertia about the IPS 
center of rotation. 

e Analysis of slewing requirements and comparison with existing IPS 
rate and torque capabilities to establish an EVAL payload envelope 
for slewing. 
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• Analysis of dynamic requirements for pointing at earth fixed 

targets and generation of an EVAL payload capability envelope 
for earth pointing. ; 

• Definition of an IPS control sys/tem suitable for the EVAL mission, 
including the attitude reference system and data flow. 

EVAL shuttle missions impose specialized, and somewhat unusual, 
pointing and control requirements on the shuttle orbiter. To minimize 
translational and rotational accelerations on the orbiter, as well as to 
reduce contamination of certain experiment payloads, the orbiter reaction 
control system (RCS) should not be used for orbiter attitude control during 
some phases of EVAL missions. The second main objective of this study is 
to analyze the control requirements for maintaining the orbiter in an earth 
viewing orientation for prolonged periods without use of the RCS, and to 
define an auxiliary orbiter control system capable of meeting these require- 
ments. Specific tasks included in this phase of the study are: 

« Evaluation of potential earth viewing orientations for the orbiter 
with respect to field of view, dynamic stability, and sensitivity 
to system and environmental variations. To minimize the control 
actuator' requirements the desirable orientations are null torque 
attitudes, meaning attitudes at which the total external disturbance 
torque on the orbiter is zero. 

• Selection of one or more orbiter null torque orientations, so that 
both earth surface and earth limb viewing can be performed. Analysis 
of the orbiter auxiliary control system requirements for the selected 
orientations based on the initial conditions imposed by the orbiter 
RCS and the external disturbance environment. 

• Analysis and definition of a suitable set of control actuators to 
stabilize the orbiter near the null torque attitude. Candidate 
actuators to be evaluated include electro-magnets and momentum ex- 
change devices (.reaction wheels or control moment gyros). 
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• Definition of an orbiter null torque attitude stabilization system 
employing the selected set of control actuators. 

Each of- the above defined tasks is addressed in one or more sections 
of the main body of this report. First, however, is a section devoted to 
summarizing the main results and conclusions of the study, and another 
section summarizing the functional requirements and assumptions that underlie 
the study. The data base and supporting details for the study are contained 
in the appendices. - 
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2.0 SUMMARY AND CONCLUSIONS 

2.1 IPS CAPABILITY ASSESSMENT 

Three aspects of the current IPS design limit its earth pointing 
and slew capabilities 

« The gimbal rate limit of P = 2.5 deg/sec imposed by the 
rate gyros. 

9 The maximum gimbal drive motor torque of T m = 20 N-m. 

• The increased value of effective inertia about the gimbal 
center of rotation (CR), due to the payload CG/gimbal CR 
offset (minimum offset is 1.19 m, of which 0.69 m results 
from inside-out gimbal arrangement and 0.50 m results from 
payload support clamp requirements). 

Of the above factors the gimbal rate limit has the smallest impact 
on performance, mainly by setting the lower altitude limit for earth- 
fixed target tracking at h . = 175 Km. The primary restraining 

factor on both tracking and slewing performance is the peak acceleration 
available from the nominal 20 N-m drive motor. As the elevation 
axis has an inertia at least as great as the cross-elevation axis 
and significantly greater rate and acceleration requirements, only 
the elevation axis need be considered in determining the current 
IPS capability envelopes. 

To realistically evaluate the time required to reposition the IPS 
between tracking targets, while including the effect of the gimbal 
angle hard limits and the angle/rate/acceleration boundary conditions 
at the start and end of tracking, the total maneuver between tracking 
targets has been divided into three submaneuvers. For the accelerate/ 
decelerate phases, the worst case condition results when the boundary 
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condition on rate, at the start or end of tracking, is the maximum 
track rate for the particular altitude in question. Table 2-1 sum- 
marizes the IPS accelerate/decelerate capability with the nominal 
motor torque T^ = 20 N-m, under the conditions of allowing a maximum 
time of 10 seconds and a maximum gimbal travel of 10 degrees to go 
from rest to the maximum track rate listed (or vice versa). Figure 2-1 
compares this capability with the required payload envelope, in terms 
of payload length and mass (dashed rectangle). To yield a two-dimensional 
plot, the payload radius has been assumed as r p = 1.0 m. Varying r p 
within the range 0.25 <_ r p <.1.5 m has little effect on the plotted 
curves, however. For each indicated altitude, the region below the 
solid curve represents the payload range that can be accelerated or 
decelerated between rest and the maximum track rate for the altitude 
in less than 10 seconds with less than 10 degrees gimbal travel. -Even 
at the lowest altitude (175 Km)', a large range of payloads can be ac- 
commodated, while at 1000 Km virtually all payloads meet the requirements. 

The nominal IPS slew capability, where the slew angle is measured betwee 
zero gimbal rate points, is summarized in Table 2-2. The requirement 
here is to achieve an average slew rate R<. = Ay/t^ >_ 1.0 deg/sec, with 
the rate limit P = 2.5 deg/sec. Figure 2-2 compares this capability 
with the required payload envelope. For the larger slew angles, where 
the slew time is most important, the nominal IPS is capable of meeting 
the average rate requirement with most payloads. 

With most combinations of slew angle and payload mass properties, 
the gimbal rate limit of P = 2.5 deg/sec is not reached. Increasing 
the value of P therefore will not result in any significant reduction 
in slew time. The only way of decreasing 'slew time while retaining 
the basic IPS structure is to increase the maximum motor torque T^ 
on the elevation drive. Slew performance is relatively satisfactory 
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Table 2-1. Nominal IPS Accelerate/Decelerate to/ from 

Maximum Track Rate Capability for t = 10 sec, Ay = 10 dea 

max 'max y 


Altitude 

h 

(Km) 

Maximum Track Rate 
Y max 

(deg/sec) 

■ 

Maximum Inertia about CR with T = 20 N-m 

m 

(Kg-m 2 ) 

175 

2.49 

3600 

200 

2.16 

4340 

300 

1.41 

8128 

500 

0.869 

13182 

700 

0.554 

20695 

1000 

0.364 

, 

31 466 
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Figure 2-1. Nominal IPS Accelerate/Decelerate Capability vs Requirements With 

t =10 sec and Ay =10 deg Allowed for Reachinq Maximum Track Rate 
max max 
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Table 2-2. Nominal IPS Slew (Zero Rate to Zero Rate)' Capability 
for Average Slew Rate R $ >_ 1 deg/sec with 
Rate Limit P = 2.5 deg/sec 


Slew Angle 
Ay 

(deg) 

Maximum Inertia about CR with T m = 20 N-m 

m 

(Kg-m 2 ) 

20 

5740 

40 

11460 

60 

17200 

80 

22920 

100 

28640 

120 

34380 • 


Table 2-3. Nominal IPS Tracking Capability 


Altitude 
h (Km) 

Maximum Inertia about CR with T = 20 N-m 

? m 

■ (Kg-nr) 

150 

11600 

175 

15911 

200 

20940 

250 

33218 

300 

48556 

500 

143043 

700 

296796 

1000 

657625 
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Figure 2-2. Nominal IPS Slew Capability v$_ Requirements for Average Slew 
Rate R s > 1 Deg/Sec and Rate Limit P = 2.5 Deg/Sec 
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with the nominal motor, and in any case slewing is a secondary considera- 
tion in the IPS design- Increasing the motor torque solely to improve 
slew performance cannot therefore- be justified. 

The primary IPS performance capability required for EVAL is tracking 
earth-fixed targets within acceptable error limits. This is obviously 
contingent on being able to meet the required tracking rates and acce- 
lerations. The tracking rate requirement can be met with the nominal 
rate limit P = 2.5 deg/sec for all altitudes above 175 Km. This- covers 
all anticipated EVAL missions. The available motor torque and required 
acceleration sets the upper limit on the inertia about the CR that 
can follow the track profile. As tracking will be done with a closed 
loop control system, some portion of the motor torque capability must 
be reserved for fine control perturbations and disturbance torques. 

The amount of torque that must be reserved for these ournoses depends 
on the control bandwidth, disturbance environment, etc., but it is in 
general small. Determination of these factors is beyond the scope of 
the present study. Lacking this data, it will be assumed here that 
the full 20 N-m torque is available for accelerating the payload. 

Table 2-3 summarizes the IPS tracking capability, in terms of maximum 
inertia about the CR. The data at h = 150 km is for reference only; 
the gyro rate limit restricts operation to altitudes above 175 km. 

Figure 2-3 translates the data in Table 2-3 into payload length/mass 
terms. The figure clearly shows that earth pointing is feasible for 
all except the largest payloads in the lowest orbits. 

The main attitude determination problem for controlling the IPS 
is that the star tracker view of the celestial sphere is highly re- 
stricted for the EVAL mission by the combined geometry of the orbiter 
payload bay doors and radiators, the orbiter wing surfaces, the or- 
biter tail, the cockpit (crew cabin), the Spacelab pressurized module, 
and the earth. There is a solution to the problem, however, and an 
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Figure 2-3. Nominal IPS Tracking Capability vs^ Requirements 
Assuming 20 N-m Available to Accelerate Payload 
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8° x 40° window of the celestial sphere is available if it is per- 
missible for the tracker field of view (FOV) to come within 20 de- 
grees of the earth's limb and 15 degrees of the orbiter surfaces. 
Depending somewhat on the star tracker aperture, a reasonably sized 
shade in the 15 - 25 inch range can be successfully utilized to per- 
mit the tracker to operate within 20 degrees of the earth's limb. 

The available 8° x 40° window always contains at least one guide 
star of +6 M v or brighter {based on Yale catalog) and with 95% con- 
fidence contains at least one star of +4 M or brighter. The base- 
line IPS star tracker (Honeywell HEAO-B) for solar and stellar mis- 
sions has insufficient track rate capability to permit "on the fly" 
star tracker readings during EVAL payload pointing. Either a tracker 
with a higher track rate is required (up to 0.068 deg/sec for limb 
viewing, up to 2.5 deg/sec to track a point on the earth surface) 
or payload pointing must be interrupted to get a star reading. 

2.2 ORBITER STABILIZATION 

When the combined effects of the external disturbance torques 
and coupled rigid body dynamics are considered, all low altitude 
(h < 400 km) and almost all high altitude null torque attitudes of 
the shuttle suitable for EVAL are unstable equilibriums. An active 
attitude stabilization system must therefore be provided, that can 
satisfy the following requirements: 

® stabilizes orbiter motion 

• able to capture from RCS limit cycle 

a able to tolerate variation and uncertainty in location of 
torque null . 

The two most generally applicable null torque orientations are 
illustrated in Figures 2-4 and 2-5, in both the nominal and maximum 
offset positions. The nose forward orientation is highly unstable 
in pitch, due to gravity gradient torques but provides the best 
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NOMINAL ORIENTATION:, h > 400 km 






Figure 2-4. 


Hose Fomard Orientation 
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view of the earth and has a relatively predictable null torque at- 
titude. The nose down orientation' permits earth limb viewing and 
is only slightly unstable in roll, but the null torque attitude is 
highly sensitive to atmospheric density variations. This is due 
not only to the large cross section of the orbiter in this orienta- 
tion, but also to a near resonance between the orbiter- gravity gra- 
dient libration frequency for the pitch axis and the predominant 
orbit frequency variation in the atmospheric density. 

With the current nominal orbiter parameters, the largest torque 
requirement in most cases is for capture from the'orbiter RCS limit 
cycle, rather than during normal steady state operation. The steady 
state torque requirement is sized largely by the variation and un- 
certainty in the location of the null torque attitude, as determined 
by the atmospheric density (the primary influence at low altitudes) 
and the orbiter mass properties. The dominant influence on the 
transient torque requirement for attitude capture from the RCS limit 
cycle is the nominal limit cycle rate of 0.01 deg/sec rather than 
the nominal limit cycle deadzone of 0.1 deg. The size of the actua- 
tors required for the null torque attitude stabilization system can 
be reduced significantly if an additional "low rate" RCS mode is made 
available. It is recommended that such a mode be provided, with the 
limit cycle rates reduced by at least a factor of 10 from the current 
nominal. A proportional increase in the RCS deadzone, if required 
to achieve the reduced rates, would be acceptable. 

Even if several optimistic assumptions are made, magnets are 
incapable of providing prolonged three axis stabilization for the 
orbiter. Among the points weighing against magnets are the following 

9 Physical limitations of magnetic control due to the relative 
orbi ter/magnetic field geometry introduce undesirable inter- 
axis coupling, with adverse effects on stability. 
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8 


7 

Very large magnets (10 pole-cm per axis) required for even 
marginal performance. 

8 RCS limit cycle rates must be reduced from current nominal 
values by a factor of 25 to 40 to permit attitude capture. 

The feasibility .of doing this is highly doubtful, 

• Lack of reserve torque capability can result in loss of control 
due to variations .in atmospheric density or orbiter mass prop- 
erties . 

• Possible adverse effects of large magnetic field on nearby 
payloads. 

Momentum exchange control, on the other hand, is capable of 
stabilizing the orbiter in either of the selected null torque orienta- 
tions over almost the entire EVAL altitude range. The relatively 
large control torque available allows a more sophisticated control 
system, including adaptive features, than is practical with magnets. 

The adaptive loops can reduce the steady state momentum storage 
requirement by constantly tracking the location of the uncertain 
and variable null torque attitude, and by using the external dis- 
turbance torques for momentum unloading. In order to obtain full 
advantage of this potential reduction in the steady state momentum 
requirement, the transient momentum requirement for capture from the 
RCS limit cycle must be reduced. At least two methods of achieving 
a reduction in the transient momentum storage requirement are possible: 

» Providing a "low rate" RCS limit cycle mode as a transition 
between the normal RCS limit cycle and EVAL operation. 

• Using short (*”50 msec) vernier thruster firings to unload 
excess momentum during the capture transient. 
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If neither of these options ‘is provided, simulation shows that the 
capture transient requires per axis capabilities on the order of 
12 N-m torque and 2700 N-m-sec momentum. 

Two idealized control strategies have been examined for orbiter 
control after the capture transient: 

1. Follow the null torque attitude exactly 

2. Hold the null torque attitude corresponding to the mean 
density. 

Analysis shows that the first strategy requires less torque and 
momentum for disturbance variations with frequencies less than the 
pitch gravity gradient libration frequency, while the second strategy 
should be chosen on this basis for higher frequency disturbance varia- 
tions. That is, the control system should attempt to follow low 
frequency null torque attitude variations, and resist being perturbed 
by high frequency variations. For all altitudes in the EVAL range for 
which aerodynamic torques are significant, the predominant orbit fre- 
quency variation in the atmospheric density falls into the range where 
following the null torque attitude is more efficient. In the nose 
forward orbiter orientation, this implies a peak-to-peak pitch motion 
of less than one degree over an orbit with a +20% density variation. 
This motion is at a frequency far below the IPS bandwidth and is un- 
likely to adversely affect IPS pointing. 

In the "nose down" orbiter orientation a +20% density variation 
results in peak-to-peak null -torque attitude variations of over 30 de- 
grees at 180 km and 10 degrees at 200 km. If the orbiter attempts to 
follow the null torque attitude, the peak pitch rate relative to the 

“ i 

orbit reference frame is 0.017 deg/sec at 180 km and 0.005 deg/sec 
at 200 km. If the orbiter attempts to maintain a fixed attitude 
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relative to the orbit reference frame instead, the peak-to-peak mo- 
mentum excursion would be over 6000 N-m-sec at 180 km and 3600 N-m 
sec at 200. km. These values exceed the momentum storage capability 
required to capture from a 0.001 deg/sec RCS limit cycle, and could 
thus size the CMG's. Expected atmospheric density variations in the 
"nose down" orientation result in either excessive attitude excur- 
sions o.r extreme momentum storage requirements, and therefore may 
make operation in this orientation unfeasible below about 250 km. 

Control compensation has been designed to use the available 
measurements to approximate the selected control strategy. The mo- 
mentum' required to accommodate the orbit period atmospheric density 
variations with the practical control law is about twice that re- 
quired with ideal control. 

A tradeoff study has indicated that double gimbal control mo- 
ment gyros (DG CMG) are the most suitable momentum exchange devices 
for the orbiter null torque attitude stabilization system. Three 
orthogonally mounted DG CMG’s, each capable of producing 2710 N-m- 
sec (2000 ft-lb-sec) momentum and 13.6 N-m (10 ft-lb) torque, will 
allow capture from the nominal 0.01 deg/sec RCS limit cycle and pro- 
vide orbiter attitude stabilization in the nose forward orientation 
above 180 km and in the nose down orientation above 200 km. Smaller 
CMG’s could be used if the recommended reduction in RCS limit cycle 
rates is provided and if operation in the nose down orientation is 
restricted to altitudes above 250 km. 

2.3 SYSTEM DEFINITION 

Figure 2-6 shows the interfaces between the major subsystems 
and components comprising the EVAL pointing and control system. 

Four subsystems participate in the pointing and control task 

© the orbiter flight control system (FCS) 

e Space! ab 

e Orbiter null torque attitude stabilization system 

• IPS control system 
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Figure 2-6. IPS/EVAL Pointing and Control System Interfaces 

















The first two of these subsystems fall, for the most part, outside 
the scope of the present study and will' be discussed only in terms of 
functional requirements implied by the study results. The latter two 
subsystems are peculiar to EVAL, and detailed data flow diagrams for 
each are presented and discussed in Sections 7 and 12 respectively. 

Several steps have been taken to simplify the interfaces, minimize 
the amount of data that must be interchanged, and minimize the compu- 
tational burden on the spacelab experiment computer. For uniformity 
in the interfaces, al 1 pointing and control data transfer is performed 
digitally over data buses. Remote acquisition units (RAU's), described 
in Reference 1, are used to perform the necessary addressing and data 
formatting functions. The pointing and control data flow through the 
experiment data bus and the computational burden on the spacelab ex- 
periment computer has been minimized by performing all high rate data 
processing locally within the subsystems, using programmable digital 
electronics (PDE). For definiteness, Figure 2-6 indicates the PDE 
as mini -computers, but a future trade study is required to choose 
between mini -computers, general purpose microprocessors, dedicated 
special purpose microprocessors, or some combination of the above. In 
any case, the experiment data bus is only required to handle the following 
low rate data 

s pointing commands (i .e. , target coordinates) 

9 orbiter attitude data 

s orbiter ephemeris data 

» mode control commands 

6 housekeeping data 
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The guidance, navigation, and control (GN&C) and reaction control 
system (RCS) elements of the orbiter flight control system are utilized 
by the EVAL pointing and control system. The required operating modes 
and GN&C system outputs are 'for the most part included in the baseline 
system described in Reference 1. The attitude and ephemeris data listed 
in Table 2-4 (or its equivalent) is required for EVAL and is assumed 
available from the GN&C computer on an external data bus. The trans- 
formation from the earth centered inertial (ECI) to orbit reference frame 
(R), can be derived from the ephemeris data if not computed in the 
GN&C computer. The orbiter null torque attitude stabilization system 
requires A RI plus the orbiter attitude reference data oig and Agj. The 
IPS control, system includes its own inertial attitude reference system 
and requires attitude data from the orbiter only for initialization. 
Ephemeris data is required for IPS pointing at and tracking of earth 
fixed targets. The orbiter null torque attitude stabilization system 
supercedes the orbiter reaction control system in controlling the orbiter 
attitude during EVAL operation requiring either low accelerations or low 
contamination levels. The RCS is, however, needed to support the null 
torque attitude stabilization system by providing suitable initial con- 
ditions on the orbiter state vector and by producing torques to unload 
excess CMG momentum. The modes of RCS operation required are: 

• The existing rotational maneuver mode is needed to initialize 
the orbiter in the estimated null torque attitude. 

9 The existing attitude hold mode should be modified to provide 
a 0.001 deg/sec limit cycle about the estimated null torque 
attitude. 

• A new momentum unloading mode is required for minimum on- 
time 40 msec) vernier thruster firings to desaturate the 
CMG cluster. This mode is used during initial capture 
from the RCS attitude hold limit cycle and for occasional 
unloading of secular torque build-up. 
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Table 2-4. Attitude and Ephemeris Data Required 
from Orbiter GN&C Computer 


Symbol 

Description 

“b 

Orbiter body components -of orbiter inertial rate 

a bi 

Direction cosine matrix indicating the attitude 
of the orbiter body (GN&C base) relative to the 
ECI frame 

a ri 

Direction cosine matrix indicating the attitude 
of the orbit reference frame relative to the 
ECI frame 

? I 

Orbiter position in ECI coordinate 

V I 

Orbiter velocity in ECI coordinate 

? I 

Orbiter acceleration in ECI coordinates 

GMT 

Greenwich mean time 






Overall operation of the EVAL pointing and control system is 
controlled and coordinated by equipment in the Spacelab module. This 
equipment consists of display and control panels for the IPS and 
for the orbiter null torque attitude stabilization system, the exper- 
iment computer, and the experiment I/O unit. To minimize the burden 
on the experiment computer, processing in the spacelab is limited to 
generation of target commands and transfer of attitude and ephemeris 
data. All high rate loops should be closed, locally in the IPS control 
and orbiter null torque attitude stabilization system computers. 


2-20 



3.0 FUNCTIONAL REQUIREMENTS 

This section briefly summarizes the EVAL functional requirements 
relevant to this study. Requirements relating to control of the IPS are 
presented first, followed by orbiter stabilization requirements.. 

3.1 IPS REQUIREMENTS 

Table 3-1 summarizes the IPS/payload constraints and requirements 
for EVAL that have been identified. The values given in the table for the 
performance/capability constraining parameters are the nominal values for 
the baseline IPS design; modifications may be necessary or desirable to 
adapt IPS to EVAL. 

The requirement side of the picture is .less clear. Data on specific 
payload/orbit/viewing requirement combinations for potential EVAL missions 
is sparse or lacking. For this reason, and also to provide results with 
wider applicability, the assumption is made that tracking of earth-fixed 
targets is required over the entire 60 degree half angle view cone, over 
the indicated orbit range, for all combinations of payload mass properties 
satisfying the baseline IPS payload constraints. No data has been avail- 
able on specific slew time/angle requirements; Table 3-1 therefore gives 
assumed values. 

3.2 ORBITER STABILIZATION REQUIREMENTS 

The orbiter attitude is normally stabilized using the vernier reac- 
tion control system (RCS) thrusters described in Appendix A. To reduce 
contamination and translational accelerations to the level required by some 
EVAL missions, the RCS must be disabled and the orbiter stabilized by auxil- 
iary means. Table 3-2 summarizes the functional requirements for orbiter 
stabilization with the RCS disabled. 
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Table 3-1. IPS Constraints and Requirements for EVAL 


" ' ' “""I 

Parameter 

Constraint/ 

Requirement 

Symbol- 

Val ue or 
Range 

Comments 

Drive motor torque 

C 

T 

m 

20 N-m 

Both motors in drive used 
simultaneously 

Maximum slew rate 

c 

P 

2.5 deg/sec 

Limited by maximum input rate 
permitted by rate gyros 

Payload CG/Gimbal CR 
Offset 

c 

- 

1 .19+3.69 m 

] Payload constraints 

Payload length 

C/R 

Up 

1 + 6 m 

( dictated by payload sup- 

Payload Mass 

C 

m p 

3000 kg 

1 port clamp requirements 


R 


<_ 2000 kg 


Payload radius 

C/R 

r P 

0.25+1.5 m 


Minimum inertia about 
CR (no payload) 

C 

T CR 

56 kg-m^ 

Estimated from outer gimbal and 
payload integration ring mass 
properties 

Gimbal motion limits 

C 

l^max, i y lmax 

70 deg 

Limited mechanically in gimbal 
assembly and electronically in 
drive electronics 

View cone half-angle 

R 


60 deg 

Target tracking required in 
view cone, centered on local 
vertical 

Target tracking error 

R 


N/A 

Largely determined by gyro, . 
star tracker and ephemeris 
errors 


* CR = center of rotation 



Table 3-1. IPS Constraints 


Parameter 

Constraint/ 

Requirement 

Shuttle orbiter altitude 

R 

Orbit inclination 

. 

R 

Average "Slew" Rate 

R 

Maximum Accelerate/ 
Decelerate Time 

R 

Maximum Accelerate/ 
Decelerate Angle 

R 


* Assumed values 




Table 3-2. Orbiter Stabilization Requirements for EVAL With RCS Disabled 


Requi remen t 

Value or Range 

Comments 

Shuttle orbiter altitude 
Orbit inclination 

150 to 1000 Km 
30 to 60 deg 

Missions below 200 Km unlikely 

Some polar orbits probable, 
equatorial orbits unlikely 

Two orbiter orientations to provide following 
fields of view from payload bay: 

e Earth surface viewing 

e Earth limb viewing 

nadir to horizon 
>120 deg sector 

Subject to limitations imposed 
by orbiter structural 'obstructions 

Candidate actuators to stabilize orbiter when 
reaction control system (RCS) cannot be used 

• Electro-magnets 
e Momentum exchange devices* 


Preferred if feasible due to lower 
cost 

Acceptable if electro-magnet con- 
trol is not feasible 

Initial conditions on magnet or momentum ex- 
change device control system 

• Maximum offset from null torque attitude 

• Maximum angular rate 

0.1 deg 
0.01 deg/sec 

Worst case vernier RCS limit cycle 
conditions 


* Reaction wheels or control moment gyros 






4.0 IPS GIMBAL/PAYLOAD MASS PROPERTIES 

This section contains a derivation and numerical results on 
the possible range of inertias, about the gimbal center of rotation 
(CR), that can be expected for IPS payloads. The IPS gimbal arrange- 
ment and allowable payload envelope, based on Appendix B, is presented 
in Figure 4-1 and Table 4-1. Figure 4-2 shows the nominal payload 
orientation relative to the shuttle orbiter payload bay. The following 
assumptions about the payload have been made to allow the results to be 
presented parametrically: 

• The payload is a uniform right circular cylinder of mass m , 
length £p, and radius r . 

• The payload CG is offset from the CR only along the z^ axis. 

A complete mass model of structures composing the outer gimbal, 
payload integration ring, gyro package, and data electronics is not 
available. This structure is therefore conservatively modeled as a 
thin circular ring of radius r^, containing the total estimated mass, 
located at the front surface of the payload integration ring. Since 
most of the. inertia about the CR is due to the payload mass and CG/CR 
offset (in most cases), the final results are relatively insensitive 
to the above assumptions. 


4.1 DERIVATION OF INERTIA ABOUT CR 

The total inertia about the CR can' be expressed in dyadic form 
(Reference 2) as' 


yCR _ yCR , -yCR 

n x p + 1 m 
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Table 4-1. Summary of Gimba4 /Payload Envelope 


Symbol 

Definition 

Range or Value 

CG 

Payload CG 


CR 

Gimbal Center of Rotation 

- 

3 

Cross-Elevation Angle 

- 70 £ 3 £ 70 deg 

Y 

Elevation Angle 

- 70 £ y £ 90 deg 

a 

0 

CR/Payload Integration Ring Offset 

0.69 m 

l P 

Payload Length 

1 + 6 m 

r P 

Payload Radius 

0.25 ->1.5 pi 

mp 

Payload Mass 

0 ^ 2000 kg 

r M 

Payload Integration Ring Radius 

0.6 m 


Total Mass of Payload Integration 
Ring, Outer Gimbal, Gyro Package 
and Data Electronics 

~ 85 kg 
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where 


l P = X P + m p ^ r CG ’ r CG^ 1 " r CG r CG^ 

^ ° T M + > [(f M • 7 H> ” - F M F M ] ' 

T = Unit dyadic 

Expressed in the p-frame, the payload inertia dyadic about its CG is 


jCGtCG +jCG + i c ® 7 2 

y L xx p p A yy y P y P zz z p p 


where for a right circular cylinder 


CG i 2 
! xx = 2 m P r P 


CG nip p P 

! yy s 12 r P + V 

CG nip p p 

I- 4(3 rf+O 


l zz 12 


The CG/CR offset vector is 


r 0G = ~ + ^p/2) x 


o A P 


and therefore 
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y€R _ /1_ J - \ " v 4. 

*p (g ^p ' 


.[ 


m 


(3^p + ^p) + nip (£ 


Lit (3fp- + ^p) + m p (» 0 + V 2)2 ] ^ ^ 
i*„ + V 2 > 2 ] 


Up z p 


r 2 , 2 


n r 

["p (-£ + 4 

+ ‘o + *o *P> 

/> A 

X E X E + r 


i-p «p 2 

'p<T + TF + *o + % l P 


.y e y E + tf ">p >-p) * E 2 e 


The "payload integration ring" inertia about its CG is 


tM. t h : + r M ; : + x m' ^ z 
*xx x m *m + L yy y M y M L ZZ Z M Z M 


where for a thin circular ring 


t M _ 2 

*xx ” 2 r M 


t M = H!m „2 

yy 2 r M 


4 = m M r M 
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The CG/CR offset is 


r M = ~ £ o Z M 


and therefore 


p o r\ y*. a. rn v. n O ^ a 0 ^ ^ 

r M = ^2“ r M + m N £ 0^ X M X M + ^2“ r M + m M £ o^ y M y M + ^ m M V Z M Z M 


- [l»H Q + l o>] *E *E + [ m M <!T + l o>] h *E + tm M r M> 


Z E Z E 


The total inertia about the CR is finally 


tCR 


A A 


A A 


lj - I 1 x E x E + 1 } y E y E + I 2 z E z E 


where 




l 2 2 m P r P + m M r M 
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The elevation motor axis is along yjr and therefore the inertia 
that must be slewed by the elevation torquer is 


I = I, 
TY 1 


The cross elevation axis is along Using the unit vector relations 
(Appendix B) 


x c = cos v x v - sin y 2 , 


y E = y X 


= sin y X x + cos Y z x 


the total inertia can be expressed in the X-frame as 


I^ R = (1-j cos 2 y + ^ sip2 t) x X X x + *1 y X y X 


( I 2 - I-j) sin y cos y x Y z Y + [l 0 - I-,) sin y cos y z v x 


X X Vi 2 A 1 


-X ~X 


2 o 

+ (I 2 cos y + I-j sin y) z x z X 


The inertia seen by the cross elevation torquer (neglecting the 
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insignificant inertia of the cross-elevation rotor) is thus 


= I, cos y + lo sin y 


4.2 NUMERICAL RESULTS 

All payloads in the allowable range result in I-j > The 
maximum value of Igg as a function of y is therefore I-j. Since I 
is also equal to I-j, only Ij need be considered further. 

Using the values in Table 4-1, the constant non-payload contribution 
to I-| is 


r 2 

'Mo o 

m D (4 1 + l d ) = 55.7685 kg-nT 
R 2 0 


while the payload contributes 


m r 


rrp ip 2 

T + T + "0 £ P + £ o 


= m, 


r r 2 ip 

+ -y + 0.69 i p + 0.4761 


Therefore I-j is, as a function of the payload properties nip, i p , and 


I-j = 55.7685 + mp 



+ 0.69 z n + 0.4761 
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The minimum value of 1-j occurs for m p = 0 

<Vn.in = 56 ks - mZ 


while the maximum value for the specified payload envelope occurs 
for 'iDp = 2000 kg, r p = 1.5m, t p = 6m 

34413 ^ 


Figures 4-3 to 4-5 contain plots of- 1^ vs. £ p with the limiting 
values of r p (.25 and 1.5 m) for three values of m p (10, 200, 3000 kg). 
It is clear from these figures that 1^ is relatively insensitive to 
r p within the allowed range. This shows that the exact radial mass 
distribution of the payload has little effect on the final inertia value 
and therefore the assumption of a uniform cylindrical payload .does not 
unduly restrict the applicability of the results. All following plots 
assume r p = 1 m. 

Figure 4-6 to 4-8 contain plots of I-j vy £ p for selected values 
of payload mass m p ranging from 10 to 3000 kg. Curves of constant I-j 
in the range 60 to 30000 kg-m 2 are plotted versus payload length £ p 
and mass m p in Figures 4-9 to 4-14. These curves are useful in 
establishing the ‘bounds on payload length and mass for a given 
inertia about the CR. 
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Figure 4-4. Inertia About Gimbal CR, m p = 200 Kg 
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Figure 4-5. Inertia About Gimbal OR, m p - 3000 Kg 
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Figure 4-6. Inertia About Gimbal 
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Figure 4-7. Inertia About Gimbal CR, m = 100, 250, 500 Kg 
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Figure 4 - 8 '. Inertia About CR, m = 1000, 2000, 3000 Kg 





PRYLORD MFISS (KG) 




PAYLOAD MASS (KG) 



Figure 4- .10. Payload Length and Mass for Inertia About CR of 100, 200, 300 Kg-rn^ 
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igure 4-12. Payload Length and Mass for Inertia About CR of 1000, 2000, 3000 Kg-m 2 
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Figure 4-13. Payload Length and Mass for Inertia About CR of 4000, 6000, 8000 Kg-m‘ 
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5.0 IPS SLEWING ANALYSIS 

The + 60 degree half-angle visibility cone about nadir of the IPS 
restricts the viewing time of a given earth-fixed target to a relati- 
vely small fraction of the time the target is above the horizon. For 
example, at 150 Km altitude, a target on the ground track is above the 
horizon for 179 seconds on each side of nadir, but is within the + 60 
degree half-angle view cone for only 35 seconds on each side of • 
nadir (see Figure 6-6). Therefore, efficient experiment utilization 
requires that the slew time between target locks be minimized. In 
general, the gimbals must already be moving when the target enters 
the field of view, and are still moving when the target leaves the 
field of view. The direction of gimbal motion may therefore have to 
reverse twi ce between loosing one target and acquiring the next. 

This situation is illustrated schematically in Figure 5-1. The ele- 
vation and cross-elevation axes are independent and do not necessa- 
rily reverse direction at the same time. 

Although the entire sequence of events between loosing track 
with one target and starting to track the next is essentially one 
maneuver, it is convenient for the present purposes to divide the 
maneuver into three phases: decelerate, slew, and accelerate. This 

not only simplifies the analysis, but also provides a means for in- 
cluding the gimbal angle hard limits. For IPS, the "buffer zone" 
between the gimbal angle hard limits and the gimbal angle limits for 
target tracking extends from 60 deg <_ jy| £70 deg and 60 deg <_ [g| £ 70 deg. 

In addition to the above, three other factors also limit the 
slew response for IPS 

# The rate gyro limit of 2.5 deg/sec. 

« The acceleration limit imposed by the available motor torque 
and slewed inertia. 

© The softmount dynamic response. 
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The softmount dynamic response is beyond the scope of the present 
phase of the study. The remaining factors are analyzed in the following 
sections, with the results presented as parametric plots. 


5.1 TIME/ANGLE RELATIONS FOR ACCELERATE/DECELERATE PHASES 

During the accelerate/ decelerate phase it is required to bring 
the gimbal axes from/to rest to/ from the required track rate at the 
start/end of a track maneuver. The motions of the two gimbal axes 
are essentially independent, but must be coordinated to reach the 
correct vector rates. The means of achieving the required coordination 
are beyond the current scope. The coordination requirement does, 
however, limit the response capability to the weaker of the two axes. 

For IPS, this is the elevation axis. Compared to the cross-elevation 
axis, the elevation axis has greater or equal slewed inertia for all 
payloads and gimbal angles, and maximum tracking rate requirements 
an order of magnitude higher. Therefore only -the elevation axis need be 
considered. 

Lack of detailed information about the IPS gimbal drives requires 
that the following reasonable assumptions be made: 

e The maximum drive motor torque T m is independent of motor 
speed over the range of interest. 

® The time required to change the drive motor torque is small 
compared to times of interest. 

Under these assumptions, the minimum time required to accelerate from 
rest to the required track rate results from application of full torque 
until the rate is reached. .The resulting acceleration is 
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v - V 1 ! 


while the rate is governed by 


' ' tf h 


and the change in angle is 


4 v - T m 


The deceleration equations are similar, with obvious minor changes. 

The maximum rate that can be reached in time t within an angular 
change Ay is 


KW = 2 4y/t 

More useful for evaluating capabilities are the relations involving 
the slewed inertia and motor torque 


^max (yTj 

1/2 
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The above two families of curves have been plotted as functions of 

the ratio (I-j/T ) in Figure 5-2. To illustrate some uses of these 

.curves, Figure 5-3 is presented as an example. The spotted region 

of this plot indicates the range of rates that can be reached in 

less than 10 seconds with a change of gimbal angle of less than 

10 degrees. If the motor torque T m is 20 N-m and i.t is' required to 

reach a rate of 2.5 deg/sec within the given constraints, the maximum 

slewed inertia is found to be 20 x 180 = 3600 kg-m . Interpolating 

between the upper curve of Figure 4-12 and the lower curve of Figure 

4-13 shows 'this to represent a moderate size payload. As another 

example, suppose again that T = 20 N-m but now there is a larger payload 

o m 

with I-j = 8000 kg-nr. The plot shows that the maximum rate achievable 
within the given constraints is 1.45 deg/sec. As will be shown in 
Figure 6-3, this represents the peak tracking rate requirement at a 
'240 km altitude. 


5.2 TIME/ANGLE RELATIONS FOR SLEWING 

■ During the slew phase, it is required to bring the gimbal angle 
from rest at y = y Q to rest at y = y^ in minimum, time subject to . 
limitations on the gimbal acceleration and rate. The gimbal acce- 
leration is limited physically by the torque/inertia ratio T^/I-p The - 
gimbal rate for IPS is limited by the drive electronics to P = 2.5 deg/sec 
so as not to exceed the allowable gyro input rates. The' rate limit 
is considered a parameter in the following analysis, however, to permit 
evaluation of the impact of the rate limit on overall performance. 

For a minimum-time reposition with rate and acceleration limits, • 
either the rate or acceleration must be at a limiting value at all 
times (Reference 3). Two cases exist here, depending on whether the . 
rate limit is reached, as shown in Figure 5-4. In, case 1, maximum 
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Figure 5-2. Time/Angle/Rate Curves for Acceleration/Deceleration Phases 
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Figure 5-3. Examples With t v = 10 sec. Ay = 10 deg 
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torque is applied until the rate limit is reached at 


4 ■ *1 ■ P V T m 


The motor is then turned off (neglecting friction), resulting in 
constant rate P until 


t “ t£ “ t^T — t-j 


when the maximum torque is applied in the opposite direction. The 
change in angle 


Ay = Yf - y c 


from t = 0 to t = tr is equal to the area under the rate curve 


Ay = Pt 2 = P(t f - t,) = P(t f - PI i /T m ) 


Solving the above for the slew time gives 


, PI, 

+ _ kL + I 

t f P + T 

m 
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Case 1 applies when 


tf >_ 2 t-j 


or 


Ay > 


iip 2 

T m 


Otherwise, the rate limit is not reached and case 2 applies, 
the torque switches from maximum one direction to maximum the 
direction at 


t = t 3 = t f /2 


The maximum rate reached is 


Y 


max 


■ T m V 1 ! 


The area under the rate curv.e is then 


Ay = 



^ax *3 * 'Vh - T m 4' 4T 1 


Here, 

opposite 
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so the slew time is 


*f ■ 2 "i ^ 


1/2 


Slew time, as a function of inertia/torque ratio, is plotted in 
Figures '5-5 to 5-10 for selected values of Ay in the range of 20 to 
120 degrees, with the' rate limit P as a parameter. It is apparent 
from these plots that the slew time is relatively insensitive to the 
rate limit for all but the smallest payloads. The IPS rate limit of 
2.5 deg/sec does not unduly restrict the slew time and any increase 
would produce only a marginal improvement. 

Figures 5-11 to 5-15 contain the same data, plotted for fixed 
values of the rate limit with the slew angle as the parameter. 
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Figure 5-6. Slew Time for 40 Degree Slew Angle 
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Figure 5-7. Slew Time for 60 Degree Slew Angle 
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Figure 5-9. Slew Time for 100 Degree Slew Angle 
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Figure 5-10. Slew Time for 120 Degree Slew Angle 
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Figure 5-11. Slew Time for Rate 
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P = 2.0 Deg/Sec 














Figure 5-12. Slew Time for Rate Limit P = 2.5 Deg/Sec 
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Figure 5-13. Slew Time for Rate Limit P = 3.0 Deg/Sec 
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Figure 5-14. Slew Time for Rate Limit P = 4.0 Deg/Sec 
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Figure 5-15. Slew Time for Rate Limit P = 5.0 Deg/Sec 
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6.0 IPS EARTH POINTING ANALYSIS 


This section contains the general relationships for the gimbal 
angles, rates, and accelerations required to maintain the LOS aimed 
at a fixed point on earth. The special case of the aim point in the 
orbit plane is considered first. This case represents the worst case 
condition for the elevation axis and yields closed form expressions 
for the maximum elevation gimbal rates and accelerations as a function 
of altitude. The case of the general aim point is considered next, 
to determine requirements on the cross-elevation axis. The complexity 
of the resulting relationships precludes an analytical determination 
of the maximum cross-elevation rates and accelerations. A computer 
program was therefore written to evaluate the gimbal angle time 
histories for the general aim point case. A number of sample runs 
indicates the peak cross-elevation rate and acceleration requirements 
are about an order of magnitude less than the elevation requirements. 
Finally, the relationship between earth pointing and payload capabi- 
lity is presented. - 


6.1 AIM POINT IN ORBIT PLANE 

A simple but enlightening case results when the aim point is on 
the intersection of the earth surface and the orbit plane and earth 
rotation is neglected. This case is illustrated in Figure 6-1. The 
shuttle orbiter is assumed to be flying in an inverted orientation, 
with the payload bay "down" (xg along the velocity vector and 
pointing "up" along the local vertical). In this configuration, the 
IPS elevation gimbal angle y, to point the LOS at the aim point, is 
as shown in the figure- The elevation angle relationship is easily 
derived by considering the two right triangles indicated by the dashed 
line. 
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r s - |r s ( = orbit radius 
r Q - |r^j = earth radius 

p = I p"I = ^ r s + r o “ 2 r S r o cos v ^ 1/2 = $/C-Aim Point Distance 

V - V + U) t 
0 0 


Figure 6-1. Aim Point in Orbit Plane 
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The odd derivatives of y are even functions of t (with v Q = 0) 
and reach their maximum values at t = 0, i.e., when pointing at nadir. 

These maximum values are 


y = y{0) 
'max 1 v ' 


o o 


Y max = - y(0) 


u n r -r 
_ 0 0 s 


(r s + r 0 ) 


h' 


where 


h = r$ - r Q = S/C altitude 


The maximum values of y and y occur at non-zero times determined 
by setting the next higher derivative to zero. These maximum values 
are 


Y max = ? 


>)] 

• -1 
= sin 

' r o ] 

> J 


Lr + h J 


y = y(t = t ) 
'max rv m y 


2 / 2 2 \ . 
oj - r) sin m t 

o o s s o o m 

(r 2 + r 2 - 2 r r„ cos or t m ) 2 
v s o o s o nr 


where 


t . -l cos’ 1 1 [ 2 ♦ (dii-JT- iiAl 

” “o U ' 4r o r s / J 4 r o V 
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When y = v , the LOS is tangent to the earth at the aim point. For 
most payloads, the target is effectively lost before this angle 
is reached. 

Figures 6-2 to 6-5 contain plots of these maximum values as' a 

function of altitude. Figures 6-6 and 6-7 show time plots of the 

elevation angle and its derivatives for a low (150 km) and high 

(700 km) altitude orbit. Note that for the 150 km case, the rate and 

acceleration are still significant when the elevation angle reaches 

60 degrees. Also note that the maximum value of jerk (y) required 

for tracking an earth fixed-target is relatively low, about 0.016 
3 . 

deg/sec . The IPS gimbal rate limit of 2.5 deg/sec is only exceeded 
for altitudes below 175 km. No earth pointing missions below 175 km 
are anticipated however. 


6.2 .GENERAL AIM POINT 

When earth rotation and aim points displaced from the orbit plane 
are considered, the complexity of the gimbal angle relationships for 
target tracking is vastly increased. These relationships are presented 
in Appendix C. A computer program for generating the gimbal angle 
time histories and peak values, for a general aim point on a rotating 
earth, is also described in Appendix C. A large number of cases were run 
with this program, one of which is contained in Appendix C. The results 
of this series of runs can be summarized as follows: 

e The relations in Section 6.1 accurately describe the elevation 
axis requirements. 

% The cross-elevation axis rate and acceleration requirements for 
target tracking are no more than 10% of those for the elevation 
axis (azimuth axis not used). 
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Figure 6-7. Altitude = 700 KM 
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TARGET ON HORIZON 
AT t = 423 SECONDS 







The second point leads to the conclusion that the cross-elevation 
axis capability is not the limiting factor in earth pointing for 
IPS and need not be further considered in this regard. 


6.3 EARTH POINTING/PAYLOAD CAPABILITY RELATIONS 

The maximum inertia about the gimbal CR that the IPS is capable 
of maintaining pointed at an earth-fixed aim point at S/C altitude 
h is 


(h) 


1 'max 


* yw h) 


where the expression for v as a function of h was presented in 

max 

Section 6.1. The maximum gimbal rate required for tracking has also 
been expressed and plotted as a function of altitude. Both of the 
above relations are combined in Figure 6-8 to yield the overall 
relationship among altitude, gimbal rate limit, motor torque, and 
maximum inertia about the gimbal CR. 

Figure 6-9 illustrates the use of these curves by indicating the 
range of allowable inertias about the CR (shaded area) using the 
nominal IPS parameter values T m = 20 N-m and P = 2.5 deg/sec. The rate 
limit sets the lower limit on altitude (175 Km) while the motor torque 
sets the upper limit on inertia. For any given altitude the maximum 
inertia about the CR can be read off. The inertia about the CR can 
then be translated into ranges of payload mass and length using Figures 4-9 
to 4-14. 
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Figure 6-8. Maximum Inertia About CR for Earth Tracking VS Altitude 
T m - 10, 15, 20 N-m 


6-12 



MRXIMUM INERTIR (KG-Mk*2) 














7.0 -IPS .CONTROL SYSTEM DEFINITION 

This section defines an IPS control system configuration suitable 
for the EVAL mission. The IPS control system, shown in block form in 
Figure 7-1, is configured to be as nearly autonomous as possible. In- 
puts required from the space! ab experiment computer include only tar- 
get commands and orb iter ephemeris data, and all control loops are 
closed locally on the pallet with programmable digital electronics (PDE). 
The computational burden on the space! ab experiment computer and the 

transfer of control related data through the experiment data bus are 
thereby minimized. 

The major blocks of the IPS control system are described in the 
following subsections. Figure 7-2, illustrates the IPS gimbal struc- 
ture and coordinate frames. 

7.1 ATTITUDE REFERENCE SYSTEM 

A functional block diagram of the IPS ARS is shown in Figure 7-3. 

The ARS inertial sensors, gyros and a star tracker, are mounted along 
with the payload on the IPS platform. Processing of the sensor data is 
performed on the pallet, with PDE. A high bandwidth, stabilized attitude 
reference is obtained by integrating processed gyro rates with a closed 
form quaternion algorithm (Reference 4). The gyro processing consists 
of correcting for gyro biases, misalignments and scale factor errors. 

The star tracker data is processed to correct for known errors, 
compared to data in the star catalog for identification, and combined 
with estimated platform attitude to yield the attitude residual. The 
attitude residual is processed by an extended Kalman filtering algorithm 
to compute optimal attitude and gyro bias, updates for long term attitude 
reference stability. The frequency of the required update depends on 
the quality of the system components and the required attitude deter- 
mination accuracy. 
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Figure 7-1. IPS Control System Functional Block Diagram 
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Figure 7-2. IPS Coordinate Frames 
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Figure 7-3. IPS Inertial Attitude Reference System functional Block Diagram 










The main attitude determination problem is that the star tracker view 
of the celestial sphere is highly restricted for the EVAL mission by the 
combined geometry of the orbiter payload bay doors and radiators, the or- 
biter wing surfaces, the orbiter tail, the cockpit (crew cabin), the 
Spacelab pressurized module, and the earth. Figure 7-4 illustrates this 
restrictive geometry. The figure also indicates that there is a solution 
to the problem, and that an 8° x 40° window of the celestial sphere is 
available if it is permissible for the tracker field of view'(FOV) to 
come within 20 degrees of the earth's limb and 15 degrees of the orbiter 
surfaces. Depending somewhat on the star tracker aperture, a reasonably 
sized shade in the 15-25 inch range can be successfully utilized to per- 
mit the tracker to operate within 20 degrees of the earth's limb. 

Since the window on the celestial sphere available to the star . 
tracker is limited to a region comprising only 8° x 40° = 320 square 
degrees, star availability must be examined. Figure 7-5 plots the width 
of one side of the square field of view required on the celestial sphere 
versus star visual magnitude to’ assure with 95% confidence the presence 
of at least one star of the indicated brightness, or brighter. For ex- 
ample, if one can detect stars up to M = +6, then with 95% confidence 
one can find a usable star within any 5x5 = 25 square degrees of the 
celestial sphere. The actual worst case. is also plotted using Yale 
star catalog data. For M y = +6 the plot indicates that one can find no 
10° x 10° area on the entire celestial sphere where there is not at 

J.L 

least one 6 U order magnitude star or brighter. Guide star availabi- 
lity within the 320 square degree window is, therefore, assured. In 
fact, with 95% confidence the guide star will be a 4 tn order magnitude 
star or brighter. 

The baseline IPS star tracker, described in Appendix B, has a maxi- 
mum tracking rate of only 0.05 deg/sec. Since, as derived in Section 6, 
earth pointing will often require higher rates, payload pointing will 
usually have to be interrupted to obtain a star reading. At update time, 
the IPS will therefore interrupt payload pointing and move, if required, 
to acquire a star with the 2°xZ° F0V tracker within the available 8°x40° 
area of the celestial sphere. 
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Figure 7-5. Guide Star Availability 



In normal IPS pointing/tracking operation, the attitude control 
loops are closed through the gyros and the IPS gimbal angle resolvers 
play no part in generating attitude errors. The gimbal angle resolvers 
are used, however, in acquisition mode for attitude control, and the 
gimbal angles and direction cosine matrix Ayj are ARS outputs. Ini- 
tializing the IPS ARS requires an attitude transfer from the orbiter 
GN&C base to the IPS platform. Potentially large errors from orbiter/ 
pallet flexibility, softmount rotation, and gimbal angle resolver in- 
accuracy enter into the transfer through the relation 


Apj “ A pu A ul A lb A bi 


where A p y involves unknown gimbal angle resolver errors, Ayy involves 
unknown softmount rotation angles, A LB involves unknown orbiter/pallet 
flexibility, and Agj is the orbiter inertial attitude (relayed through 
the spacelab computer). The initial errors in the attitude transfer 
are eventually reduced by using the IPS star tracker for attitude up- 
dates. A visual star tracker, such as described in Reference 5, would - 
improve the attitude transfer and speed convergence to the required 
pointing accuracy by providing a direct visual indication of the IPS 
platform attitude. 


7.2 COMMAND AND ERROR PROCESSING 

The IPS command and error processing data flow is shown in Figure 7-6. 
The target command and orbiter ephemeris are the only inputs required 
from the spacelab computer. The target command is given in a user 
convenient coordinate frame that is identified as part of the command. 
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Figure 7-6. 










Using the orbiter ephemeris, the target command is converted to inertial 
coordinates in the command preprocessing block. This allows simple and 
uniform command processings whether the target is the earth limb, a 
point fixed to the earth surface, or an inertially fixed target, because 
the ARS provides an inertial reference. The command preprocessing also 
computes the commanded inertial rates w and angular accelerations w 
required to maintain the platform pointed at the target. The commanded 
rate is compared to the processed gyro rates to yield an inertial rate 
error o> e for the servo control law. The commanded angular acceleration 
Is used as a- feedforward input to the servo control law for improved 
performance in tracking earth-fixed targets. 

There are two modes of operation, one for acquiring targets and one 
for tracking and holding. When a new target is acquired the IPS gimbal 
angles are commanded and controlled relative to the gimbal base (upper 
part of softmount) and the resolvers are usedas position sensors. 

When the IPS is pointing at or tracking a target, the attitude error is 
computed directly from the inertial attitude reference and resolver 
inaccuracies do not influence fine pointing. 

The ultimate IPS pointing accuracy for EVAL depends not only on the 
accuracy of the ARS and IPS gimbal servos but also on the quality of avail 
able ephemeris data. It is interesting to note that ephemeris errors in- 
fluence the pointing accuracy differently, depending upon the type of 
pointing being considered. For example, space pointing of experiments 
to targets defined on the celestial sphere is independent of ephemeris. 
Furthermore, the accuracy of pointing relative to local vertical is much 
less sensitive to ephemeris errors (by the ratio of altitude to orbital 
radius) than pointing at landmark targets (i.e., latitude, longitude, 
radius). This latter point becomes apparent from the geometry shown in 
Figure 7-7, where this somewhat simplified case considers- only in-track 
ephemeris errors. Ax. To demonstrate the effects which ephemeris errors 
contribute, relative to ARS errors, it is most meaningful to consider 
angular resolution. This is summarized for a 200 Km orbit in Figure 7-8. 
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Figure 7-8. EFFECT OF EPHEMERIS ERROR ON POINTING 
ACCURACY AT 200 KM ALTITUDE 
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Note that for pointing relative to the local vertical frame the down 
track ephemeris errors given in Reference 1 add less than 2 arc sec 
to the pointing error. For landmark pointing, however, the ephemeris 
error contribute from 10 to 65 arc sec to the pointing error, depend- 
ing on the ephemeris source used and the time since the last tracking 
pass. 


7.3 IPS SERVO CONTROL LAM/SHAPING 

The IPS servo control law/shaping block diagram is shown is Figure 7-9. 
As with the command and error processing, there are two modes of opera- 
tion. In the acquisition mode, the IPS gimbal angle errors are applied 
to a proportional plus integral plus rate shaping function (implemented 
with PDE) to yield commanded gimbal drive torques.’ 

In track/hold mode, the attitude error e is processed (again 
digitally) with a proportional plus integral compensation function and 
rate feedback is provided with the computed rate error u . An important 
feature of the track/hold control law is the feedforward of commanded 
angular acceleration to improve pointing performance without excessive 
bandwidth in the feedback path. The reduced feedback bandwith minimizes 
the effect of sensor noise on attitude stability. Feedforward control 
is particularly important for raster scanning and tracking earth fixed 
targets. The commanded torque in track/hold mode is in platform fixed 
coordinates and is converted to gimbal axis coordinates j through a 
gimbal angle dependent transformation matrix. 

Selection of control law gains, especially feedforward gains and 
compensation for gimbal compliance poles, is strongly dependent on 
payload mass properties and mission requirements. This is a direct 
result of the IPS design, in particular the softmount and large center 
of mass/center of rotation offset. 
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IPS GIMBAL ANGLES FROM RESOLVERS 


NOTATION: 


t c = COMMANDED GIMBAL DRIVE TORQUES 

T = ACTUAL GIMBAL DRIVE TORQUES 

M «= TRANSFORMATION MATRIX FOR RESOLVING PLATFORM 
TORQUES INTO GIMBAL AXES 
(FUNCTION OF v» p ) 


Figure 7-9. IPS Servo Control Law Block Diagram 








8.0 EVALUATION OF ORBITER EARTH VIEWING ORIENTATIONS 

Some EVAL missions may require that the shuttle orbiter RCS be 
disabled to eliminate either the high frequency disturbances or con- 
tamination caused by firing the RCS thrusters. An earth pointing 
orientation can be maintained without the RCS by flying the orbiter 
at a null torque attitude, that is an attitude for which the total 
external torque on the orbiter vanishes. For the EVAL altitude range, 
150 to 1000 km, the predominant external torques are gravity gradient 
and aerodynamics. As the gravity gradient nulls are in general unstable 
equilibriums, the sensitivity of the torque to small perturbations 
in attitude from the null are equal in interest to the location of the 
null itself. 

The following topics relating to the existance and characteriza- 
tion of null torque attitudes and their suitability for EVAL are covered 
in this section: 

• definition of coordinate systems and transformations, 

6 definition of potential earth viewing orientations. 

0 aerodynamic and gravity gradient disturbance torque models 
for the shuttle orbiter. 

e discussion of- method of determining. null torque attitudes and 
incremental sensitivity matrices. 

.# discussion of each potential orientation, including existence 
and location of torque nulls, stability and sensitivity, and 
relative merits for EVAL. 

9 selection of two baseline earth viewing orientations, one for 
viewing the earth surface and the other for viewing the earth 
1 imb. 

e discussion of the effect of atmospheric density variations on 
the null torque attitude. 



All of the null torque attitudes found are either unstable or 
highly sensitive to atmospheric density variations at low altitude, 
and most are unstable at all altitudes. Some form of active fine 
control will therefore be required to stabilize the orbiter near the 
null torque attitude. The design of such an orbiter null torque atti- 
tude stabilization system is discussed in following sections of this 
report. 

8.1 DEFINITION OF POTENTIAL EARTH VIEWING ORIENTATIONS 

In this section potential orbiter orientations for earth viewing 
are defined. The necessary coordinate systems and transformations 
for describing and analyzing these orientations are presented first. 

8.1.1 Coordinate Systems and Transformations 

The orientation of the orbiter body axes relative to' the orbit 
reference frame is described by the 3x3 transformation matrix 


A BR a bo a on a nr 


(8-1) 


The orbit reference axes are defined in Figure 8-1 and the orbiter 
body axes are defined in Figure 8-2. 

Each nominal orbiter orientation relative to the orbit reference 
frame i-s defined by an A^ R matrix. As the nominal orientations 
involve rotations in multiples of 90 degrees about one or more prin- 
cipal axes, the A., ^ matrix contains a single + 1 or - 1 in each row 
and column and zeros elsewhere. 


The Aq^ matrix defines the rotation from the nominal orientation 
(N-frame) to the offset attitude required to achieve a torque null 
(0-frame). In general, the required offset angles will be large and 


vary with altitude. In terns of the offset angles c c 
A^ matrix is 


o’ 


the 
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Since the nominal orientations are 90 degrees apart, the maximum 
total (eigenaxis) offset rotation allowed for is 45 degrees. This 
condition can be expressed as 


Tr 



> 1 + /~2 


(8-3) 


where Tr [-] is the sum of the diagonal elements of the matrix. If 
the disturbance torque T^ cannot be made to equal (T without violating 
relation (8-3) then a null torque attitude does not exist for the • 
given nominal orientation. 

The A bq matrix describes the small' anole perturbations of the 
orbiter body from the null torque attitude, as illustrated in Figure 8-2 


4 > 

1 * 
♦ 1 

This matrix is used for stability .analysis . 
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8.1.2 Potential Earth Viewing Orientations 

. The potential earth viewing orientations considered in this 
study are listed in Table 8-1 and illustrated in Figures 8-3 to 
8-8. The orientations consist of all combinations of orbiter ro- 
tations, in multiples of 90 degrees, that satisfy both of the 
following conditions 

e some part of the earth's surface is visible from the 
payload bay. 

a the orbiter "tail" is not between tne payload bay and the 
earth. 

Each orientation is categorized by an alphanumeric "orientation 
number". 'The six possible combinations of axes perpendicular to the 
orbit plane (POP) and to nadir are denoted by the numerical part, as 
indicated in Table 8-1. The axis to nadir determines the gravity 
gradient properties. The letter suffix distinguishes the two possible 
orientations the orbiter may have with respect to its orbit velocity, 
with given axes POP and to nadir, as follows 

® A = payload bay "forward" 

9 B = payload bay "back" 
e C = orbiter nose "forward" 
o D = orbiter tail "forward" 

The axis along the velocity influences the orbiter aerodynamic 
properties. The orbiter Xg - Zg plane is a plane of symmetry, and 
therefore no distinction need be made between the + y B axes in de- 
fining the orientations. 
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Table 8-1. Definition of Potential Earth Viewing Orientations 


Orientation 

Number 

Axis POP* 

Axis to Nadir 

Axis Along Velocity 

1A 

IB 

a 

*B 

t 

+ Yg (wing "down") 

- Zg (bay "forward") 
Zg (bay "back") 

2 

- Zg (bay "down") 

+ y g (wina "forward") 

3C 

3D 

A 

* 

- z B (bay "down") 

* ^ 

Xg (nose "forv/ard") 
- Xg (tail "forv/ard") 

4A 

4B 


x B (nose "down") 

- Zg (bay. "forward") 
z R (bay "back”) 

5 

A 

Z B 

Xg (nose "down") 

+ y B (wing "forward") 

6C 

6D 

+ y B (wing "down") 

Xg (nose "forward") 
- Xg (tail "forward") 

t . 


* POP = Perpendicular to orbit plane. 
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EARTH 



Figure 8-4. Orientation 2 - XPOP, Bay Down 
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4A: 


Bay Forward 


4B: Bay Back 



Figure 8-6. Orientation 4 - YPOP, Nose Down 





Figure 8-7. Orientation 5 - ZPOP, Nose Down 
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8.2 DISTURBANCE TORQUES 


A null torque attitude is, by definition, one at which the total 
external disturbance torque is zero 

f - f A + V f S + ? M - ° (8 - 4 > 

where the four external disturbance torques acting on the shuttle orbiter 
are 

• aerodynamics (f A ) 

« gravity gradient (Tg) 

• solar (T s ) 

t residual magnetism (T^) 

The disturbance torque models adopted for this study are discussed in 
detail in Appendix D, and only the major results will be summarized 
here. 

Over the shuttle orbiter altitude range expected for EVAL missions, 

150 to 1000 km, the disturbance torques have the following relative im- 
portance 

s gravity gradient - significant at all altitudes 

« aerodynamic - largest disturbance below 200 km, negligible above 
400 km 

6 residual magnetism - somewhat uncertain, but estimated to be 
at least an order of magnitude less than gravity gradient 
and aerodynamic torques 

e solar pressure - negligible at all altitudes 
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The residual magnetism and solar pressure torques can be safely re- 
flected in this preliminary study on the ground that the mission-to-mis- 
sion variation in the gravity gradient and aerodynamic torques exceed the 
neglected terms. Orbiter fuel and payload variations result in changes 
in the orbiter inertias and CGlocation. The inertia variations, in 
particular the products of inertia, influence the gravity gradient 
torques.' while CG shifts affect the aerodynamic moment coefficients. . 

An even greater uncertainty in the aerodynamic torques is produced by 
the large (i.e., factor of 5- or more) variation in dynamic pressure at 
a given altitude due to envirionmental factors. The nominal parameter 
values used in this study are the best estimates currently available 
and are believed to be sufficiently representative to allow a valid 
comparison of the candidate null torque orientations. 

8.3 .NULL TORQUE ATTITUDES AND SENSITIVITIES 

At altitudes above 400 km, gravity gradient torque is the only 
significant external disturbance, and with the assumed zero products 
of inertia, all of the proposed nominal orientations are gravity' gra- 
dient null attitudes with zero offset angles. Table 8-2, based on 
Table D-2, summarizes the gravity gradient stability properties of the 
candidate orientations. All except orientations 4 and 5 (orbiter .nose 
"down") are unstable for small, perturbations from the null attitude, 
however, unless a control system is used. In particular, the orienta- 
tions with the most favorable fields of view towards the earth, orienta- 
tions 2 and 3, are unstable in both roll and pitch, with sensitivities 
of 1.32 N-m/rad = 0.023 N-m/deg in roll and 33.9 N-m/rad = 0.59 N-m/deg 
in pitch. Even orientations 4 and 5 may possess long- term instabilities 
when^ the complete nonlinear coupled dynamics/ kinemati cs are considered 
along with the residual aerodynamic, solar and magnetic disturbance 
torques (Reference 6). As an- example, the destabilizing effect of 
periodic atmospheric density variations on orientation 4 is discussed 
in Section 8.4. 
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Table 8-2. Gravity Gradient Stability Seminary 


Orientation Number 

Att 

itude 


Gravity gradient Stability* 

Axis- POP 

Axis 

to Nadir 

In Orbit Plane 

Out of Orbit Plane 

Both 

1 

X B 


a 

y B 

S 

U 

U 

2 • 



Z B 

U 

U 

U 

3 

y B 


Z B • 

U 

U 

U 

4 



X B 

•. s , 

S 

s 

5 

A 

Z B 


A 

X B 

s 

S 

s 

6 



y B ' 

tr 

S 

u 


* S '.= STABLE 


U = UNSTABLE 



As the orblter altitude is decreased, the aerodynamic torque 
increases rapidly, and the null torque attitudes become offset from 
the nominal orientations. Indeed, in many cases no null torque atti- 
tude exists "near 11 (i.e., within a 45 degree eigenaxis rotation) the 
nominal orientation. By definition, a null torque attitude is one for 
which - • 

T(Aq R ) = T(A qn A nr ) = 0 (8-5) 

where 

T = + T G = disturbance torque 

T a = aerodynamic disturbance torque 

Tq = gravity gradient disturbance torque 

A^ r = nominal orientation direction cosine matrix 

A()N = Aqn U 0 » 8 q , * ) = null offset direction cosine matrix 

a , e , ill = offset angles 

O 0 r 0 


Equation (8-5), when expanded with equations (8-2), (D-4), (D-10), (D-ll), 
and (D-17), is a nonlinear function of the unknowns 8 , and 
and involves the tabulated aerodynamic moment coefficients C. , C , 

Hi 

and C n . In general, numerical search techniques must be used to find 
the null attitude, if it exists, as a function of altitude. 

The sensitivity of the disturbance torque to incremental attitude 
perturbations from the null torque attitude can be expressed by the 
matrix of partial derivatives. 
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^T x / a 4 ) 

3T x /ae 

3T x /3* 

3T / 01 } 

3 T /ae 

3T /a# 

y 

y 

y 

3T 2 /3$ 

olJoQ 

3T z /3'4) 


a br a or 


The partials are evaluated at the null torque attitude. The negative 
of the above matrix has the form of a compliance or spring matrix * with 
the usual convention, that a positive spring constant K indicates a 
restoring torque 



K = 


K 


K 


K 


K 

K 

X9 

Xy 

K . 

' K . 

y> 

VP 

K , 

K 

20 

2 'P 


(8-7) 


The K matrix can be used in deriving the incremental equations of 
motion about the null torque attitude. This will be deferred until 
the next phase of the study. Instead, attention will be focused here 
on the diagonal elements of K to obtain estimates of the comparative 
sensitivity of the proposed orientations. 

Since the gravity gradient torques are analytical functions, the 
required sensitivity partials are straightforward. For example 





a 23 a 33 
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and therefore 


aT Gx 

3$ 


- 3 “o (I zz 


9a 


33 


! yy )(a 23 3* 




i 

i 




The aerodynamic torques, on the other hand, are functions of the ta- 
bulated moment coefficients C^, C m and C . For example 


T Ax = ^ ^ b c ; (“> 


and therefore 



a$ 


q S b 



da + 3C £ (c *- 6 ) 
d(p 3S 



The required partials of the moment coefficients are obtained 
by fitting the tabulated data over a rectangular region with corners 
(a , B ), (a , B i ) » (a-j , B Q ), (ci-j , ) with the interpolation function 


C(a, B) « C Q + C a (a 


a 0 ) + S (B " 6 0 ) + C a6 (a ' a 0 K ' S ' B 0 } 


(8-8) 
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where 


c o ■ C( V s o> 

C Ca-| j 3 0 ) - C (v Q , 3 q ) 
a (a } - a Q ; 


C (a-j > S-j ) + B Q ) - C(-i Q , s-j) C(a-|, “ 0 ) 

C a8 ~ _ 1 Ca 1 " 


and C(^ 0> £ q ), C(a Q , 8-j), C(c^ , 3 q ), C(a-j , Sj ) are tabulated data 
points. This interpolation function is also used in computing the 
null torque attitudes. The parti als of the moment coefficients 
are then simply 


(8-9) 


( 8 - 10 ) 

( 8 - 11 ) 


( 8 - 12 ) 


3C(o J B) 

Dot 


c + 

ct 


'aB 


(e - s 0 ) 


(8-13) 


3C(a > B) 
9B 



C «B ( “ 



(8-14) 


In the following subsections each potential orientation is 
discussed in turn. Only the more promising and generally applicable 
orientations are analyzed in depth. Conclusions are contained in 
Section 8.3.7. 
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8.3.1 Orientation 1 


Orientation 1 (Figure 8-3) has one of the orbiter wings "down" 
and the payload bay either along (1 A) or opposite (IB) the velocity 
vector. In general, no null torque attitude exists near this orien- 
tation, due to the imbalance between the aerodynamic and gravity 
gradient torques about the axis. For small offsets from the 
nominal orientation 


T y = 33.92 * 0 t 0 *qSc ( 0 Q + \ , + * Q ) 

where the top set of signs applies to 1A and the lower to IB. At an 
altitude of 200 km, for example, q S c - 30 N-m. For the range of 
offsets je Q j <_ 10 degrees, j$ | <_ 15 degrees the moment coefficient 
is in the range for orientation 1A ' 


- °- 31 f. c m < - 0.18 
and for orientation IB 

0.22 < C < 0.34 
— m — 

The magnitude of aerodynamic torque about y^ is thus at least 5.4 N-m, 
while there is no first order gravity gradient torque for cancellation.. 


8.3.2 Orientation 2 


Orientation 2 has one of the wings along the velocity vector and 
the payload bay pointing towards nadir. Again, no null torque atti- 
tude exists in general because of torque imbalances about the and 
Zg axes. The asymmetry of the orbiter with respect to the airflow 
in this orientation (mainly due to the vertical stabilizer) results 
in large aerodynamic torques around the axes (x B and z B ) that- have 
the smallest gravity gradient torques in this orientation. 
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8.3.3 Orientation 3 


Orientation 3 has the payload bay. pointing "down" and ei.ther the 
orbiter nose (3C) or tail (3D) along the velocity. Since the orbiter 
is aerodynami cal ly symmetric about the airflow in this orientation, 

= 0 and the only offset required for a torque null is about 
y g . The required offset angle e Q is plotted in Figure 8-9 for orien- 
tation 3C and in Figure 8-10 for orientation 3D. The primary abcissa 
scale is the aerodynamic constant q S c, because this is the actual 
parameter that influences the offset angle. The nominal altitude 
corresponding to the aerodynamic constant is also indicated, but 
it should be- recalled that wide variations in q at a given altitude 
occur in practice. 

Also shown on Figure 8-9 and 8-10 are the torque "compliance" 



30 


In both cases, K wC < 0, and the orbiter is unstable .in Ditch. Orien- 
tation 3C is less unstable than orientation 3D, however, at- low alti- 
tudes because the slope of the aerodynamic torque tends to oppose’ 
rather than add to the slope of the destabilizing gravity gradient 
torque. Orientation 3C is therefore preferable to orientation 3D. 
Table 8-3 summarizes the torque compliances for small perturbation ; 
from toe null torque attitude of orientation 3C. Only the yaw 
axis has a restoring torcue. Roll and pitch are unstable at all 
altitudes and a control system is therefore essential to remain near 
the null torque attitude. 
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Figure 8-9. Pitch Offset and Sensitivity for Orientation 3C 


8-23 














At t i tude 

(km) 
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300 


Torque Compliances (N-m/deq) 
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K xe 
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K 
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K z» 
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-0.023 

0 

-0.047 

o 

-0.50 

0 

-0.019 

0 

0.22 

-0.023 

0 

-0.029 

0 

-0.54 

0 

-0,0099 

0 

0.13 

-0.023 

0 

-0.0082 

0 

-0.58 

0 


0 

0.036 

'-0.023 

0 



-0.0022 

0 

. .... 

-0.59 

0 

-0.0007 

, 0 
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8.3.4 Orientation 4 


Orientation 4 has the orbiter nose "down" and the payload bay 
either along (4A) or opposed (4B) to the velocity vector. As with 
orientation 3, only a pitch offset is required to achieve. a null 
torque attitude. The required offsets are plotted in Figures 8-11 
and 8-12. It is obvious that much larger offsets are required as 
compared to orientation 3. In fact, the offset angle for orientation 
4B exceeds 45 degrees for altitudes below 180 km. This is not ne- 
cessarily a disadvantage, however, because the offset tilts the pay- 
load bay towards the earth, increasing its usable field of view for 
EVAL. The low altitude null attitudes of orientation 4A, on the other 
hand, tilt the payload bay away from the earth and are generally un- 
usable for EVAL. 


Figure 8-12 also shows that the pitch torque compliance K A is 

ye 

positive and therefore pitch is incrementally stable around the null. 
The compliance decreases rapidly at low altitude, changing sign for 
an offset angle of approximately e o ~ 60 degrees. Because of the 
variability of q with environmental factors, it would be risky to 
assume stability at altitudes below 200 km. 


Table 8-4 summarizes the torque compliances for orientation 4B. 

It is clear from the- table that the large offsets caused by the aero- 
dynamic torques at low altitudes has a destabilizing effect on the 
system and that active control is required to stay near the torque 
null. 


8.3.5 Orientation 5 

Orientation 5 nas the orbiter nose “down" and one of the wings 
in the direction of the velocity vector. There is, in general, no 
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null torque attitude near this orientation because there is no first 
order gravity gradient torque about Xg to oppose the large aerody- 
namic torque about this axis caused by. the vertical stabilizer-. 


8.3.6 Orientation 6 

This orientation has a wing pointing towards nadir and either the 
nose (6C) or tail (6D) directed along the velocity. For this orien- . 
tation there is no first order gravity gradient torque to cancel the 
aerodynamic torque about yg. The aerodynamic torque about yg is rela- 
tively small however and a torque null about, yg may exist for some com- 
bination of offset angles. It is extremely unlikely, however, that 
nulls will exist about the three axes simultaneously. Besides, the 
null, if it exists, is certain to be unstable and this orientation 
has no particular field of view advantage for EVAL. Further analysis 
of this orientation cannot therefore be justified. 

8.3.7 Cone! us ions 

In general, null torque attitudes exist at the lower altitude range 
only for those orientations having the orbiter plane of symmetry in the 
orbit plane. These are the YPOP orientations, namely orientations 3 and 
4. Of the two variations of orientation 3, the more favorable is orien- 
tation 3C with the orbiter nose "forward" and payload bay "down". This 
orientation is unstable with a negative pitch “spring constant" Ky 0 *-O.6 
N-m/deg, but provides the best view of the earth and has a small and 
relatively insensitive offset angle. For convenience, orientation 3C 
will be referred to as the "nose forward" orientation in the sequel. 

Of the two variations of orientation 4, only orientation 4B with 
the orbiter nose nominally "down" and the payload bay nominally "back", 
allows viewing of the earth limb and a significant portion of the earth's 
surface at low altitudes. The pitch offset angle for orientation 4B 
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increases rapidly with decreasing altitude, reaching 50 degrees at h 
* 180 km. However, pitch is incrementally stable and the offset increases 
the area of the earth visible from the payload bay. The roll axis is 
unstable with a maximum spring constant K x ^ “0.1 N-m/deg and again a 
fine pointing control system is required to maintain the orbiter near 
the null attitude. Orientation 4B will be referred to as the "nose 
down" orientation. ' 


8.4 EFFECT OF ATMOSPHERIC DENSITY VARIATIONS ON ALL TORQUE ATTITUDES 


Variations in atmospheric density, at a given altitude, produce 
variations in the null torque attitude through changes in the aerodynamic 
disturbance torque. The effect of these density variations on the orbiter 
motion depends strongly.on the frequency content, as can be seen with a 
linear analysis. . - 

Linearizing the orbiter pitch equation of motion about the nominal 
null torque offset angle e Q , corresponding to the nominal atmospheric 
density p Q , yields the perturbation equation 


where 


I 60 + K 50 = 6T 

y ye 


Ap 


(8-15) 


y 

66 




orbiter pitch inertia 


attitude perturbation from nominal null torque 
offset angle 6 Q 

* nominal disturbance torque 
9 0 - e o , p' = p o compliance 


_ 

. ee 


8T A y 
3p ' 


(p - p 0 ) = aerodynamic disturbance 
e , p = p 0 torque perturbation due 

0 0 to atmospheric density 

variation 
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For a sinusoidal atmospheric density variation of magnitude &p at 
frequency u> p , 

p - p n = Ap sin to t 
o p 

the magnitude of the orbiter response is for K > 0 

ye 


and for K,. < 0 

yu 


66 


1 

1 


6T Ao 

K ye 

l - 

K /to e )2 


60 


1 



1 

_1 

6T Ap 


K ye 


1 + («i p /u e ) 2 


where the orbiter incremental pitch motion natural frequency is 


“9 ’ (|K ye |/I y ) 1/2 

The normalized responses, |Ky 0 6e/ST^|, are plotted in Figure 8-13 
versus the normalized frequency ratio oj /w c for the two cases. 

p 0 

For frequency ratios far removed from unity, the response is in- 
dependent of the sign of K _. When the atmospheric density variation 

yu 

is at a low frequency compared to a) Q , the orbiter tends to follow the 
resulting motion of the null torque attitude. In contrast, the orbiter 
shows little response to density variations with frequencies that are 
high compared to co,. For density variations with frequencies near ui_, 

o O 

the response depends drastically on the sign of Ky 0 . If K y0 >0, as 
in the nose down orientation, the response shows a sharp resonant peak 
around . Thus, although the nose down orientation is stable in 


(8-16) 


(3-17) 


(8-18) 


(8-19) 
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pitch in a constant atmospheric density environment, atmospheric density 
variations at w p = oj q can have- a strong destabilizing effect by producing 
orbiter attitude motions far larger than the change in null torque atti- 
tude. On the other hand, if K. < 0, as in the nose forward orientation, 
pitch is unstable regardless of the density variations and density varia- 
tions at frequencies near m have no particular additional effect on the 
system stability properties. 

The effect of atmospheric density variations on the orbiter null 
torque attitude is discussed in quantitative terms in the following sub- 
sections. 

8.4.1 Long Term Variations 

Long term atmospheric density variations are those with frequency 

components u p < < oj 0 . Over the altitude range 180 <_ h £ 300 km. Tables 

8-3 and 8-4 show a range of 0.25 ] K^ 0 1 _< 0.59 N-m/deg or 14 £ |K Q | 

< 34 N-m/rad. With the nominal value 1^ = 9.39 x 10 6 Kg - m 2 . Equation 

(8-19} yields the range of orbiter incremental pitch motion natural 

-3 -3 

frequencies 1.2 x 10 £n) 0 <_1.9 x 10 rad/sec. The corresponding 

periods for one cycle of pitch motion are 3300 < T < 5150 sec, or on 

0 

the order of one half to one orbit period. Therefore any density 
variation with a period greater than roughly four orbits or about six 
hours can be considered long term. 

The pitch null torque offset angle is plotted as a function of 
altitude and percent deviation from nominal atmospheric density in Fig- 
ure 8-14 for the nose forward orientation and in Figure 8-15 for the 
nose down orientation. These results apply directly in a constant 
atmospheric density environment or when the density varies with a period 
of greater than about six hours. It is apparent from these figures that 
both the nominal offset angles and the variations, at a given altitude, 
are more than an order of magnitude larger in the "nose down" orientation 
than in the "nose forward" orientation. This is due to the larger aero- 
dynamic torque in the "nose down" orientation, resulting from the larger 
surface area intercepting the air stream. Figure 8-16 compares the nominal 
aerodynamic torques in- the two orientations. 
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Figure 8-14. Variation of Pitch Null Torque Offset Angle 
With Atmospheric Density In "Nose Forward" 
Orientation 
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Figure 8-1 5, Variation of Pitch Null Torque Offset Angle 
With Atmospheric Density In "Nose Down" 
Orientation 
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Figure 8-16. Pitch Aerodynamic Torque at Nominal 
Null Torque Offset Angle 
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8-4.2 Short Term Variations 


The orbiter dynamics attenuates the orbiter attitude response to, 
null torque attitude variations due to atmospheric density variations 
occurring at frequecies w >> to.. In practical terms, this means that 

P O 

the orbiter will not respond significantly to disturbance variations with 

a) >3 to., that is, with periods of less than about twenty minutes. 

P o 

In the nose down orientation, the linear analysis indicated a 

resonant type of response for w. =co Q . As explained in Appendix D, a 

P o 

significant atmospheric density variation can occur at orbit frequency 
due to differences in density on the day/night sides of the earth. As 
mentioned previously, w Q is also nearly equal to the orbit frequency 

u 

for altitudes below 300 km. This unfortunate coincidence leads to the 
conclusion that atmospheric density variations may be highly disruptive 
to the otherwise stable nose down orientation at low altitudes. 

As the linear analysis is strictly valid only for small perturba- 
tions, simulation is required to obtain meaningful results in the near 
resonant case. A digital simulation was therefore written to simulate 
the orbiter rigid body pitch dynamics along with the complete disturbance 
torque model described in Appendix D. The case run simulates the orbiter 
in the nose down orientation at 200 km altitude with a +20% sinusoidal 
variation in atmospheric density at orbit frequency. Time plots of 
the results, spanning slightly more than two orbits, are presented in 
Figures 8-17 and 8-18. The variable e g = 9 0 + 0 in Figure 8-17 repre- 
sents the total pitch angle from the nominal nose down attitude. For 
this case, 6g is initially equal to the nominal offset angle for the 
nominal atmospheric density at this altitude. In other words, the orbiter 
starts at the nominal null torque attitude with the nominal atmospheric 
density. The atmospheric density variation induces a divergent pitch 
motion, reaching a peak-to-peak amplitude of over 57 degrees. In com- 
parison, Figure 8-15 shows that the null torque attitude has a peak-to- 
peak variation of only 10 degrees under the stated conditions. Figure 8-17 
also shows that the pitch rate and total disturbance torque are also diver- 
gent. 
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TIME (SEC) 

Figure 8-18. Disturbance Torques on Orbiter With +20% Atmospheric 
Density Variation at 200 Km in Nose Down Orientation 








Figure 8-18 expands the total disturbance torque into its gravity 
gradient and aerodynamic components. The gravity gradient and aerodynamic 
torques are initially equal in magnitude and opposite in sign, signify- 
ing that the orbiter is initially at the null torque attitude. The aero- 
dynamic torque variation is nearly sinusoidal at orbit frequency, at 
least initially, -reflecting the fact that the aerodynamic torque is 
proportional to atmospheric density but only weakly dependent on attitude. 
The gravity gradient torque is a nonlinear function of attitude only, 
and shows a complicated variation with large amplitude. 
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9.0 ORB ITER CONTROL REQUIREMENT ANALYSIS 

* 

The analysis in the preceding section has demonstrated that the 
external disturbance torques cause all potential low altitude (h< 400 
Km), and most of the high altitude, orbiter null torque attitudes suit- 
able for earth surface or limb viewing to be unstable. An active atti- 
tude stabilization system is therefore required to keep the orbiter near 
the null torque attitude. 

At a minimum, the oroiter null torque attitude control system should - 
meec the following requi reme'nts : 

e incremental stability - orbiter stays near torque null if 
started nearby. 

e ability to capture from orbiter vernier RCS limit cycle rates 
and attitudes ’n a reasonable time. 

o relative insensitivity to such uncertain Quantities as the 
location of the torque null and the slopes of the torques 
about the null (K matrix) 

o capable of controlling the orbiter in either the nose toward (3C) 
or nose down (4B) orientation, with at most a change of gains. 

These requirements are examined in general terms in this section, 
without regard to the type of control actuator used. The particular impli- 
cations of using either magnets cr momentum exchange devices for control 
are discussed in detail in following sections. 


9.1 STABILITY ANALYSIS 


The incremental staoility of the orbiter about the null torque at- 
titude can be determineo from the eigenvalues of the matrix A in the 
state equation 

x = A x 


(9-1) 
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The state vector is 


*' 

0 


X = 


. Tp 


to 

x 




( 9 - 2 ) 


where <j>, 6, ij> are the small angle attitude perturbations from the torque 
null, -o , <o , to are the body components of the inertial rates, and to 

X jr Z 0 

is the orbit rate. It is important to realize that the location of the 
torque null is uncertain and variable; the stability analysis is, how- 
ever, conducted about the true null. The 6x6 matrix A can be expressed 
as the sum of three matrices 


A " Dynamics + Disturbance + Dontrol 


( 9 - 3 ) 


For either the nose forward or nose down orientation, the dynamics/kine- 
matics are described by 
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The slopes of the disturbance torques about the null enter in the 
partitioned matrix 
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0 

0 

A 



^DISTURBANCE = 




- I ‘ K 

0 


— 



(9-5) 


where 0_ is the 3x3 zero matrix, I is the diagonal inertia matrix for 
the orbiter 


xx 
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and K is the disturbance compliance matrix defined in Section 8.3 
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The nominal values of K for orientations 3C and 4B are given in Section 8.3 
for several alti.rudes*;the units must be converted to il-m/rad for use here. 

The control law is assumed for now to be state variable feedback, 
leading to 

A C0NTR0L = 


I -1 G 


(9-8) 


where G is the 3 

x 6 feedback control 

gain 

matrix 
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The control torque vector T c is 


T c = 


T 


Cx 


Cy 


Cz 


= G x 


(9-10)' 


9.1.1 Open Loop Stability 

With the control loops open, stability is determined by the eigen- 
values of 


Dynamics + Disturbance 


Evaluating these eigenvalues for the nose forward and noise down orientations 
over the altitude range 180 to. 300 km reveals the following: 

e roll - unstable in both orientations at all altitudes 

e pitch - unstable in the nose forward orientation at all 
altitudes, pure imaginary eigenvalues (oscillatory with 
no damping) in the nose down orientation at all altitudes. 

$ yaw - marginally stable (oscillatory with damping ratio 
? ^O.Q2)in both orientations at all altitudes 

Both orientations at all altitudes thus have at least one unstable 
axis, and long term earth viewing with the RCS disabled will require 
active control. 
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9.1.2 Closed Loop Stability 


In general. Equation (9-1) represents a coupled sixth order system. 
In the orientations under consideration, however, near the null torque 
attitude 


= K. = K... = K = 0 


X9 


yo 


y-p 


Z6 


Therefore pitch is uncoupled from roll/yaw and there is no reason not 
to choose 


G 2x ‘ G 5x = G ly G 3y ~ G 4y G 6y " G 2z ~ G 5z 0 


The remaining ten gains are available for achieving the following 
desired dynamic characteristics 

© stability - all eigenvalues in left half plane. 

o decoupling - each axis to behave as an (approximately) uncoupled 
second order system. 

o insensitivity. - both of the above conditions to be met in both 
orientations at all altitudes. 


Greater insight into achieving these goals can be gained by using 
the relations 
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* * “x ' “d * 

(9-11) 


9 = » y - „o 

(9-12) 


<l> ~ + ^ 

(9-13) 

to change variables. The dynamic equations become 



Z xx 4 + B Rx * + B Px * = C Rx * + C Px 0 

(9-14) 


! yy 6 + B Ry ® + B Py 6=0 

(9-15) 


X zz * + B Rz * + B Pz * = C Rz * + C Pz 4 

(9-16) 

where the 

effective rate and position gains are 


B Rx 

= ' G 4x 

(9-17) 

B Px 

- - G lx + “o V “o <‘zz - 'yy> + K Xj , 

(9-18) 

CO 

<< 

= ' G 5y 

(9-19) 

B Py 

" ‘ s 2y + V 

(9-20) 

8 Rz 

= ‘ G 6z 

(9-21) 

B Pz 

* ' G 3z ‘ "o a 4z + “o {I yy " ! xx> + K zj. 

(9-22) 
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and the interaxis coupling gains are 


C Rx ' G 6x ' “o 
C Px = S 3x + “o 

C Rz = G 4 z + “o 
C Pz = G lz " “o 


f r zz + "xx 


V 


G 4x ' K x„ 


("zz + "xx 


V 


G, - K 
6z z<p 


(9-23) 

(9-24) 

(9-25) 

(9-26) 


The design goals can be met, by applying classical control 
theory, in the following ways: 

o stability: select the control gains, G, to yield positive 

"rate gains" 3 R . and "position gains", Bp.. , i = x, y, z. 

a decoupling: select the control gains, G, to yield (nearly) 
zero "rate coupling gains" and "position coupling gains" 

Cp.j , i — x , z . 

a insensitivity to disturbance variations: select the control 
gains, G, to insure the previous two conditions are met over 
the expected range of disturbances, K. 

If adequate control torque is available and the state variables 
can be measured, the control gains can be selected to yield a closed 
loop system with virtually any desired response. With the limited con- 
trol torque available from momentum exchange devices, and more particu- 
larly magnets, tight control of the orbiter around the null torque atti- 
tude is not feasible. A more realistic goal is to merely stabilize the 
orbiter, with relatively low gain loops, yielding low bandwidth, low damp- 
ing responses. For example, the orbiter can be stabilized about the null 
torque attitude in a constant atmospheric density environment, in either 
the nose forward or nose down orientation, at all altitudes above 180 km 
with the gain matrix 
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The variation of the closed loop roots over the altitude range 180 to 
300 km, using the associated nominal disturbance compliances from Sec- 
tion 8.3 is shown in Figure 9-1. The relatively large difference in 
pitch roots between the two orientation can be reduced by using a dif- 
ferent value of ^ or eac * 1 orientation. 

9.2 ACTUATOR REQUIREMENT ESTIMATE 

In this section, estimates of the control actuator requirements are 
obtained. Torque requirements can be better defined without extensive 
analysis than can momentum requirements, and apply to both magnets and 
momentum exchange devices. Therefore emphasis is placed on torque re- 
quirements although momentum storage requirements are addressed where 
appropriate. Initial capture fronvthe orbiter RCS limit cycle is con- 
sidered first, followed by an analysis of normal operations, and finally 
concl usions . 

9.2.1 Capture From RCS Limit Cycle 

An approximate analysis of actuator torque requirements for initial 
capture from the orbiter vernier RCS limit cycle is presented here. 
Momentum requirements for momentum exchange actuators are considered in 
Sections 11.3 and 11.4. 
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For the second order system 


with feedback control law 

In + Kn * T 

n 

(9-27) 

*r 

r 

n 

" % " + G P„ n 

(9-28) 

and initial conditions 

O 

II 

o 

(9-29) 


n(0) = n Q 

(9-30) 


the peak torque required for control is approximately, for small values 
of the damping ratio r. 


IT 


n ‘peak 



I (1 + 2 d - Kl 
0 - C 2 ) 1/2 


2 + , n 0 2 

"o 11 ' 
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(9-31 ) 


where the feedback gains are related to the bandwidth a> and damping 
ratio r by 


G r = - 2 03 5 I (9-32) 

G p = - (ui 2 I - K) (9-33) 


The estimate in equation (9-31) assumes the worst case combination 
of the signs of the initial conditions and is exact in the limit as 
^ •+ 0 . 



For the orbiter, the largest torque requirement is likely to be 
along pitch in the nose forward orientation. For this case 


K = K.„ = - 33.8 N-m/rad 

ye 

I = I, = 9.39 x 10 6 kg-m 2 

yy 


and the worst case initial conditions are the estimated attitude error 
and rate at the time the RCS thrusters are disabled 


n Q = 6 (0) = 0.1 deg = 1.745 x 10 2 rad 
n Q =6(0) = 0.01 deg/sec = 1.745 x 10 ^ rad/sec 


With these parameters, the peak torque in equation (9-31) is minimized 



for t ~ 0, 03~2 x 10 rad/sec, resulting in 


IT 


Cy 1 Peak 


« 6.23 N-m 


1 T i 

Minimizing 1 Cy ' p ^ leads to completely unsatisfactory performance-* 
however, with a continuous pitch oscillation with amplitude 




a 


8.73 x 10'" rad = 5.00 deg 
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It can therefore be concluded that a control torque greater than 6.23 
N-m will be required for satisfactory attitude capture from the .RCS 
limit cycle with the 'assumed error and rate at the transition time. - 

An order of magnitude reduction in orbiter RCS limit cycle rate 
may be possible through modification of the RCS. This would reduce the 
rate initial condition to 

n Q = e(0) = 0.001 deg/sec = 1.745xl0~ 5 rad/sec 
The peak torque in Equation. (9-31) is minimized for ? = 0, 0.0018 

_3 

xlO rad/sec, resulting in 


Vpeak “ 0- 633 N ’ m 


9-2. 


For future reference, l^Qylpg^ is plotted against e(Q) in Figure 


9.2.2 Normal Operation 

Following, the initial capture transient, the. period o'f normal opera- 
tion is entered. If the null torque attitude was constant and known 
exactly, essentially no control torque would be-irequired to remain at . 
the torque null once the initial transient was damped out. Realistically, 
however, the null torque attitude is neither constant over the period ■■ 
the orbiter is to be stabilized nor readily determinable in real time. 

The largest source of variation and uncertainty in the location of 
the null torque attitude at low altitudes is temporal variations in at- 
mospheric density. In the two selected YPOP orbiter orientations, the 
gravity gradient and aerodynamic components of the total disturbance torque 
individually vanish at the null torque attitude, for both the orbiter 
roll and yaw axes. Therefore atmospheric density variations have no signi- 
ficant effect on the roll and yaw axes, and actuator sizing is based on 
the initial capture transient. Along pitch, however, both the gravity 
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• Figure 9-2. Minimum Peak Torque Required for Pitch Capture From RCS Limit Cycle 
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gradient and aerodynamic components of the total disturbance torque may 
have large magnitudes (Figure 8-16) at the null torque attitude, with 
cancellation resulting from opposite signs. Sensitivity to atmospheric 
density variations can therefore be of sufficient magnitude to influence 
actuator sizing. 

A linearized incremental model of the pitch control system is illus- 
trated in Figure 9-3. The model is linearized about the null torque off- 

i 

set angle corresponding to the nominal aerodynamic disturbance torque. 

The dynamics are described in general by the pair of equations 


where 


ly 80 = sT c + 5T e + 5T Ap 





60 


Sh = change in momentum exchange device stored momentum 
6T c = incremental control torque 


and the remaining quantities are as defined in Section 8.4. Equation 
(9-35) applies only when a momentum exchange device is used as the con- 
trol actuator. 


Assume that the aerodynamic torque variation is at frequency co with 
amplitude 6T^ p (u), that is 


6T^ = (co) sin cot 


Two idealized control strategies will be considered: 

(1) Follow the null torque attitude exactly at all times 

(2) Hold the nominal null torque attitude exactly at all times 


(9-34) 

(9-35) 

(9-36) 


(9-37) 
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MOMENTUM EXCHANGE 
DEVICE DYNAMICS 


Figure 9- 



Incremental Pitch Model for Analysis of Response to Aerodynamic 

Disturbance Variations 









With the first stragery, the total external disturbance torque pertur- 
bation vanishes, that is 


ST e + 6T Ap = 0 

This imp! ies 


I y 69 - 6T C 


6T Ao 

60 = — sin toT 

K ye 

Since 66 is given explicitly by Equation (9-40), derivatives can be 
taken to yiel d 

6T = - 6h = - I. to — ^ sin (oT 


6T„ (w) 

6h = -I,, to — rr- COS tot 

y ye 


Therefore the ratios of peak torque and momentum variation to aerodynamic 
disturbance variation are 


G l<“> - 6T^W 


T 2 
i. . to 


H-j (w) = 


(9-44) 
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With the second strategy, the attitude perturbation is zero 


66 = 0 (9-45) 

resulting in 

6T C = - ST Ap (9-46) 

Therefore the ratios are in this case 


(9-47) 


(9-48) 


Comparison of equations (9-43) and (9-44) with Equations (9-47) 

and (9-4 8) shows that the first strategy requires smaller actuators 

if co < yl Ky 0 !/Iy while the second strategy requires the smaller 

actuators if u> > 1~* For EVAL, the predominant component of 6T Ap 

is at a = 1.1 x 10' 3 rad/sec and I = 9.39 x 10® Kg-m 2 . The value 

of | K ] decreases with increasing offset angles but using the maximum 

value [K 1= 33.92 N-m/rad yields 

ys 
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H-j (w 0 ) = 304 N-m-sec/N-m 


6 2 (w 0 ) = 1 N-m/N-m 
H 2 {u 0 ) = 909 N-m-sec/N-m 

Up to three times as much momentum and torque is thus required to 
hold a constant attitude compared to following the null torque at- 
titude motion at orbit frequency. 

Fiqure 9-4, based on Figure 8-16, is a Mot of 6T c = G 2 (aj Q ) 6T Ap 5 
or the pitch control torque required to hold the orbiter at the nominal 
null torque attitude in the face of a +20% variation in atmospheric 
density at orbit frequency. The torque values are in all cases less 
than that required to capture from the nominal 0.01 deg/sec RCS limit 
cycle (6.23 N-m). However, in the nose down orientation below 235 km, 
the indicated control torques are greater than that required to capture 
from the reduced rate (0.001 deg/sec) RCS limit cycle (0.633 N-m). 

Following the null torque attitude may lead to large orbiter 
motions at low altitude, and the linear analysis is not strictly valid 
for determining the control torque requirements. Under most conditions, 
however, the control torque requirement for following the null torque 
attitude is approximately one-third of that required to hold a constant 
attitude at the same altitude. 

With the current nominal orbiter parameters, the largest torque 
requirement is for capture from the orbiter RCS limit cycle, rather 
than during normal steady state operation. The steady state torque 
requirement is sized largely by the variation and uncertainty in the 
location of the null torque attitude, as determined by the atmospheric 
denisty (the primary influence at low altitudes) and the orbiter mass 
properties. The dominant influence on the transient torque require- 
ment for attitude capture from the RCS limit cycle is the nominal 
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Figure 9-4. Pitch Control Torque Required to Hold at Nominal Null 
Torque Attitude With 20% Atmospheric Density Variation 
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limit cycle rate of 0.01 deg/sec rather than the nominal limit cycle 
dead 2 one of 0,1 deg. The' size of the actuators required, for the null 
torque attitude stabilization system can be reduced significantly if 
an additional "Tow rate" RCS mode is made available. It is recommended 
that such a mode be provided, with the limit cycle rates reduced by 
at least a factor of 10 from the current nominal. A proportional in- 
crease in the RCS deadzone, if required to achieve the reduced rates, 
would be acceptable. 

.It .should be emphasized that idealized control has been assumed. 
Means of implementing practical approximations to these idealized con- 
trol strategies are described in the following sections. 


9-21 



10.0 MAGNETIC CONTROL OF ORBITER 


If electromagnets are used for control , the magnetic torque is 
related to the magnetic moment of the magnet bars II and the earth's 
magnetic field B 7 by 


T., = M x B 


( 10 - 1 ) 


From the properties of the cross product, it is clear that no torque 
can be generated in the direction of B and in general there will be 
a torque error. T g between T^ and the desired control torque That 
is. 




( 10 - 2 ) 


_ I/O 

It can be shown that the magnitude of the error torque, (T * T ) ' 
is minimized if the magnetic moment is chosen to be 


B x T r 

M = ~ — zr (10-3) 

B • B 
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The minimum magnitude error torque is along B and proportional to 
the component of T c along B, that is 


(T r * B) 

(T ) . = £ 

e min B « B 


B 


The relation between these vectors is illustrated in Figure 10-1. 

\ 

In a practical implementation, two additional complication 
enter. The ambient magnetic field must be measured with a three- 
axis magnetometer to implement equation (10-3), leading to the mag- 
netic moment command 



B x T c 

A A 

B • B 


where B is the measured field. In addition, a given magnet can only 
produce a limited magnetic moment. A simple but realistic model for 
the magnets. is therefore 


M = SAT (M c , m) 


where 

M = actual magnetic moment 

SAT(», •) = vector saturation function 

m = saturation magnetic moment 


(10-4) 


(10-5) 


( 10 - 6 ) 
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M = MAGNETIC MOMENT 



B = MAGNETIC 
FIELD 


B x T. 
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Figure 10-1. VECTOR RELATIONSHIPS FOR MAGNETIC 
' CONTROL SYSTEM 
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A magnet control system incorporating these" features is illustrated 
in Figure 10-2. 

10.1 CONTROL DESIGN 

The error torques in Equation (10-4) introduce undesirable time- 
varying coupling between the control axes that may intefere with con- 
trol performance. Since control torque is not available in. the di- 
rection of B, three-axis stability cannot be assured, especially if 
the open loop system is unstable. The "most unstable" axis for the 
null torque orientations under consideration is pitch in the nose 
forward orientation and therefore attention will be focused on pitch. 

Figure 10-3 shows a pitch control system using magnets. The dis- 
turbance torque model explicitly includes the offset of the null 
torque attitude from the nominal attitude, e , and the attitude error 
from the torque null e. The total pitch angle from the nominal at- 
titude (i.e., nose along orbit velocity in nose forward orientation, nose 
towards nadir in nose down orientation is 


e B = e o + e 


(10-7) 


while the rate variable is 


Sn = 


u. 


(10-8) 
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NOTATION: 

T = Command Torque 

T^ = Actual Magnetic Torque on Orbiter 
13 - Actual Magnetic Field 

B = Measured Magnetic Field 
M = Commanded Magnetic Moment 
M = Actual Magnetic Moment 

= Saturation Magnetic Moment, i = x, y, z 


Figure 10-2. Magnet Control System 
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Figure 10-3. Pitch 
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System Using Magnets 









A command input @ c is shown entering the attitude control block of 
Figure 10-3. For present purposes it is sufficient to assume that this 
input will be used to implement one of the idealized control strategies 
presented in Section 9.2.2. That is, e c is used as follows in each of 
the control strategies 

(1) e = actual null torque attitude at all times, to 

follow the null torque attitude 

(2) e = nominal null torque attitude at all times, to 

c hold a fixed attitude 

Practical means of approximating these strategies will be discussed in 
a later section. 


10.2 IMPLICATIONS OF MAGNETIC CONTROL 

In the previous section it was assumed that the required pitch 
control torque could be produced by the magnets. This assumption will 
now be examined, along with the inter-axis coupling effects introduced 
by the magnets. 


10.2.1 Interaxis Coupling 

Commanding a torque about one body axis with a magnetic control 
system will in general produce magnetic torques in all three axes. 

In particular, if the commanded torque is about pitch only 
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Equations (10-2) and (10-4) give the magnetic torque as 



( T C • I ) B 
B • B 


( 10 - 10 ) 


or in this case 


b x b y / b 2 


T M =I 


0 - b 2 )/b 2 


T 


Cy 


b y V b 


( 10 - 11 ) 


where 





(10-12) 


It is clear that in this case if and only if b = 0. 

Using the tilted dipole magnetic field model described in Appendix E, 
Figures 10-4 to 10-6 contain plots of the normalized magnetic torques 
V T Cy for orbit inclinations of 30 s 60 and 90 degrees. These plots 
are representative over the entire 150 to 1000 km EVAL altitude range. 
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It is obvious from these figures that in low inclination orbits i.t is 
impossible to generate a magnetic torque solely, or even predominantly, 
along pitch. As has been previously stated, the control law embodied 
in Equation (10-10) is optimal with respect to minimizing the error 
between and and therefore no improvement is possible. 

Figure 10-7 summarizes the ranges of the normalized magnetic 
torques as a function of orbit inclination. The maximum and minimum, 
values are generally reached at least once each day. Therefore even 
at the "best" orbit inclination, i = 90 degrees, the coupling between 
the pitch torque command and roll magnetic torque will reach 


T Mx /T Cy 


0.19 


and the coupling between pitch torque command and yaw magnetic torque 
command and yaw magnetic torque will reach 


T Mz /T Cy 


0.11 


part of each day. Since pitch is unstable' (at all altitudes) in the nose 
forward orientation, pitch control torque must be applied continuously 
and not just when the magnetic field conditions are favorable. In- 
terference with ‘roll and yaw control is therefore inevitable and, 
considering that roll is also open loop unstable, three-axis stability 
cannot be guaranteed. 
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10.2.2 Magnet Sizing 


To successfully control an unstable axis, the magnets must be 
sized to produce sufficient control torque at all points of the 
orbit. From Equation (10-3), the magnetic moments implied by a pitch 
torque command are 


M x » - (b z /b 2 )T Cy 

(10-13) 

v° 

(10-14) 

M z = (b x /b 2 )T Cy 

(10-15) 


The normalized magnetic moments M x / T Cy anci M z /T^ are plotted in Figures . 
10-8 to 10-10 for orbit inclinations of 30, 60 and 90 degrees at 200 km 
altitude. The peak values from these, and other similar plots, are sum- 
marized in Figure 10-11. It is apparent that the normalized magnetic 
moments range from about 2 x TO' 7 to 5 x 10 7 pole-cm/N-m over the EVAL 
.altitude range. Recall (figure lOr-7), however, that only near 90 degree • 
inclination does the actual magnetic torque along y, T^, come close 
to equalling the command torque Tj. . When. this fact is taken into 
consideration, the magnetic moment required to produce a given actual 
magnetic torque along pitch is minimized by assuming a 90 degree in- 
clination orbit. This highly optimistic .assumption is made in the 
following analysis. 

4 

Magnets with magnetic moments in the 10 pole-cm range have been 
used for momentum management on such spacecraft as TIROS, 0S0, 0A0, 


1 0-14 
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Figure 10-11. Peak Normalized Magnetic Moments 



6 

and SAS. Larger magnets, in the 10 pole-cm range, have been proposed 
for LST and other large spacecraft. It will be assumed here that 
the current practical limit is 10 7 pole-cm per axis. The estimated 
physical properties of such a set of magnets are listed in Table 
10 - 1 . - 


Another important consideration, that will not be addressed 
here, is the effect such large magnets will have on the operation of 
nearby payloads. 

The maximum RCS limi't cycle rates in pitch, from which attitude 
capture is possible with the 10 7 pole-cm magnets, are listed as a 
function of altitude Tn Table 10-2. The table is based on Figures 9-2 
and 10-11 with a 90 degree inclination orbit. The table shows that 
capture from the nominal 0.01 deg/sec limit cycle is not possible, even 
with the optimistic assumptions made. In fact, the limit cycle rates 
must be reduced by a factor ranging from 23 at 200 km, to 40 at 1000 km. 
Even if a "low rate" RCS mode were provided, the variability and un- 
certainty inherent in the RCS would almost certainly preclude reliably 
achieving the rates required. 

Assuming that the limit cycle rates can, somehow, be reduced to the 
required levels for attitude capture, there remains the need to supply 
continuous control torques to either follow the atmospheric density varia- 
tion induced changes in the null torque attitude or to hold a constant 
atttitude. Figure 10-12 compares the pitch torque required to hold a con- 
stant attitude with the torque available using 10 7 pole-cm magnets. It 
is clear from the figure that even 10 7 pole-cm magnets are incapable of 
stabilizing the orbiter against short term atmospheric density variations 
in the nose down orientation below 265 km altitude. 

10.3 CONCLUSIONS 

Even if several optimistic assumptions are made, magnets are 
incapable of providing prolonged three axis stabilization for the orbiter. 
Among the points weighing against magnets are the following 
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Table 10-1. Estimated Physical Properties of 10 7 
Pole-cm Per Axis Magnets 


Parameter 

Per Axis 

Total 3-Axis 

Magnetic Moment 
(Pole-cm) 

!0 7 


Mass (kg) 

146 

438 

Length (m) 

7.70 

- 

Diameter (m) 

0.065 

- 

Volume (m^) 

0.026 

0.077 

Power (Watts) 

43.3 

130 




• 
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Table 10-2. Maximum RCS Limit Cycle Rate for Pitch Attitude Capture With 10 7 Pole-cm 

per Axis Magnets in 90 Degree inclination Orbit 


Altitude 

(Km) 

Normalized Magnetic Moment 
(Pole-cm/N-m) 

. 

Pitch Torque With' 

7 •' 1 

m =m =1Q- , Pole-cm 

, A 4. 

Maximum RCS Limit 
Cycle Rate for At- 
titude Capture 

Ratio of Nominal 
to maximum RCS Li- 
mit Cycle Rate 

; 

V T Cy 

V T Cy 

■ T Cy (N-m) 

W ( de 0/sec) 

• . • 

0 Nonr e Max 

180 

2.01' x, 10 7 

- 

3.49 x 10 7 

0.287 

4.3 x 10" 4 

i • 

23 

200 

2.03 x 10 7 

3.51 x 10 7 - • • 

0.285 

4.3 x TO" 4 

23 

•250 

2.08 x 10 7 

3.60 x TO' 7 . 

0.278 ’ 

• 4.0 x 10" 4 

25 . 

300 

2.13 x 10 7 

*. . 3.69 x iq 7 ;. “ ; 

0.271 

. 3.9.. x 10“ 4 

26 

500 

2.32 x 10 7 

4.02 x 10 7 ‘ : 

' ’ 0.249 

3.5 x 10~ 4 '" 

29 

1000 

2.87 x 10 7 

4.97 x 10 7 

0.201 

2.5 x 10" 4 ‘ 

40 


* 


» 

e 


Nom 


= 0.01 deg/sec 
















Torque (N-m) 


Torque available with 10 
pole-cm magnets in 90 degree 
inclination orbit 



Altitude (Km) 

Figure 10-12. Comparison of Required Control Torque 
With Available Magnetic Torque 


9 Physical limitations of magnetic control due to the relative 
orbi ter/magnetic field geometry introduce undesirable inter- 
axis coupling, with adverse effects on stability. 

a very large magnets (10^ pole-on per axis) required for even 
marginal performance. 

• RCS limit cycle rates must be reduced from current' nominal values 
by a factor of 25 to 40 to permit attitude capture. The feasibi- 
lity of doing this is highly doubtful.- 

e -lack of reserve torque capability can result in’ loss of control 

■' due to- variations in atmospheric density or orbi ter mass pro- 
perti es . 

a possible adverse effects of large magnetic field on nearby payloads 
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11.0 MOMENTUM EXCHANGE CONTROL OF ORBITER 


Compared to magnets, momentum exchange devices offer the possibility 
of more accurate and sophisticated control, due to the higher magnitude 
and arbitrary direction of the available control torque. Potentially ap- 
plicable momentum exchange devices- include single gimbal and double gimbal 
CMG's (control moment gyros) and RW's (reaction wheels). For brevity, 
the discussion that follows is stated in terms of CMG's; the analysis ap- 
plies to RW's as well, however. 


11.1 SYSTEM EQUATIONS 

Rather than assuming a particular CMG configuration, the discussion 
will be in terms of the orb iter body axis components of the CMG cluster 
momentum 



The state vector x will be defined as the 12- vector 


6 



where. h was defined above and 



6 


(11-3) 


* 


* 


= small angle perturbations from 
null torque attitude. 


= to - 


“o y B 


U> "CO 

1 y o 


to. 


inertial rates of body relative to 
orbit reference frame. 


(11-4) 


n = 


= time integrals of wheel momenta 


(11-5) 


The reason for .including n will become apparent later. 

With these definitions, the small angle dynamic and kinematic 
equations become 


- _ r 1 

“R • 1 


[to x (I * to + h) + Ks + 


" + ‘V 


( 11 - 6 ) 


5 = tOj^ + to <5 


(11-7) 
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where I is the orbiter inertia dyadic, K is the disturbance torque com- 
pliance matrix, h is the CMG torque vector, ST^ is the aerodynamic dis- 
turbance perturbation due to atmospheric density variations and 


0 


0 



(U = 


0 0 


O) 


0 


0 


0 


( 11 - 8 ) 


The differential equation for n is, by assumption. 


n 


= h 


(11-9) 


There remains to determine the differential equation for h, defining 
-the CMG control law. The form of this equation can be selected to yield 
the desired control characteristics, and will be assumed here to be a 
combination of linear state variable feedback and a term to cancel the 
nonlinear coupling terms in the dynamics 


h s S x - ux (I • #+ h) 


(11-10) 


where G is in general a 3 x 12 gain matrix composed of four 3x3 
submatrices 


11-3 



( 11 - 11 ) 



Since the exact location of the null torque attitude is not known 
a priori, a more realistic form of the control Taw is 


h = Gx-cdx{I • to + h) (11-12) 


where x is obtained from x by substituting for the actual 6 vector 
the estimated small angle perturbation from the null torque attitude 

« = 6 “ B ~ « 0 (11 ' 13) 

where 

9g = Actual euler angles in (11-14) 

A bm matrix 
*B 


= Estimated euler angles in (11-15) 

A qn matrix 


and the Ag^ and Aq^ matrices are as defined in Section 8.1.1. 
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Similarly, 


6 = 5g - 6 q (11-16) 

where 


= Actual euler angles in Aq^ matrix (1.1-17) 


The implicit assumption is made in Equations (11-13) and (11-17) that at 
most one element of 6 Q or is a large angle; this is true for the nose 
forward and nose down orientations under consideration, with <f> 0 = = 0. 

With the assumed CMG control law the orbiter dynamics are described 
by . 

w = - r 1 ■ [Gx + K<5 + <5T Ap ] (11-18) 

or in terms of the euler angles 

6 = at 5 - I [Gx + K6 + ST Ap ] (11-19) 



11.2 CONTROL DESIGN 

As the pitch axis has the largest potential control problems (highly 
unstable in the nose forward orientation, large variation in null torque 
attitude in the nose down orientation, this- axis will be treated here. 
Figure 11-1 contains a block diagram of a" control system of the type dis- 
cussed in the previous section, with the notation slightly changed to 
emphasize the physical structure. The nonlinear coupling compensation 
terms in Equation (11-12) do not appear in the single axis analysis and 
are in any case relatively small. The control system can bo nndpr^tnnd 
in terms of a conventional attitude loop and two ada 
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Figure 11-1. CMG Pitch Control System 















e attitude loop (K p and K R )- tends to drive e g . to e c = ,e + 
and 0 g to zero. 

© null trim loop (K c ) - provides adaptive attitude command 
correction based on momentum stored. in CMG to drive e g to 
the true null torque attitude despite errors in the es- 
timate of its location e . 

o 

9 momentum unloading loop (Kj)' - provides offset, from adaptively 
determined null torque attitude, based on integral of wheel 
momentum, to utilize disturbance torque for unloading the 
CMG. (Most effective when |Ky Q | is large.) 

When both adaptive loops are used, and the gains are selected to 
yield a stable system, the final steady state values will be e g = e o> 
h = 0 in a constant atmospheric density environment. If the Kj loop is 
omitted, the stored CMG momentum will reach a steady state value different 

from zero, unless 0 = e . In this case the RCS thrusters could be fired 

..oo- 

briefly for momentum unloading. The adaptive loop feedback limiter improves 
the large signal transient response- during initial captive from the RCS 
limit cycle. 

The gains required to yield a particular set of closed loop roots 
can be obtained by matching coefficients in the characteristic equation. 

In the case where .two, roots are complex conjugates, and" two are real, 
the characteristic equation is 


d(s) -= 


(s 2 + 2u ? s + w 2 )(s 


+ X)(s + n) 


( 11 - 20 ) 


and the required control gains are 
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( 11 - 21 ) 


( 11 - 22 ) 


(11-23) 


(11-24) 


Different sets of gains will in general be required depending on 
whether the orbiter is in the nose forward or nose down orientation, and 
whether the operating mode is initial capture from the orbiter RCS limit 
cycle or normal long term stabilization. The considerations entering the 
design for each of these conditions and sample designs are presented in 
the following subsection. 

11.3 INITIAL CAPTURE FROM RCS LIMIT CYCLE 

In general, the largest torque and momentum requirements will be 
during the initial capture transient from the RCS limit cycle rather than 
during steady state operation. Because of the complexity of the system, 
accurate sizing can only be done by simulation. Analysis is useful, how- 
ever, for obtaining a rough sizing estimate and establishing the relative 
importance of the various factors influencing the torque and momentum re- 
quirements. The analysis here assumes a constant atmospheric density 
over the capture interval; density variations are considered in the normal 
operation section. 
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11.3.1 Momentum Estimate 

If no unloading takes place during the transient, the stored momentum 
consists of two terms 

a the excess momentum initially stored in the orbiter due to 
its RCS limit cycle rate. 

s the net momentum absorbed from the disturbance torque due to 
being off the null torque attitude. 

The first term is, assuming the nominal RCS limit, cycle rates of 
‘0.01 deg/sec. 



(11-25) 


— 

■ 

1. 

(1.24 x 10 6 kg-m 2 )(1.75 x 10~ 4 rad/sec) 


216 N-m-sec 

(9.39 x IQ 6 kg-m 2 ) (1 .75 x 10~ 4 rad/sec) 

- 

1639 N-m-sec 

(9.72 x 10 6 kg-m 2 ) (1.75 x 10~ 4 rad/sec) 


1696 N-m-sec 


The peak value of the second term depends on the details of the 
transient response and is therefore more difficult to determine. As a 
rough estimate, the steady state value for the control system describee 
in the previous section will be examined. If Kj = 0, the steady state 


11-9 



momentum stored is, for pitch 




or in terms of the closed loop roots 




2 Z X 
to 




r) 

yy 




For example if 

to = 0.002 rad/sec 

5 = 0.5 

A = 0.001 rad/sec 

K,. = - 33.8 N-m/rad 
ye 

I = 9.39 x 10 6 kg-m 2 

J s/ 

A 

6 - e = 0.2 deg 


(11-26) 


(11-27) 
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the steady state momentum would be 


h 


ss 

Dy 


283 N-m-sec 


The peak -stored momentum due to the disturbance torque is actually 
several times as larger as will be shown by simulation. However } the 
momentum due to the limit cycle rate is. still the largest term. 


11.3.2 Torque Estimate 

The largest control torque is required in pitch in the unstable. 

nose forward orientation. Equation 11-12 shows that the pitch control 

torque (T r = -h ) consists of a feedback term 

by y 


( T Cy^FB - Gx - Kp (e Q 0 B + K c h y + K I f^y dt ^ K R ®B ( "J 1-28) 


and a decoupling term. 


^Cy^DEC ~ ^‘ x ^ 


co + h)] =u) to, (I - I ) + (h 

/ Jy V X)( \ J 


<u 2 - h^ to- x ) (11-29 


As with the momentum, the actual peak torque depends on details of the 
transient behavior -and is best determined by simulation. If the damping- 
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ratio r, is relatively high, however, the peak torque will occur near 
t = 0 and can be approximated by 


^ T Cy)pEAlTM 0 o " " K R 0 B^ + “x^ “z^^xx “ W • (11-30) 

For example, suppose the initial conditions are 
0 Q - 0 B = 0.2 deg 

9g(0) = fc> x (0) = to z (0 ) = 0.01 deg/sec 
and the closed loop roots are selected as 

00 = 0.002 rad/sec 
c = 0.5 

X = n - 0.001 rad/sec 

wi th 

I xx = 1 -24x1 0 S kg-m 2 

1 = 9. 39x1 0 6 kg-m 2 

yy 

I z2 = 9.72xl0 6 kg-m 2 

K = - 33.8 N-m/rad 
yy 
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The position gain is then 


’ K p = 129 N-m/ rad 
and : the rate gain is' 

Kp = 63600 N-m/ (rad/sec) 

The feedback components of the control "torque are in turn - 


‘ M 6 o " 0 b ( 0 ) ). = 0.45 N-m 
• K r e B (0) = 11.1 N 

while the decoupling term, is ■ 

m x ( 0) co z (0) (I xx - I zz )= 0.26 N-m 

The total torque, 11.8 N-m, is clearly dominated by the rate feedback 
term, while the decoupling term gives the smallest contribution. 


11.3.3 Simulation 


To better understand the control system operation and verify 
the sizing estimates of the previous sections, a single-axis simulation 
of the pitch control system shown in Figure 1-1-1 was conducted. The 
parameter values used are listed in Table 11-1, the results are 
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Table 11-1. • Initial Capture Transient Simulation, Parameter Summary 


Symbol 

" — - ■ - 

• Parameter 

Value 

Units 

03 

Bandwidth 1 

> Of dominant roots 

0.002 

rad/sec 

X, 

Damping ; 

0.5 


X 

Null Trim Loop Root 

0.001 

rad/sec 

n 

Momentum Unloading Loop Root 

O.'OOl 

rad/sec 

K P 

Position Gain 

129 

N-m/rad 

k r 

Rate Gain 

6.36 x 10 4 

M-m/( rad/sec) 

K c 

Momentum Feedback Gain 

-2.16 x 10” 5 

rad/(N-m-sec) 

; h 

Momentum Integral Feedback 
Gain 

-8.63 x 10" 9 

rad/ (N-m) 

r yy 

Pitch Inertia 

9.39 x 10 6 

kg-m 2 

K ye 

Pitch Disturbance Compliance 

-33.8 

N-m/rad 

h UM 

Adaptive Loop Feedback Limit 

1500 

N-m-sec 
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summarized in Table 11-2, and time plots- are presented in Figures 11-2 
to 11-4. The two cases run represent the open loop unstable nose for- 
ward orientation, and differ only in the value of the initial rate. For 
added realism, the initial attitude 6^(0) = 0.1 degree differs from both 
the actual null torque attitude 6 q = -0.2 degree and the estimated null 
torque attitude e Q = 0. Both cases cover 6000 seconds, or slightly more 
than one orbit. 

In case 1, the rate initial condition is the nominal RCS limit cycle 
rate of 0.01 deg/sec. As shown in Figure 11-2, the torque reaches a peak 
value of 11.3 N-m at t = 0. This is in reasonable agreement with the 
estimate of the feedback terms in the previous section. The stored mo- 
mentum reaches a peak value of 2580. N-m-sec, of which approximately 1-640 
N-m-sec can be attributed to the transfer of the .limit cycle motion mo- 
mentum from the orbiter'body to the CMG's. The stored momentum is almost 
completely unloaded at the end of the run, while the attitude is slowly 
converging to the true torque null. The nonlinear effect of the adaptive 
loop feedback limiter ;is apparent in the plots'; use of this limiter re- 
duces the peak momentum requirement by 410 N,-m-sec. 

In case 2, the initial rate was reduced by a factor of 10. compared 

to case 1, to 0.001 deg/sec. As shown in Figures 11-3 and 11-4, the con- 
vergence to the true torque null is considerably more rapid, the momentum 

is unloaded sooner, and the peak torque is reduced to 1.34 N-m, while 

the peak momentum is reduced to 429 N-m-sec. 

11.4 NORMAL OPERATION 

In normal operation following the -initial capture transient, the 
major influence on the orbiter null torque, attitude stabilization system 
performance at low al titudes is variations in atmospheric density. The 
design example presented in the previous subsection provides adequate per- 
formance during the capture transient, but is excessively sensitive to dis- 
turbance torque variations with frequency components near orbit frequency. 
Designs with lower sens iti.vi ties will now be, described. 
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Table 11-2. Initial Capture Transient Simulation Result Summary 




Initial Conditions 

Peak Values 

'Case 

Figures 

0 B 

(deg) 

' 0 B 

(deg/sec) 

' 6 o 
(deg) 

A 

■ e 0 ' 
(deg) 

h y 

(N-m-sec) 

% 

(N-m) 

h y , 

(N-m-sec) 

1 

11-2 

0.1- 

• 0.01- 

-0.2 

0. 

0 

11.3 

2580 

2 

11-3, 11-4 

0.1 

0.001 

' -0.2 

0 

.. 

0 

1.34 

429 
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Figure 11-3. ■ Case 2: Pitch Convergence From 0.1 x Nominal Rate 
RCS Limit Cycle in Nose Forward Orientation 
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Figure 11-4. Case 2: Pitch Convergence From 0.1 x Nominal Rate 
RCS Limit Cycle in Nose Forward Orientation 


n_iQ 





11.4.1 Selection of Control Law 


Two idealized control strategies for dealing with periodic distur- 
bance torque variations were described in Section 9.2.2. These were: 

(1) Follow the null torque attitude exactly at all time 

(2) Hold the nominal null torque attitude exactly at all times 

The first strategy requires a precise knowledge of the location of thq instan- 
taneous orbi ter null torque attitude, a condition. that is all but impossible to 
meet in practice. The second strategy only requires the mean or nominal null 
torque attitude, but any error in its estimate results in a secular momen- 
tum buildup in the CMG's. To achieve a practical control design, a com- 
promise must be made between holding the orbi ter attitude steady and 
minimizing the required CMG momentum storage capability. 

The nature of the design depends strongly on the orbiter orientation. 

In the nose toward orientation, even a +50% variation in atmospheric 
density at h = 180 Km results ideally in a momentum perturbation <5h of 
only +167 N-m-sec and an attitude perturbation <$6g of +}.2 deg if' the 
null torque is followed, or a momentum perturbation of +500 N-m-sec and 
no attitude perturbation if the attitude is held. In either case, the 
momentum change is less thanthat required to capture from a 0.001 deg/sec 
limit cycle. On the other hand, in the nose down orientation very large 
attitude and momentum perturbations are possible at altitudes below 250 Km. 
Figure 11-5 shows the peak-to-peak momentum and attitude perturbations 
resulting from +20% variations in atmospheric density using either of 
the control strategies previously described. For all altitudes above 185 
Km less momentum is required to follow the null torque attitude variations 
at orbit frequency than -to hold a constant attitude. However,- this ap- 
proach results in large attitude excursions at low altitudes. 

It was demonstrated in Section 9.2.2 that the momentum storage require- 
ment , is. mipimized by following variations in the null torque attitude with 
frequencies less that / K, J/I w and resisting attitude perturbations caused 
by null torque attitude variations with frequencies greater than / I K y e I / ^ - 
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ANGLE (DEG) OR MOMENTUM (N-m-sec) 



Figure 11-5. Attitude and Momentum Excursions Due to + 20 % Atmospheric 
Density Variation at Orbit Frequency in Nose Down Orien- 
tation With Idealized Control 
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This control strategy, which approximates each of the two idealized con- 
trol strategies over different frequency bands, is adopted for normal 
pitch control-. 

The assumed structure of the pitch control law is shown schematically 
in Figure 11-6. The compensation is more general than in Figure 11-1, with 
a lag (Tp) added to the momentum feedback shaping and integral compensation 
(Kp) added to the attitude loop. Since the null torque attitude is not 
known a priori , the momentun feedback shaping is used to attempt to adap- 
tively seek the true null torque attitude and also unload excess momentum. 
Each of the two orbiter orientations is treated separetely, due to the open 
loop stability differences. 

11.4.2 Nose Forward Orientation 

In this orientation the disturbance compliance K . is dominated by 

y “ 

the gravity gradient slope of 33.92 N-m/rad for small offset angles,' and 
pitch is open loop unstable. A large number of cases were analyzed with 
various combinations of non-zero gains. . Five of the better cases are 



ure 11-7. In case 1, K c = Kj = 0 and there is no adaptive action to seek 
the null torque attitude. The attitude therefore remains nearly constant 
and the momentum grows without bound, indicated by the fact that H(0) = 

In case 2, K c f 0, and a step change in T^ results in a finite net change 
in h because H(0) is finite. Making Kj f 0 in case 3 causes H(0) = 0, 
resulting in no net momentum change for a step change in V The response 
now parallels the desired function at low frequencies, but is much larger 
at high frequencies. In case 4, the lag in the momentum feedback shaping 
is made non-zero. This “rolls off" the adaptive loop resulting in lowered 
momentum perturbation from high frequency disturbance variations, compared 
to case 3. The low frequency response is significantly degraded however. 

In case 5, Tp = 0 again, but the integral gain in the attitude loop, Kp, 
is non zero. Only a small improvement in high frequency response results. 
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Figure 11-6. Pitch Control Law Structure 







Table 11-3. Compensation Parameter Summary - Nose Forward Orientation 


Case 

K p 

N-m/rad 

k r 

N-m/{ rad/sec) 

K C 

rad/(N-m-sec) 

K I 

(rad/sec)/ 

(N-m-sec) 

t l 

sec 

k f 

N-m/{rad-sec) 

«<»„> 

db 

H ( 2w q ) 
db 

Desired Response 

- 


- 

- 


- 

49.9 

53.0 1 

1 

9. 73x1 0 2 
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Figure 11-7. Momentum Response to Disturbance Variations in Nose Forward Orientation 






It is obvious that no one case matches the desired function closely 
over the entire frequency range. Since the disturbance variation is 
basically at a = a) Q , the momentum response should be minimized at a) 0 - 
Case 3 is therefore selected as the nominal gain set, for having the 
smallest H (co ) at w and for all to < 10 <u . H(w Q ) is, however, 55.8 db 
= 615 N-m-sec/N-m or about twice the ideal value. 

Figure 11-8 contains a plot of the attitude response ratio <$e B /ST Ap 
for gain set 3. At low frequencies the ratio closely approximates the 
value 1/K that would result from exactly following the'null torque 

y “ 

attitude. The response begins to roll off at the desired break frequency, 
initially with a slope of -1 and then with a slope of -2. The attitude 
response at o) Q is 1.4 db or a factor of 1.18 greater than in the ideal 
case. As the attitude variations are small in the nose forward orienta- 
tion, the linear analysis is valid over the entire EVAL altitude range. 


11.4.3 Nose Down Orientation 


In this orientation the gravity gradient slope is Kg = - 33.92 for 
small offset angles ( h > 250 Km), and pitch is open loop stable. A 
sai table set of gain is 

Kp = 74.7 N-m-rad 

K r = 2400 N-m/( rad/sec) 

K ■ 6.94 x 10~ 5 rad/{N-m-sec) 
c 

Kj = 3.5 x 10~ 8 ( rad/sec}/ (N-m-sec) 


These gains are much lower than in the nose forward orientation because 
of the open loop stability, but H (o> 0 ) - 55.8 db again, as shown in 
Figure 11-9. At high frequencies, the response is actually better than 
that which would result from holding a constant attitude. 
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.The attitude response ratio is plotted in Figure 11-10.- 

At low frequencies, the ratio again approximates 1 /K^q, but the initi 
slope after the break is now -2, resulting in improved rejection of 
high frequency disturbance variations. The attitude response at ui Q i 
3.2 db or a factor of 1.45 greater than in the idealized case. 


At low altitudes, large momentum and attitude variations may 
occur, and due to the nonlinear nature of the disturbance torques the 
linear analysis is not strictly valid. To verify .performance at low alti- 
tudes, several simulation cases were run with a large angle single axis 
simulation. The results of these simulations are summarized in Figure 
11-11. Above 250 Km, the simulation results, agree closely with the 
linear analysis. Below 250 Km, deviations from the linear analysis 
appear, as expected. The combination of large attitude variations 
and momentum storage requirements may make operation in the nose down 
orientation unfeasible at the low end of the altitude range. 


11.5 MOMENTUM EXCHANGE DEVICE SELECTION 

A. tradeoff analysis of momentum exchange devices capable of 
meeting the control requirements developed in the previous subsection 
is contained in Appendix F. Or the basis of size, weight, power, 
cost and feasibility, double gimbal control moment gyros .rate as the 
most suitable- control actuators for the orbiter null torque attitude 
stabilization system. 

Table 11-4 summarizes the characteristics of the required CMG's 
under two sets of assumptions on allowable operational restrictions. 
The first CMG set allows capture from the nominal RCS limit cycle 
rates and can stabilize the orbiter in the nose forward orientation 
at all altitudes above 180 Km and in the nose down orientation above 
200 Km. The second. CMG set can only be used if the RCS limit cycle 
rates are reduced’ by a factor of 10 and if operation in the nose down 
orientation is not required below 230 Km. 
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Table 11-4. Summary of Double Gimbal CMG Characteristics 


CMG 

Set 

CMG Characteristics 

Operational Restrictions 

Per Unit 

Total for Cluster of 3 Units 

Maximum RCS 
Limit Cycle Rate 
(Deg/Sec) 

Minimum Altitude (Km) 

for Orientation ** 



■ Torque 
(N-m) 

Momentum 

(N-m-Sec) 


Power* 

(Watts) 

Volume 

(m 3 ) 

Nose Forward 

Nose Down 

1 

13.6 

(10 ft-lb) 

2710 

(2000 ft-lb-sec) 

281 

93 

2.70 

0.01 (Nominal) 

180 

200 

2 

1.36 
(1 ft-lb) 

678 

(500 ft-lb-sec) 

136 

45 

0.97 

0.001 (0.1 Nominal) 

180 

230 


' * Peak power at maximum torque output. 

** +20% Atmospheric Density Variation at Orbit Frequency Assumed. 






12.0 ORBITER NULL TORQUE ATTITUDE STABILIZATION SYSTEM DEFINITION 

An overall design of the orbiter null torque atttitude stabiliza- 
tion system is presented in this 'section. The system concept stresses 
autonomous operation, with nearly all control function performed by 
dedicated programmable digital electronics (PDE). 

The orbiter null torque attitude stabilization system maintains 
the orbiter attitude near the null torque attitude, that is, the atti- 
tude at which the total external torque on the orbiter vanishes. The 
system attempts to hold the orbiter roll and yaw attitudes at the esti- 
mated null torque attitude, whi.le adaptively seeking the true null torque 
attitude in pitch. The orbiter RCS is disabled during operation of the 
null torque attitude stabilization system, except for occasional brief 
firings of the vernier RCS thrusters to unload CMG momentum when the 
CMG cluster is near saturation. 

A functional block diagram, showing the major blocks and data flow 
of the orbiter null torque attitutde stabilization system, is presented 
in Figure 12-1. Descriptions and more detailed data flow diagrams of 
the major blocks are contained in the following subsections. 

12.1 COMMAND AND ERROR PROCESSING 

Figure 12-2 shows a block diagram of the command and error process- 
ing performed by the orbiter stabilization system PDE. The only inputs 
required from the spacelab computer are orbiter attitude and ephemeris 
data, an indication of the orbiter orientation number, and a vector of 
estimated null torque offset angles. The orbiter attitude and ephemeris 
data originates in the orbiter GN&C computer, but is relayed through the 
spacelab to maintain a uniform interface. The orbiter orientation number 
indicates whether the "nose forward" or "nose down" nominal orientation 
is desired. The A^ matrix is fixed for each orientation and contains 
only three nonzero (+1) elements. The estimated null torque offset 
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Figure 12-1. Orbiter Null Torque Attitude Stabilization System Functional 

Block Diagram 
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NOTATION: 

A.. = a 3 x 3 DIRECTION COSINE MATRIX 
1J INDICATING THE ATTITUDE OF 

COORDINATE FRAME "i" RELATIVE 
TO FRAME "j" 


SUBSCRIPTS : 

I = INERTIAL (EC1) 

B = ORBITER BODY (GN&C BASE) 

R = ORBIT REFERENCE 
N = NOMINAL ORBITER ORIENTATION 
• T = ESTIMATED NULL TORQUE ORIENTATION 


Figure 12-2. Orbiter Stabilization System Command and Error Processing 


ORIGINAL PAGE lb 
OF POOR QUALITY 









angles are precomputed for each of the two nominal orbiter orientations, 
based on the estimated atmospheric density and orbiter mass properties 
for the particular EVAL mission being conducted. 

The command and error processing outputs are the orbiter rate and 
attitude errors. The rate error output serves to cause the orbiter to 
rotate about the orbit normal, with respect to the inertial (I) frame, 
at orbit rate. This aids in keeping the orbiter locked to the earth 
pointing orbit reference (R) frame. The attitude error output consists 
of the small angle errors between the actual orbiter attitude and the 
estimated null torque attitude. 

12.2 CONTROL LAW 

A block diagram of the orbiter stabilization system control law 
data flow is shown in Figure 12-3. The control law, derived in Section 11.1 
contains two branches. The primary attitude stabilization task falls 
upon the feedback control law. The momentum coupling compensation 
path assists in attitude control by attempting to cancel undesirable 
inter-axis coupling torques that result from the interaction of the 
orbiter inertial rate and system momentum vector. The output of the 
control law is a CMG torque command vector. 

Figure 12-4 shows the' data flow structure of the feedback control 
law. The upper two paths in the Figure implement a conventional 
proportional plus rate attitude control law and operate on all three 
(roll, pitch, yaw) control channels. The lower two paths use propor- 
tional plus integral feedback of the CMG cluster momentum to implement 
loops that adaptively seek the null torque attitude and use the gravity 
gradient torque to unload excess CMG momentum. These adaptive loops 
are used only on the pitch channel. The numerical values of the gain 
matrices (G £ , G^, G^, G^) are different in the "nose forward" and 
"nose down" orientations. 
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Figure 12-3. Orbiter Stabilization System Control law 
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Figure 12-4. Orbiter Stabilization System Feedback Control Law 









12.3 CMG CONTROL AND MOMENTUM MANAGEMENT 

The data flow for 'CMG control and momentum management is shown in 
Figure 12-5. The CMG cluster consists of three double gimballed CMG's. 

As illustrated in Figure 12-6, the CMG's are mounted to have, with zero 
gimbal deflections, mutually orthogonal momentum vectors. Control torques 
on the orbiter are generated by commanding gimbal rates, while momentum 
is stored as gimbal, deflections. 

With all three CMG's operating, there are six degrees of freedom 
in the CMG cluster, and the CMG steering law perforins three functions. 

The primary function, using three degrees of freedom, is to command 
combinations of gimbal rates such that the actual torque on the orbiter 
equals the commanded control torque. The remaining three degrees of 
freedom are used to distribute the individual CMG momentum vectors to 
avoid undesirable momentum configurations and to stay away from gimbal 
mechanical stops. These secondary functions are programmed in such a' 
manner that no net torque on the orbiter results. 

With one CMG failed, there are four degrees of freedom in the CMG 
cluster. Three of these degrees of freedom are used to generate control 
torques, while the remaining degree of freedom is used to. reduce the 
gimbal angle of the gimbal axis that is nearest to a gimbal stop. Mo- 
mentum storage and control torque capability- is reduced with one CMG 

failed, but continued operation, with some performance degradation, is 
possible. 

The CMG cluster momentum is computed from the CMG gimbal angles, 
as measured by the CMG gimbal angle resolvers; Logic is provided to 
indicate when the CMG cluster momentum is approaching a saturated 
condition. To maintain control the cluster must be desaturated with 
torques provided by the orbiter RCS. To minimize the impact of the 
momentum unloading on experiment pointing, minimum on-time 40 msec) 
vernier RCS thruster firings will be used, and transmission of thruster 
commands to the RCS can be inhibited, on a priority basis, by the 
space! ab computer. 
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Figure 12-5. Orbiter Stabilization System CMG Control 
and Momentum Management 
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Figure 12-6. Orthogonal Double Gimballed CMG Cluster 
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APPENDIX A 

SPACE SHUTTLE ORBITER DATA “ 


A. 0 INTRODUCTION 

This appendix briefly summarizes basic space shuttle orbiter 
data that will be required for the Earth Viewing Applications Laboratory 
(EVAL) study. The areas covered include coordinate systems, payload bay 
characterization, internal disturbance sources, and mass properties. The 
MKS system has been adopted for the EVAL study and therefore units are 
converted to MKS where necessary-. 

A.l COORDINATE -SYSTEMS 

Most shuttle orbiter source data is specified in the "orbiter co- 
ordinate system", (x q , y , z ) , illustrated and described in Figure A-l. 
The "station numbers" appearing on dimensional drawings are the (x o> y Q , 
z ) coordinates, in inches. This frame is convenient for specifying lo- 
cations on the orbiter because the (x Q , y Q , z Q ) coordinates are always 
fixed in the orbiter. 

A more useful frame for dynamics is the orbiter body coordinate 
system (x g , y g , z g ) illustrated in Figure A-2. This frame has its origin 
at the orbiter CG and body-fixed axes oriented in the standard "airplane" 
arrangement. The orbiter body coordinates are related to the orbiter 
coordinates by 
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where (x Q , y Q , z q ) is the CG location in orbiter coordinates. 
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TYPE: ROTATING, ORBITER REFERENCED 

ORIGIN: APPROXIMATELY 200 INCHES AHEAD OF THE NOSE AND APPROXIMATELY 

400 INCHES BELOW THE CENTERLINE OF THE PAYLOAD BAY 

ORIENTATION AND LABELING: 

THE X AXIS IS PARALLEL TO THE CENTERLINE OF THE PAYLOAD BAY, 
NEGATIVE IN THE DIRECTION OF LAUNCH 

THE Z AXIS IS POSITIVE UPWARD IN LANDING ATTITUDE 

THE Y COMPLETES THE RIGHT-HANDED SYSTEM 

THE STANDARD SUBSCRIPT IS 0 

Figure A-l. ORBITER COORDINATE SYSTEM 
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NAME: ORBITER BODY COORDINATE SYSTEM 


ORIGIN: ’ - ’ORBITER CG 

ORIENTATION: ' THE Xr’ AXIS IS PARALLEL TO THE 

CENTERLINE OF THE PAYLOAD' 

BAY ; POSITIVE IN THE DIRECTION 
OF THE NOSE OF THE ORBITER' ' 

THE Z B AXIS LIES IN A PLANE PAR- 
ALLEL TO THE PLANE OF SYMME- 
TRY AND IS PERPENDICULAR TO 
THE AXIS AND IS POSITIVE DOWN 
IN THELANDING ATTITUDE 

THE Y r AXIS COMPLETES THE 
RIGHT-HANDED SYSTEM 

CHARACTERISTICS: ROTATING RIGHT-HANDED SYSTEM ■. 


Figure A-2. ORBITER BODY COORDINATE SYSTEM 
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A. 2 PAYLOAD BAY 


The instrument pointing system (IPS) will be located in the payload 
bay, illustrated in Figure A-3. The payload envelope is restricted to 
a cylinder centered at 


x - 942 in ^ 23.9268 meters 
o 

2 . = 400 in a 10.16 meters 
o 

with overall length 720 in = 18.288 meters and radius 90 in = 2.288 
meters. A 180 degree lateral field of view is available at z Q = 427 in * 
10.8458 meters. The longitudinal fields of view at the payload bay center 
point are illustrated at the top of Figure A-3. 


A. 3 INTERNAL DISTURBANCE SOURCES 

The two largest disturbance sources originating within the orbiter 
during attitude hold operation are crew motion and vernier reaction control 
system (RCS) thruster firings. External disturbances are mission dependent 
and are discussed in Appendix D. 


A. 3.1 Crew Motion 

The worst case crew motion force profile is illustrated in Figure A-4. 
In general, this force profile can be applied at any point within the 
orbiter that the crew has access to. 
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Figure a-3. Payload bay clearance and view angles 
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A. 3.2 Vernier RCS Thrusters 


Fine orbiter attitude control is achieved with six vernier thrusters 
with thrust 25 LB f = 111.2 N, specific impulse I = 228 sec, and a minimum 
on-time of 0.04 sec. The thruster locations and force components are listed 
in Table A-l. Thruster torques depend on. the overall orbi ter/payload CG loca- 
tion and are therefore somewhat configuration dependent. The orbiter motion 
about each axis is a limit' cycle with a selectable deadband of either +0.1, 
+ 0.5 or +1.0 degrees. 


A. 4 MASS PROPERTIES 

Orbiter mass properties are highly configuration dependent. Table A-2 
lists typical weights, less the payload. Definitive inertia values are 
not readily available. A typical set of principal axes inertias, including 
the payload, are 

I = 1.24 x 10 6 kg-m 2 

AA 

I = 9.39 x 10 6 kg-m 2 

I = 9.72 x 10 6 kg-m 2 

The corresponding mass is 

m = 97000 kg 

A. 5 REFERENCES 

A-l. NASA TM X-64978, "Experiment Pointing Subsystems (EPS) 

Requirements for Spacelab Missions," Marshall Space Flight 
Center, December 1975. 

A-2. JSC 07700, Volume XIV, Revision D, Change 13, "Space Shuttle 

System Payload Accommodations" , Johnson Space Center, 26 November 
1975. 



Table A-1. Vernier RCS Thruster Locations and Force Components 


Thruster # 

Coordinates * (M) 

FORCE COMPONENTS 

(N) 

x o 

-V 

z o 

F x 

F y 

F z 

1 

8.260 

- 1.228 

9.166 

0 

72.95 

- 84.96 

2 

8.260 

1.228 

9.166 

0 

- 72.95 

- 84.96 

3 

39.751 

- 3.641 

11.659 

0 

111.2 ' 

0 

4 

39.751 

3.641 

11.659 

0 

-111.2 

0 

5 

39.751 ' 

- 2.870 

11.568 

0 

0 

111.2 

6 

39.751 

2.870 

11.568 

0 

0 

- 111.2 


* In Orbiter Coordinate System, CG is approximately 
(x 0 , y Q , z Q ) CG = (27.4, 0, 9.5) meters 








Table A-2. 


Typical Orbiter Height Suawary 



WJUjt-lbSi 


Inert 

150,000 

68, 039 

Non-Propulsive Pluid 

4,700 

2,132 

Personnel Group 

2,411 

1, 094 

RCS Peopellant (Fully loaded) 

7,391 

3,353 

OMS Propellant (Fully loaded) 

_25 x _1 04 

11,377 

Orbiter Ht./Less .Cargo 

189,608 

86,004 

ORS Payload Bay Kits (1) 

14,255 

6,466 

(2) 

27,631 

12,533 

(3) 

41 ,009 

.18,601 

Cryogenic Wt. Per Kit 

1 ,800 

861 
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APPENDIX B 
IPS DATA 


B.O INTRODUCTION 


This appendix briefly summarizes the currently available Instrument 
Pointing System (IPS) data that is relevant to the Earth Viewing Applications 
Laboratory (EVAL) study, IPS, under development by ESA for the Spacelab 
program, is a precision, softmounted, three-axis gimballed platform primarily 
intended for pointing experiments at inertially fixed (stellar or solar) 
targets. The EVAL study assesses the suitability of the IPS for earth view- 
ing applications. 

Data on the following IPS characteristics is included: 


e mechanical structure- 
8 coordinate systems 
© mass properties 
© payload characteristics 
o drive characteristics 
e softmount characteristics 
e gyro characteristics 
© star tracker characteristics 


B.l MECHANICAL STRUCTURE 

The overall IPS mechanical configuration, including the pallet and 
payload, is shown in Figure B-l . Much of the equipment illustrated (payload/ 
gimbal separation mechanism, payload clamp, jettison device, softmount 
clamp) is used only during shuttle ascent/descent or emergency* operation, and 
is not directly related to payload pointing. The star tracker assembly is 
not shown because its location and orientation for earth viewing applications 
has not been determined. 
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Figure B-l . IPS Mechanical Configuration 



Figure B-.2 is a more detailed illustration of the gimbal and soft- 
mount configuration. The ginibal is shown in the stowed position., with 
the payload integration ring normal to the pallet floor. For pointing 
operation, the nominal gimbal orientation is rotated 90 degrees about the 
elevation axis (payload integration ring parallel to pallet floor). During 
shuttle ascent and descent, the softmount is .locked closed with the soft- 
mount clamp and the payload i's locked to the pallet with a pay-load clamp. 

In order to center various size payloads properly in the payload clamp, 
the softmount/gimbal assembly position must be adjusted for each payload 
by sliding along a rail on the gimbal bottom structure and using a replace- 
able column for height adjustment. 


B. 2 COORDINATE SYSTEMS 

Figure B-3 contains a simplified, exploded view of the IPS gimbal/ 
softmount in the null pointing orientation. .A coordinate system is 
defined in each of the seven bodies shown. The payload attitude relative 
to the orbiter is described by the three gimbal angles 


a = azimuth angle (rotation about z^) 

B = cross-elevation angle (rotation about x^) 
Y - elevation angle (rotation about y £ ) 

and the three (small) softmount angles 

e x = x sofcmount angle 
6 y ~ y softmount angle 
e z = z softmount angle. 
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Figure B-2. IPS 
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Figure B-3. Definition of Coordinate Systems 
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The payload unit vectors are related to' the orbiter body unit vectors by 
the transformation 


r p = A PB r B 

where .the transformation Apg pan be factored into the product of transformations - 

A PB = A PE a ex a xa a au a ul a lb 

Assuming zero IPS/orbiter mounting misalignments, the orbiter 
body (Appendix A) to softmount base transformation is 



The softmount base to softmount upper assembly rotational transformation 
is the small angle transformation • 



The softmount also allows three-axis translation between the L and U frames. 
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The gimbal angle transformations are, for azimuth 


W 


cos a sin a 0 


-sin a COS a 


for cross-elevation 


A XA 


cos 3 sin 3 


-sin .3 cos 3 


and for elevation 


EX 


cos y 
0 

sin y 


1 

0 


-sin y I 

0 

COS Y 


Assuming zero payload mounting misalignment, the final transformation is 




0 

1 

0 . 


1 

0 

0 
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This transformation is for agreement with the IPS documentation (References B-l 

A 

B-2) convention of defining Xp as the payload LOS. 

The gimbal angles have the ranges (Reference B-2) 


azimuth: 

cross elevation: 
elevation: 


180 < « < 180 degrees 

- 70 <_ g <_ 70 degrees 

- 70 <_ y _< 90 degrees 


The gimbal axes are specified to intersect within a 1.0 cm radius sphere, 
and are required to be aligned orthogonally to within 20 arc min with 
an uncertainty of less than 1 arc min. The softmount displacement limits 
are presented in Section B.7. 


B. 3 MASS 'PROPERTIES 

A preliminary mass summary, from Reference B-l, is contained in 
Table B-l with further details in Tables B-2 and B-3. Reliable inertia data is 
not available. Preliminary inertia estimates for some of the structural 
sections illustrated in Figure B-3 are contained in Reference B-3. 


B.4 PAYLOAD CHARACTERISTICS 

Figure B-4 illustrates the payload size/mass envelope derived 
from Reference B-2. The limits given for payload radius, length, and CG 
location are largely based on the requirement to center the payload CG 
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Table B-l. Mass Summary 


Assembly 

Subassembly/Unit 

No. 

Mass 

(Kg) 

Comments 

Gimbal 



256 

See Table B~2 for 

Structure 

Payload/Gimbal 

Separation 

1 


breakdown 


Oettison Device 

n 


- 


Softmound Clamp 

H 



Drive 




See Table B-3 for 
breakdown 


Torque Motor 

B 




load By-Pass Drive 

B 




Resol ver 

1 



Thermal 



10 

Distributed throughout 

Control 




qimbal system 

Payload 



144 

Mounted on pallet 

Cl amp 

. ‘ 

■B 

• 

Attitude 

Measurement 

• 

i 

r 

57 

— 


Star Sensor 

3 

27 

\ May not apply for 


Solar Sensor 

(1) 

2 

[ earth viewing 


Optical Sensor Housings 

1 

18 

J missions 


Gyro Package 

1 

12 

Mounted on payload 
integration rinq 

Power 

Electroni cs 


1 

27 



PE Unit 

2 

21 

Mounted on pallet 


Control Panel 

2 

6 

External to pallet 


B-9 
















Table B-l. Mass Summary (Cont'd) 


Assembly 

Subassembly/Unit 

No, 


Comments 

Data 

Electronics 


i 

24 



Data Distribution Unit 

.1 

5 

Mounted on inner 
gimbal 


Harness 

i 

19 

Distributed through- 
out gimbal system 

MARGIN 



106 


TOTAL 

• 

750 
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Table B-2, Gimbal Structure Estimated Mass Breakdown' 


Integration Ring ‘ 

20 

kg. 

Outer Gimbal 

35 

Kg 

Payload/Gimbal Structure - Separation 
Mechanism 

12 

kg 

Inner Gimbal 

10 

kg 

Yoke 

15 

kg 

Softmount Clamp 

50 

kg 

Softmount Dampers (6) 

7 

kg 

Jettison Device with Harness- 
Separator 

35 

kg 

Plate 

25 

kg 

Gimbal Bottom Structure Including 
Rails 

40 

kg 

Miscellaneous 

7 



256 kg- 


Table B-3. Drive Assembly Estimated Mass Breakdown 


Item 

Mass/Assembly 

2 Torquers 

6 kg . . 

2 Resolvers 

4 kg 

2 Angular contact bearings 
1 Radial contact bearing 

1.5 kg 

1 Motor for load by-pass mechanism 

4.5 kg 

1 Case with part of by-pass 

12 kg' 

1 Shaft with part of- by-pass 

9 kg 

Miscellaneous 

5 kg 

Total per assembly 

42 kg 

Total for 3 Drive Assemblies 

126 kg 
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Integration 

Ring 


Quantity 

Definition 

Range 

Units 

nip 

Payload Mass 

. 0 •* 3000 

Kg 

■ 

Payload length 

1 6 

m. 

Y 

Payload radius 

0.25 ->1.5 

m 

*x 

Payload eg location along payload 
integration ring center line 

0.5 -> 3 

m 

h 

Payload eg location normal to 
payload integration ring center 

line = Uy + £^) 1/2 

(0 -> 0.3)* 

m 


* Reference B-2 is inconsistent on this point. Value given is from 
Table 3.2-1 of Reference, while Section 3. 3. 4-1 implies maximum 
is 0,05 m. 

Figure B-4. Payload Envelope 
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in the payload clamp assembly during orbiter ascent/descent. This re- 
quirement imposes minimum as well as maximum payload dimensions. 

Table B-4 contains physical data on the two reference payloads 
defined in Reference B-2, along with the computed payload inertias about 
the gimbal center of rotation (CR) expressed in r^ coordinates. The 
large CG/CR offsets greatly increase the inertias about the CR. 


B.5 DRIVE CHARACTERISTICS 

Each of the three gimbal drives is identical with the following 
characteristics (Reference B-2). 

® 2 Torque motors for fine pointing and slewing 


- type 

- stall torque 

- weight 

- ripple 

e 2 Resolvers 

- type 

- accuracy 

e 3 Ball Bearings 

- maximum friction torque 

0 Cable follow-up 

- flex lead spring torque 

- hysteresis 


brushless DC 
10 N-m 
3 kg 

< 7 % of mean torque 

single speed/multispeed 
3 arc min 

0.03 N-m per axis 
/ 

0.5 N-m for 90 deg rotation 
~0.005 N-m for 1 deg limit cycle 


The two torque motors can be operated in parallel to yield -up to 20 Nm 
torque. The drives can be locked with a self-contained load by-pass 
mechanism. 
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Table B-4. Definition of Reference Payload Physical Data 



Large Payload 

Small Payload 

Mass 

2000 kg 

200 kg 

Dimensions : 



Length 

4 m 

1.5 m 

Diameter 

2 m 

1.0 m 

Moment of inertia about 



payload CG: 

about LOS axis 

1000 kgm 2 

25 kgm 2 

about axes perp. to LOS 

1200 kgm 2 

20 kgm 2 

CG offset from CR of 



gimbal axes: 
along LOS 

2.50 m 

1 .50 m 

perp. to LOS 

0.30 m 

0.1 0 m 

Moment of inertia about 
CR in r F - coordinates 



i» ' 

13880 

472 

I „ 

0 

0 

xy 



*xz 

0 

0 

V 

13700 

470 

*yz 

-1500 

-30 

I,_ 

1180 

27 

zz 


* Assumes payload CR/CG offset perpendicular to LOS is along y^. 
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Figure B-5 contains a block diagram of the torquer drive system 
(Reference B-l). The primary input is the torque command from the DDU 
(Data Distribution Unit), and the motor torque nominally equals the 
commmand torque. The single speed resolver is used to derive gimbal 
angle and rate (with an accuracy of +_ 2 deg and + 0.15 deg/sec) for 
use in the angle and speed limit circuits. The speed limit is linearly 
decreased as the gimbal angle reaches its motion limits. The multi- 
speed resolver is used to commutate- the brushless DC motor. 


B.6. SOFTMOUNT CHARACTERISTICS 

The softmount consists of six spring/dampers. Its purpose is to 
support and isolate the IPS gimbal system from high frequency shuttle 
disturbances, such as thruster firings and crew motion. Figure B-2 shows 
the location of the softmount in the IPS base, while Figure B-6 shows the 
arrangement of the individual spring/dampers. The softmount attachments 
use spherical bearings. Figure B-7 shows the details of a single spring/ 
damper. 


The rotational limits of the spring damper assembly are (Reference B-2) 


e + 5 mm translation (3-axis) 

0 + -3 deg rotation (3-axis) 

The translational spring constants for the assembly are (Reference B-3) 

K Tx = 2000 N/m 
K Jy = 2000 N/m 
K Tz = 800 N/m 
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Figure B- 


Q & 


Softmount 

Base 


. Gimbal Structure/Softmount Suspension Arrangement 
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The translational damping constants are 


B Tx = 600 N/On/sec) 

= 600 N/(m/sec) 

B.^ = 240 N/ (m/sec) 

The rotational spring constants are 

K Qv - 25 N-m/rad 
K H 25 N-m/rad 

°y 

K a ?250 N-m/rad' 

OZ 

The rotational damping constants are 

B ov s 7.5 N-m/(rad/sec) 

B ey - 7,5 N-m (rad/sec) 

B - 75 N-m (rad/see) 

The z-pxis spring/damping constants depend on the skew -angle shown 
in Figure B-6. The value of used for the above figures is not 
available. 


B.7 GYRO CHARACTERISTICS 

The inertial - attitude and rate of the IPS payload is measured with 
a rate integrating gyro (RIG) package mounted on the payload integration 
ring. The gyro package has the following characteristics: 
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Type of gyros 


Four strap-down pulse rebalanced, single 
DOF RIG 

Hamilton Standard RT^- 1139 
Bandwidth 30 Hz 

Gyro Output Noise 


Frequency Band 
(Hz) 


RMS Noise 
(arc sec/sec) 


0,5 to 3 
0.5 to 10 
0,5 to 50 


0.09 found by graphical 
0,43 integration of PSD 
2.99 


Short Term Random Drift (10 min): 0,005 deg/hr 


Max Input Rate (operational) 
Pulse weight (low mode) 

Pylse weight (high mode) 
Configuration of input axes 

Reliability 


2,5 deg/sec 
0.0088 arc sec/pulse 
0,07 arc sec/pulse 
3 orthogonal 
1 eqyiangul ar skewed 
0,9999 (7 days) 
0.9330 (1 year) 


Rate is calculated from the incremental change of attitude (number of 
pulses) per 0.01 second sampling period. Rate is integrated to obtain 
attitude data. Attitude is updated periodically using star tracker data. 


B. 8 STAR TRACKER CHARACTERISTICS 

Reference B-l (Section 3. 4. 2. 2.1) states that attitude information for 
gyro updates for earth missions will be provided by the payload. For 
completeness, however, data on the star tracker proposed for stellar 
missions is presented below. 
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Proposed FHST supplier: 

Type of Sensor 

FOV 

IFOV 

Lens aperture 
Acquisition Probability 
False Acquisition 
Bias Calibrated Accuracy 
Uncalibrated Accuracy 
Noise equivalent angle 

7 threshold levels from +2.5 Mv to 9.3 

Maximum tracking rate 

Mass per Set' (1 FHST, 1 BOD, 1 FEB, 

1 IC) 

Power per Set 
Size 


Honeywell Inc. 

HEAO-B with additional clock network 
2° x 2° 

2 arc min x 2 arc min 
11.68 cm 

98.7% ifor 8 Mv star with 
2.7%^ + 10 Mv background 

0,75 arc sec 
0.5 arc min 

0.76 arc sec for 8 Mv star with + 10 Mv 
background averaged over 0.64 sec 

3 arc min/sec 

9 kg 

4 W 

60 x 21 x 21 cm 
17.5 x 10 x 10 cm 


FHST 

BOD 


(FHST = fixed head star tracker] 
(BOD = bright object detector) 
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APPENDIX C 

DERIVATION OF EARTH POINTING RELATIONS 


C.O INTRODUCTION 

This appendix contains a derivation of the gimbal angles, rates 
and accelerations required to point an orbiting payload at an earth- 
fixed aim point. The derivation is divided into two independent parts: 
■ orbit relations and gimbal relations. The orbit relations yield 
the aim vector in attitude reference coordinates, and do not involve 
the gimbal angles. The gimbal relations depend only on the orientation 
and order of rotation of the gimbal axes with respect to the attitude 
reference frame. For the TPS, it is assumed that the elevation and 
cross-elevation axes will be used for earth pointing. The derivation 
therefore assumes zero azimuth angle. The softmount angles are also 
assumed to be zero. Even with these simplifying assumptions, the 
resulting equations are sufficiently complex to preclude hand calcu- 
lationof numerical results. A computer program for performing these 
calculations is described in Section C.3. 


C.l ORBIT RELATIONS 

The position of the earth-fixed aim point, with respect to the 
spacecraft, is in vector notation 


P = 


where 


P = vector from spacecraft to aim point 
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' r^ = vector from geocenter to aim point 
Fj = vector from geocenter to spacecraft 

It is required here to find the attitude reference frame components 
of p\ and their first and second time derivatives, as explicit time 
functions. As a preliminary step, five coordinate frames are defined 
in Table C-l and. Figures C-l and C-2. 

The vector F_ in attitude reference coordinates fs simply 

O 

F $ = - r S Z R ' - (r 0 + h) Z R 


where 


r<. = radius of S/C orbit 
r Q = radius of earth 
h = s/c altitude 


The vector r^ is, for a given aim point, fixed in the geographic frame 



To express this vector in attitude reference coordinate components, 
the following series of transformations is required 
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Table C-l . Definition of Coordinate frames 


1 

Frame 

Subscript 

Origin 

Type 

* 

Orientation 




rotating 
(fixed in earth) 


Geographi c 

9 

geocenter 

A ‘ 

x in equatorial plane along prime meridian 
a y 

Zq along earth rotation axis 

Earth Centered 
Inertial 

i 

geocenter 

inertial 

A ' 

x. in equatorial plane in direction of vernal equinox 
2 . along earth rotation axis 

Nodal 

n 

geocenter 

inertial 

in equatorial plane along line of ascending node 
of S/C orbit 

A 

z normal to orbit plane 
n 

Orbital 

0 

geocenter 

rotating 

- x in orbit plane pointing towards S/C 

aO 

z q normal to orbit plane 

Attitude 

Reference 

R 

S/C 

rotating 

A *' 

x R in orbit plane along velocity vector 

z pointed at geocenter 
K 


A 


* y completes right handed set 




Symbol 


i (subscript) 

EC I frame 

g (subscript) 

■ Geographic 

X A 

Geographica 

^A 

Geographica 

x e 

= X_ + to_ t 
o e 

x o 

Angle betwe< 

“e 

Earth rotat 


Figure C-l . 


orth 

Pole 

A 



inition of Aim Point Unit Vector 
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Symbol ■ , 

Definition 

i (subscript) 

ECI frame 

n (subscript) 

'Nodal frame 

o (subscript) 

Orbital frame 

. R (subscript) 

Attitude Reference frame 


Longitude of ascending node 

i 

Orbit inclination 

•V 

= V + 61 t 


0 0 

v o 

orbit angle at t = 0 

w b 

orbit rate 












"^Ro ^"on 




In matrix form, these transformations are 


eg 


COS X, 


sin X, 


- sin x. 0 
e 


cos X. 


0 

1 


ne 


cos 

ft 



sin 

ft 


sin 

ft 

cos 

i 

cos 

ft 

cos i 

sin 

ft 

sin 

i 

- cos 

ft 

s i n ■ i 


on 


cos v 


sin v 


sin v 
cos v 
0 


Ro 


0 

0 

1 


1 

0 

0 


0 

- 1 
0 


0 

sin i 
cos i 
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The vector P can now be expressed in attitude reference coordinate 
components as 


“ T Rg ^g ' ^ r S^ R 


Performing the indicated operations, yields the result 


p = r 
x o 


cos sin v cos (si - X a - A fl ) - cos <p n cos i cos v sin (si - x p - A fl ) 


+sin ^ sin i cos v 


e A A 


e A A' 


(C-l) 


p = r 

y o 


-cos sin i sin-(ft - A. o - A fl ) - sin <J> fl cos i 


(C-2) 


e A A ; 


p z = r o 


-cos <j> fl cos v cos (ft - X Q - \ A ) + cos <j> fl cos 1 sin v sin (si - x p - A fl ) 


e A A' 


-sin <p A sin i sin v 


+ r c 


(C-3) 


where 


. v = v + w t 
o o 


A ~ a +w t 
e o e 


(C-4) 
C C —5 } 
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and 


w o = t GM /^ V2 (C-6) 

6 3 ? 

GM f 0.3986 x 10 Km /sec = gravi.tat.ional constant (C-7) 

r s » r 0 + h (C-8) 

r Q = 6371 Km ~ mean earth radius (C-9) 

-5 

(o = 7.27 x 10 rad/sec = earth rate (C-10) 


Substituting for v and \ & and taking derivatives yields the first 
derivatives of the aim vector components 


p x = r Q [-(« - co Q cos i) cos 4 > a sin v sin (n - - X A ) 

+ (o) e COS i - u) Q ) cos <f> A cos v cos (q - - \ A ) 

- u 0 sin <j> A sin i sin v] (C-ll ) 

Py = r 0 [w e . cos <j> A sin i cos (fi - x g - x A )] (C-12) 


P z = r Q E-(“ e " » 0 cos i) cos 4 > a cos v sin (fi - X g -- * A ) 

-(oj 6 cos i - a) Q ) cos (j > A sin v cos (q - - X A ) 

- a) Q sin <j> A sin i cos v] (C-13) 
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and the second derivatives of the aim vector components 


P„ = 


r 0 j C (» 0 + " 2 w o w e cos ^ cos S1 * n v cos ^ ' X e " x fl ) 


2 2 

+ [(w Q + u e ) cos i ■ 2 ii o to e ) cos <}>£ cos v. si : n (si - x 


V 


- w Q sm. ^ sin i cos v 


(C-l 4 


p v = r 0 ^e cos h sin 1 sin ( fl " K - * fl )3 


• A' 
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\ Z 2 

= r 0 [(o) 0 + <D e ) - 2 w Q co g cos i] cos 4 > a cos v cos (n 
2 2 

“•[(« 0 + C 0 g) cos i - 2 a u ] cos <j> A sin v sin (n 


A e “ x fi) 
X e " V 


+ o> Q sin 4 > a sin i sin v 


(C-16) 


The distance between the S/C and the aim point and its derivatives, are 


= |p| - fo - i + 2 * s f z 3 


1/2 


(C-17! 


p = r s p z /p 


(C-l 8) 


»* •• t 

P = r s [p p 2 - p 2 p]/ p ‘ 


(C-l 9) 
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C.2 GIMBAL RELATIONS 


With the shuttle orbiter flying in an inverted attitude, the 
attitude reference and orbiter body frames are related by (Appendix A) 

0 
-1 
0 

With zero softmount and azimuth gimbal angles, the payload LOS unit 

A • 

vector Xp has the attitude reference coordinate components (Appendix B) 

A A 

n x * Up * x R ) - - -sin Y 

' A A 

n y * (x p • y R ) = - sin 3 cos Y 
n z = (x p * z R ) = cos 3 cos y 

The condition for alignment of the payload LOS along the aim vector is 

A A 

P = p/p = Xp 
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Resolved into attitude reference components, this vector equation is 
equivalent to the three scalar equations 


p x /p = = ~ s "* n T 

p y / p = n y “ - sin 8 cos y 
p z /p = n 2 = cos $ cos y 


The required gimbal angles are therefore 


y = sin"^ [- p v /p] 


3 = tan -1 [- P y / p z l 


The first derivatives of the gimbal angles are 


y =. 


p p x _ p p x 


p COS Y 


• 2 r p y P z p v p Z t 

3 = cos^ 3 [ V 2 ^ J 


(C-20) 


(C-21) 


(C-22) 


(C-23) 
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while the second derivatives are 


^ *• *» ^ , 

. ^ p cos y[p,P x - p P X 1 - [p P x - P p ] [2 P pcos Y + p p x y] 


Y “ 


X r r X- 

4 2 

p COS Y 


(C-24) 


f 6 = - 2 8 sin 8 cosg ( - y Pg ■ “ ^ P - ^ 


2 / p y p z " p y p z 
+ cos 6 - — 2 ”^ — - 

p z 


j Li 


z ^ p y p z ~ p .y p z ^ 
3 
p z 


( C— 25 ) 


C.3 PROGRAM AIM 


Program AIM, written in FORTRAN for timeshare use, performs the 
calculations represented by 'Equations (C-l) to (C-25). The program 
listing, contained in Section C.3. 2, is intended to be largely self- 
explanatory. An option exists to limit printout and computation of 
peak values to the time period in which the aim point is within a 
cone of half angle CONE around the local vertical. It should be 
noted that the assumption X Q = 0 does not really restrict the 
generality of the program, as A Q always appears together with A^ and 
n in the grouping a - A g - a^ = n - t - a q - a^. 


C.3.1 Example 

To illustrate the use of the program and the type of results 
obtained, an example case is presented. The example uses the following 
parameter values 
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ALT 

RLAhOA 

Phi 

= 2. Lit + L2> 

- -i . b t 

s 0 » .» 

ALPriAG 

• C « y 

Oh EG AO 

= ' o.. 

P INC ■ 

=- 6 . 0 E + 0 1 * 

C unt 

= o • Ot + f 1 j 

ISTEP • 

= it 

$ F'NO 



Physically, the S/C is in a 200 Km orbit with 60 deg inclination, 'at 
its ascending node at t = 0. The aim point is on the equator, sepa- 
rated 1.5 deg in longitude from the S/C ascending node at t = 0. 


Figures C-3 and C-4 contain time history plots of the elevation 

and cross-elevation angles, rates, and accelerations (units are deg, 

2 

deg/sec, and deg/sec ) obtained from the data on TAPE13.with ICONE = 
60 deg, points are only plotted for times when 


Xp * = cos y cos 3 _> cos (60 deg) 


that is, for 


- 55 < t < 35 sec 


Action must be taken before the aim point enters the 60 degree half- 
angle cone to ensure the proper combination of angles, rates, and 
accelerations at t = - 55. 

The peak value output on TAPE6, with explanatory annotation, is 
presented below. 
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-MINIMUM 

-5?.00G0G00C 


MAXIMUM 

•35.00000000 t 


-4 9.94&34430 

'53. fc9-5443ol 

Y 

*■ 

*57495533 

1.67338761 


Y 

-**03238377 

.03277787 

• * 
Y 


3 1.7 8995201 

3 t » 3 o :/5t 0 84 

S 


12547415 

— . 0 1 7 i t 667 


£ 

- *00120794 

COGvi 2 52 

3 


250.67326736 

42U.C3602473 

P 



Note that the peak values of g and 8 are more than an order of magnitude 
less than the peak values of y and y 


C-3.2 Listing of Program AIM 

A listing of program ArM follows. 
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00100 PROGRAM AIM t INPUT, TAPE5= INPUT, TAP Eb>TAP El 3) 

OOliG CUMMON/R HJBUFK/BUM 241 ) 

00120 C 

00130 C **** PROGRAM TO COMPUTE GIMBAL ANGLES, RATES, AND A CC fc L F p A T I ONS 
00140 C REQUIRED TO POINT IPS PLATFORM AT E ARTH-F IXtD AIM POINT 

00150 C 

00160 c »*** computes time history ion tapei3) fop plotting 

00170 C AND PEAK MAXIMUM AND MINIMUM VALUES (UN FAPEb) 

00160 C 

00190 C **** REQUIRED INPUT DAI A : (ANGLES IN DEGREES) 

00200 C ALT' 3 ORBIT ALTITUDE (KM) 

00210 C R LAND A * LONGITUDE OF AIM POINT (IN GEOGRAPHIC COORDINATES) 

00220 C PHI « LATI1U0E OF AIN POINT (IN GEOG R APHIC COORDINATES) 

00230 C ALPHAO = URBI1 ANGLE FROM ASCENDING NDufc AT T=0 

00240 C UhEGAO = LONGITUDE OF ASCENDING NUDE (ECI COORDINATES) 

00250 C P INC = ORBIT INCLINATION . ‘ ‘ 

00260 C CUKE * HALF CONE ANGLt AROUND LOCAL VERTICAL FOR WHICH 

00270 C DATA IS TO BE CALCULATED 

0C28Q C 1STEP » COMPUTATION TIKE STEP (SEC) 

00290 C 

00300 C +*** FOLLOWING DATA IS PRESET FOR EARTH ORBIT 

00310 C RO = EARTH RADIUS (KM) , __ 

00320 C • WE ' EARTH ROTATION RAPE (RAD/SEC) 

00330 C CM 55 GRAVITATIONAL CONSTANT ( K + 3 / S EC **2 ) 

00340 C 

0 0 3 60 C **** THESE ASSUMPTIONS HAVE BEEN MADE * 

00360 C - CIRCULAR ORBIT 

00370 C - SPHERICAL LARTH 

00380 C - LAMD AG*0 (!.£♦ PRIME MERIDIAN ALONG VERNAL EOUINOX 

00390 C ' AT T*0) 

00400 C - ZERO AZIMUTH AN D SOFTMUUNT ANGLES 

00410 C ’ 

00420 CALL LINESIZ< 24C ) 

00430 MAMELIST/LlST/AlT,PLAMDA,PHi,ALPHAO,OMLGAO,KINC,COMF,ISTEP 

00440 DA Fa RO, WE , GM / 63 71 . , 7. 2722032E-5, 0.3966E6/ 

0 04 50 DATA GM IN , GDM I N, GDDN I N, BM 1 N , BOM I N, B ODM I N / 6* 1 . E 5 u / 

0 0460 DATA G M A X , G D H A x , G 0 D M A X, BM AX, BL'M A X , H 9 DM AX / 6*-l . E 60/ 

0G47U DATA KM I N , TMI N , k M A X , T HA X / 2 * 1 . l u * -UfSO/ 

00460 DAfA PI/3. 141^926635897^3/ 

00490 DA ( A CGufc , 1 S T E P / oO ♦ , 1 / 

00500 P6M>l) 13 " 

0051G P ADDE G s 160 . /PI 


IS O 

teg 

If 

£ 

fc c 
sHS*- 


ORIGINAL PAGE g» 

OF POOR QUALITY 

* 


' 00520 D E CP A 0* 1 » /P. ADO if G 


00530 

‘ 1 1 

>• U-J.1 / r. jjU CM 

FfcAD(5>LlST) 

00 5 90 


WPITE(6,LIST) 

OC 5 50 


WRIT E (13, 101 ) A L f jRLAMOAjPHIjALPHAO»OMEGAOjP INC 

OG&faO 

101 FURKAT(5X,6E12.9) 

00570 

C 

» 

00580 

C 

**** COMPUTE CONSTANTS 

■ 00590 
’ 00600 

C 

KS=R04AlT 

' 00610 


RGS=RG*RO-PSYRS 

00620 


W0*S0RT<GM/PS**3) 

00630 


woe*wo*we 

' 00690 


l*i02"W0*W0 

' 00650 



00660 


N2 = WG2 + UiE2 

. 00670 


Rl a DEGRAD*RL AMDA 

' 00680 


PH*OFGPAD+PHI 

00690 


AL0*DEGRAD*ALPHA3 

• 00700 


O.MO=DEGkAD*OMEGAO 

00710 


RI=UFG»AD*R1nC 

' 00720 


CMCS(RI) 

00730 


SI-SIN (RI ) 

00790 


CP*C OS ( PH ) 

00750 


SP»SIMPH) 

00760 


Xll»CP*(WE“toO*CI) 

OG 770 


X12-CP* ( WE*CI-wG ) 

00780 


X13=k0*SP*$I 

00790 


X21=(w2-2.+W0£*Cl )*CP 

00800 


X22={W2*CI-2#*WCt)*CP 

00810 


X23 = U'02*SP*SI 

00820 


Yll a WE*CP*SI 

00830 


Y2J. a WE*Yl l 

00890 


CCONF*COS{D6GKAO+CONE ) 

00650 


NSTEP*1+600/ISTLP 

00860 


KSTEP*l+(NS-TtP-X )/2 

00870 

00880 

_C 

J J*0 

00690 

r 

**** SIAM l-F C jMPLTATlUNAL LuUP 

OC 900 

c 


00910 


DU 10 I = 1 » N S l F. P 


00920 T a IiTF?M J-KSTEP) 

00930 AL*Atu+WO*l 
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00940 ■ CAL 3 C. Q S ( A L ) ! 

00950 SAL=SIN(AL) 

00 9b0 UMEG=OMO-WE*T 

00970 CAf-,= ONEG-RL 

0C98C CGAM=COS { GAM ) 

• 00990 SGAM*SIN-(GAM) 

01000 C . L_ 

01010 C **** S/C-AIM POINT VFCTDR IN ATTITUDE REFERENCE CULKoMATES 
01020 C 

; 01030 RXgPC*(-CP*SAL*LGAl ,, i-CP*CI*CAL*SGAH+SP*SI*CAL ) 

01040 F<Y*R0*(-CP*SI*56Af'<-SP*Cn 

0105 0 RZ=PO*{-CP*CAL*CGAfHCP*Cl*SAL*SGAH“SP*Sl*SAL )+RS 

Q1Q60 C , . 

01070 C **** FIRST DERIVATIVE IF S/C AIM POINT VECTOR 
■ 01080 C 

01090 RXD=R0’M-Xll*SAL»SGAM + Xi2*CAL*CGAfl-xl3»bAL) 

01100 RYD-E0*Y11*CGAM 

OHIO RZD*R0*(-X11*CAI*SGAM-X12*SAI*CGAM-X13*CAL ) 

01120 C ' 

•01130 C **** SECOND DERIVATIVE OF A/C-AIN POINT VECTOk 
. 01140 C 

01150 RXDD = RQ*( X 2 1 * S A L *C GA M+ X22*C A L * SG A H-X2 3 *C AL ) 

01160 RYDD = R0*Y21*SGAM 
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APPENDIX D 

ORBITER DISTURBANCE MODEL 


D.O INTRODUCTION 

Four types of external disturbance torques act on the shutt 
orbitor 

& aerodynamics (T^) 

« gravity gradient (7g) 
e solar (T<.) 

a residual magnetism (T^) 

A null torque attitude is, by definition, one at which the total 
external disturbance torque is zero 


T T A + T 6 + T S + T M 0 


(D-l) 


At the expected range of orbiter altitudes for EVAL (150 to 1000 km) 
the gravity gradient and aerodynamic torques predominate. Therefore 
the total torque will be approximated by the aerodynamic and gravity 
gradient torques alone 


1 '“ T A + T G 


(D-2) 


In general, the mission-to-mission and time-to-time variability 
of the orbiter and environmental properties exceeds the neglected 
disturbance sources. For example, the orbiter is estimated to 
have a residual magnetic moment density of from ? to 5 pole-cm 
(0.002 to 0.005 Anp-m) per kilogram mass. With a typical orbiter 
mass of 97000 kg the total residual moment would be in the range 


D-l 
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184 to 485 Amp-m . The worst case condition occurs with the residual 
magnetic moment normal to the earths B field with the orbiter 150 km 
over the magnetic pole. The magnitude of the resulting torque on 
the crbiter is 


T 


M 


2 M M 



(D-3) 


3 2 

where M g /r^ = 2.9 x 10 weber/m at 150 km. The residual magnetism 
torque is thus no more than 0.0113 to 0.0281 N-m, under worst case 
conditions. This is at least an order of magnitude less than the 
nominal gravity gradient and aerodynamic torques. 

The aerodynamic and gravity gradient disturbance models used in 
this study are discussed in the following sections. 


D.l AERODYNAMIC TORQUES 

The aerodynamic torque model is based on Reference D-l , with the 
torque in body coordinates expressed as 



|qSb Cj,(a, a) 

I 

qSc C m (a, S) 

qsb C n (a, B) 


(D-4) 
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The factors entering into Equation (D-4) can be conveniently divided 
into three groupings. The influence of the external environment on 
the aerodynamic torque is summarized by the single factor 

q = dynamic pressure. 

Orbiter dimensions enter through the quantities 

? 2 

S = reference area = 2690 ft = 249.9 m 

b = wing span = 936.68 in = 23.792 m 

c = mean aerodynamic chord = 474.81 in = 12.060 m 

Finally, the influence of the orbiter attitude on the aerodynamic torque 
is summarized by the moment coefficients 

C-j = rolling-moment coefficient 
C m = pitching-moment coefficient 
C n = yawing-moment coefficient 

and the associated attitude angles with respect to the air stream 

■ a = angle of attack 
S = angle of side slip 

Each of these groupings is discussed, in turn, in the subsections that 
follow. 


D.1.1 Environmental Factors 

The dynamic pressure q is dependent on the spacecraft environ- 
ment through the relation 


q 



( D— 5) 
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where 

p = atmospheric density 

v r = relative velocity between orbiter and atmosphere 

Although dynamic pressure is nominally a function of altitude, many 
other factors can produce significant long and short term variations. 

Atmospheric Density Model 

This section contains a brief discussion of an atmospheric density 
model and an example case of density variations over an orbit. Details 
can be found in References D-l to D-6. 

Recognizing that the atmospheric density at a given altitude can 
undergo large (factor of 5 or more) variations, it is necessary to 
adopt a density model that is both realistic and tractible. The pur- 
poses of the study are served by adopting the density model in Reference 
D-l as the nominal altitude/ density model. This model, converted to 
metric units and extrapolated for altitude above 240 Km, is plotted 
in Figure D-l. Data on density variations is provided by the MSFC Modi- 
fied Jacchia Atmospheric Density Model (References D~5 and D-6). 

Measurements made with density gauges, mass spectrometers, and 
derived from satellite drag data, have revealed a number of effects 
other than altitude that result in variations in upper atmospheric 
density. All, except the last, of the following effects are included 
in the MSFC Modified Jacchia Model Atmosphere: 

« Variations with solar activity 

• diurnal variation 

» variations with geomagnetic activity 

• semi annual variation 

• seasonal-latitudinal variations 

o rapid density variations (probably associated with tidal 
and gravity waves) 
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The primary physical mechanism affecting upper atmospheric density 
is the heating and dissociation of the atmosphere caused by solar ultra- 
violet radiation. The intensity and spectral content of the radiation 
varies with solar activity with periods ranging from less than a day to the 
11 year solar cycle. Diurnal variations are caused by the rotation of the 
earth, resulting' in a density bulge around 1400 solar time at a lati- 
tude approximately equal to that of the sub solar point, and an anti- 
bulge centered around 0300 hours at about the same latitude in the op- 
posite hemisphere. The ratio of day over night density increases with 
altitude, and is largest at lower latitudes. Geomagnetic activity, 
caused by solar flares, can be correlated with large short term density 
increases. There is an average delay of 6.7 hours from the start 
of the geomagnetic storm to the time of the density perturbation. Semi- 
annual variations are not well understood, but observations indicate 
the density shows a high maximum in October, a secondary minimum in 
January, a secondary maximum in April, and a primary minimum in July. 
Seasonal-Latitudinal variations are small over the EVAL altitude 
range. 

The Quick-Look Density Model" in Reference D-6 incorporates all the 
effects in the MSFC Modified Jacchia Model in the form of equations 
and tables. The complexity of the model precludes determining the 
"worst case” variation over an orbit for EVAL. Instead, a "typical" 
low altitude (200 Km) case was worked to obtain an estimate of the 
short term density variation that can be expected on a day with low 
geomagnetic activity. The results are plotted in Figure D-2, which 
also contains a listing of the parameters used. 

^The density varies from approximately 1.8 x icf" 10 to 2.8 x 10"^° 

Kg/m , or a variation about the mean of + 20%. The variation is also 
almost sinusoidal at orbit frequency. 
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Altitude: 200 Km 

Orbit Inclination: 90 Deg (Polar) 

Orbit Hour: 1400 (Local time at ascending node) 
Sun Declination: 23.45 Deg (Summer Solstice) 
Geomagnetic Index: a^ = 6 
Solar Flux: F = T = F(t - 400) = 100 








Note that the mean density under the conditions in Figure D-2, 
about 2.3x10"^° Kg/m 3 , is somewhat less than the 3.3x10”^° K g /m 3 den- 
sity at 200 Km implied by Figure D-l . Rather than being an error, this 
discrepancy illustrates the wide range of values possible under varying 
environmental conditions. For purposes of consistency, the analysis in 
the main body of this report always assumes the short term density varia- 
tions occur about the nominal density values shown in Figure D-l. That 
is, the assumed atmospheric density model has the form 

p = p Q + Ap (D-6) 

where 

p Q = nominal density (Figure D-l) 

Ap = density variation 


Relative Velocity 

The relative velocity between the orbiter and atmosphere is in 
vector form 


where 


v 


r 



(D-7) 


v Q = orbiter inertial velocity vector 
= earth rotation rate vector 
? s = orbiter position vector 


The magnitude and direction of 7 thus depends on the orbiter altitude, 
orbit inclination, and orbit position. As an approximation, the ro- 
tation of the atmosphere with the earth will be neglected, resulting 
in 


v„ = r x D 
o o s R 


(D-8) 
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where 


1/2 

= orbit rate 

5 3 2 

GM = earth gravitational constant = 3.986 x 10 km /sec 

r s = r 0 + h 

h = orbiter altitude 

r Q - earth radius = 6371 km 

x R = unit vector in direction of orbit velocity (Appendix C) 

The maximum error in the direction of v introduced by this approximation 
i s 



tan 


-1 




7.27 x 10 5 
1.1 x 10~ 3 


= 3.8 degrees 


The maximum relative error in the magnitude of v r> over the altitude 
range where aerodynamic torques are significant, is w e /u Q = 0.066 = 6.6°/. 

Dynamic Pressure Model 

The nominal dynamic pressure, with the above approximations, is 
plotted as a function of altitude in Figure D-3. As a consequence of 
the atmospheric density model and Equation (D-5), the dynamic pressure 
model has the form 

9 = q 0 + Aq (D-9) 

where 

q Q = nominal dynamic pressure (Figure D-3) 

Aq = dynamic pressure variation 


D-9 








D.1.2 Dimensional Factors 


The orbiter dimensional factors, obtained from Reference D-l are 

S = reference area = 2690 ft 2 = 249.9 m 2 

b = wing span = 936.68 in = 23.792 m 

c = mean aerodynamic chord = 474.81 in = 12.060 m l 


Combining the dynamic pressure with the appropriate orbiter dimensions 
yields the nominal lateral (roll/yaw.) aerodynamic torque constant qSb 
and the longitudinal (pitch) constant qSc. These' are plotted as a func- 
tion of altitude in Figure D-4. 


D.1.3 Attitude Factors 

The attitude dependence of the aerodynamic torques appears through 

the moment coefficients C (a,e), C (a,e) and C (a,e). The correspondence 

J6 in n 

of these coefficients to the orbiter body axes and the direction of the 
relative wind is illustrated in Figure D-5. With the approximation of 
Equation (D-7), the relative wind is along Xp and the angle of attack is 


a = tan 


-1 


X R • 2 B 


u 


X R * X B 



L_ 


(D-10) 


while the angle of side slip is 


6 = sin 


[x D • y R ] = sin 


-1 


0 21 ] 


(D-l 1 ) 


where a-j-j, a^, and are elements of the A^p matrix. Table D-l 
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Figure D-4. Nominal Orbiter Aerodynamic Constants 
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Note: Senses indicated 
are positive. 




DIRECTION OF 
RELATIVE WIND 


Figure D-5. Definition of Aerodynamic Moment Coefficients 





Table D-l. Angles of Attack and Sideslip for Small 
Perturbations from Nominal Orientations 
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lists the correspondence between « and 3 and the small perturbation 
angles 0 , ? with zero offset angles ^ , 0 o , for the nominal 
orientations outlined' in Section 8. 

The values of the moment coefficients are plotted and tabulated 
in Reference D-l for 0 £ a <. 350 degrees in 10 degree increments 
and for 0 £ 3 £ 180 degrees in 15 degree increments. The plots are 
reproduced here as Figures D-6 to D-8. Because the orbiter is sym- 
metrical with respect to the - Zg plane, the moment coefficients 
for - 180 <_ 3 <_ 0 are given by 


C z (ct C z (a, 3) (D-l 2) 

C m (a * ’ S) = C m (a * S) (D-13) 

C n (a , -B) “ - C n (a, 5) ' (D-l 4) 


It should be noted that the moment coefficients in Reference D-l 
were computed with the located 65t of the orbiter length back from 
the nose and with the payload bay doors closed. Data with the doors 
open is not available at this time. As shown in Figure D-9 the 
payload bay doors fold down over the wings with 1 ittle- protrusion . 

The main effect of opening the doors is therefore to replace the 
smooth upper shell of the orbiter with the contents of the payload 
bay. The mission dependence and complex shape of the payload 
bay contents makes a general evaluation of the impact on the moment 
coefficients a difficult task. 
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(b) [3 = 90° to 180° 


Figure D-6. Variation of Rolling Moment Coefficient With Angle of Attack 
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D. 2 GRAVITY GRADIENT TORQUES 


The gravity gradient disturbance torque model is 



r x I • r 


(D-l 5) 


where 

to = orbit rate 
o 

r = unit vector from orbiter to geocenter 
I = orbiter inertia dyadic 


In terms of the body axis inertias and elements of the Ag^ 
the torque is 


matrix 






Z yy) a 23 a 33 + T xy a 13 a 33 " *xz a l3 a 23 + l y z ' a 33 
l zz^ a '[3 a 33 + l yz a l3 a 23 “ T xy a 23 a 33 + l xz (a 13 
r xx^ a 13 a 23 + T xz a 23 a 33 " l yz a 13 a 33 + T xy ' a 23 


a 2 ' 
a 23' 

a 2 ' 
a 33 


The orbiter mass properties, in particular the products of inert - ' a, 
are mission dependent. For this study the typical values from Appendix A 
will be assumed. 


(D-l 6) 
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I = 1.24 x TO 6 Kg-m 2 

A A 

I = 9.39 x 10 6 Kg-m 2 

I = 9.72 x 10^ Kg-m 2 

**y = T xz " V " 0 

The gravity gradient torque therefore reduces to 

( ! zz ~ T yy ^ a 23 a 33 
^ 1 xx “ ! zz^ a l 3 a 33 
^yy " W a l 3 a 23 



(D-17) 


The factor 3^ is related to the altitude h in km by 


* 

3 '2 = 3 GM = 1.1958 x 10 6 

~° (r t h) 3 (6371 + h) 3 


rad 


sec 


(D-18) 


This relation is plotted in Figure D-10. Whereas doubling the altitude 
from 180 km to 360 km reduces the nominal aerodynamic constants by more 
than two orders of magnitude (Figure D-4) the gravity gradient constant 
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3u is reduced by only 7.8%. In view of this, and the fact that the 
o . 

aerodynamic constants are subject to wide variations from environmental - 

factors, the altitude dependence of the gravity gradient torque can 

2 -6 

be safely neglected at low altitudes. The value 3j = 4 x 10 

2 0 
rad/sec , corresponding to h * 300 km is therefore assumed in deriving 

the gravity gradient torques for small perturbations frormthe nominal 

orientations listed in Table D-2. 
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Gravity Gradient Torque (N-m) 
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Torque values given hold within +_ 10% for altitudes in the range 150 £ h <_ 550 km 
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APPENDIX E 


TILTED DIPOLE MAGNETIC FIELD MODEL 


For present purposes, an adequate description of the earth's 
magnetic field is the tilted dipole model. In attitude reference 
coordinates, the field is 





(r + h) 3 


-cos s sin i cos a + sin e (sin a 
cos e cos i + sin e sin i cos u 


sin u + cos i cos a 


-2cos s sin i sin a -2 sin £ (cos a sin u -cos i sin 


where 

M e = 8.056 x 10 10 gauss - km 3 
r Q = 6371 km 
h = orbiter altitude 
e = dipole tilt »11 deg 
i - orbit inclination 
a = « 0 t 
u = « e t 

3 1 

0 ) = [GM/(r 0 +h) d ]2 

GM = 3.986 x 10 5 km 3 /sec 2 

ti) e - 7; 27 x 10 -3 rad/sec 
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cos u) 


a COS UX 



The field in body coordinates is found by transforming through the 



matrix 


A BR b r 


(E-2) 


The "most unstable" axis for control purposes is the y-body 
axis in Orientation 3C, and this will be adopted for sizing purposes. 
The A BR matrix for orientation 3C is approximately 



Using this transformation, the body components of B* are plotted in 

Figures E-l and E-2 for h = 2Q0 km and i = 30 and i = 90 degrees. 

The corresponding plots of vs b z are in Figures E-3 and E-4. The 

field components are plotted in the common units of Gauss. The more 

proper MKS units of weber/m can be obtained by multiplying the field 

-4 

strength in Gauss by 10 . 
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APPENDIX F 


MOMENTUM EXCHANGE DEVICE TRADEOFF 


F.O INTRODUCTION 

This appendix documents the tradeoff study of momentum exchange 
actuators for the orbiter null torque attitude stabilization system for 
EVAL. The candidate actuator types are the reaction wheel (RW), single 
gimbal control moment gyro (SG CMG) and double gimbal control moment gyro 
(DG CMG). 

F. 1 ACTUATOR REQUIREMENTS 

The actuator requirements fall into two main categories, namely, con- 
trol requirements and accommodation requirements. Each of these groups will 
be discussed below. 

t 

F . 1 . 1 Control Requirements 

The primary control requirements on the momentum exchange actuators are 
adequate control torque and momentum storage capability for capture from the 
orbiter RCS limit cycle and during extended operation in a variable distur- 
bance environment. Table F-l, based on the analysis and simulation in 
the main body of this report, summarizes the estimated control torque 
and momentum requirements under a variety of conditions. The largest 
requirements, by far, occur for capture from the nominal orbiter RCS 
limit cycle rates. Note that the capture requirements imply an essentially 
spherical momentum envelope. 

An additional important control requirement is continued, although pos- 
si'bily degraded, three -axis control capability with one unit failed. This 
implies that a minimum of four RWs or SG CMG’s, or three DG CMG's are re- 
quired. Further desirable control features are simple actuator control 
(steering) laws and freedom from control singularities. 
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Table F-l.. Estimated Control Torque and Momentum Requirements 


Case 

Control Requirements 

Mode 

Orientation 

Altitude 

(Km) 

Conditions 

Torque (N-m) 


Pitch 

Roll /Yaw 

Pitch 

Roll /Yaw 

Capture 
from RCS 
Limit Cycle 

Nose 

Forward 

(worst 

case) 

180 to 
1000 

Nominal Rate Limit Cycle 
(0.01 deg/sec) 

12 

8 

2700 

2700 

0.1 x Nominal Rate .Limit 
Cycle (0.001 deg/sec) 

1.3 

1 

500 

' 500 

Follow 

Null 

Torque 

Attitude 

(practical 

control 

law) 

Nose 

Forward 

200 

+ 20 % Atmospheric 
' Density Variation 

at Orbit Frequency 

0.12 

★ 

105 

* 

Nose 

Down 

200 

1.9 

* 

1740 

★ 

250 

0.34 

* 

288 

* 



*Nomina11y zero; depends on unloading policy, products of inertia, pitch motion, etc. 























F.1.2 Accommodation Requirements 

* The orbiter payload size, weight, and power accommodations are subject 
to the limitations outlined in Reference F-2. Since the size/weight/power 
used by the momentum exchange actuators must be deducted from the experiment 
allowance, it is essential that these be minimized. 

Size 

Depending on the spacelab configuration used on a particular mission, 
the volume available for payload, outside the spacelab module, ranges from 
about 32 meter (1 pallet segment) to 184 meter {5 pallet segments with 
overhang). As a goal, the actuator volume should be limited to about 10% 

3 

of a pallet segment or 3 meter . 

Weight 

Again depending on the spacelab configuration, the mission dependent 
payload weight allowance ranges from about 5800 to 9350 kg (12,760 to 20,570 
pounds). As a goal, the actuator weight should be limited to about 10% of 
the minimum weight allowance or 580 kg (1276 pounds). 

Power 

Electrical power for the orbiter is provided by hydrogen/oxygen fuel 
cells. The orbiter baseline provides only 50 KWh of electrical energy for 
spacelab; however, an additional 840 KWh "energy kit" is included in the 
spacelab weight, providing a total of 890 KWh for the spacelab and its pay- 
loads. The heat rejection capability of the orbiter limits the total 
space! ab/payload electrical power dissipation to 7 KW continuous, and 12 KW 
peak (15 minutes out of 3 hours). Of this total, the power available to the 
payloads and mission dependent equipment ranges from 2.57 to 5.37 KW con- 
tinuous and 7.0u to 10.2 KW peak. As a goal, the actuator power, at peak 
torque, should be limited to about 10% of the minimum continuous available 
payload power, or about 250 watts. 
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F. 2 SIZING ANALYSIS 


Estimates of the weight, size, and power of each type of actuator (RW, 

SG CMG, DG CMG) are derived in this section, using data in Reference F-l. 
Three combinations of torque and momentum are included in the analysis to 
span the requirement range implied by various sets of operational conditions, 
as summarized in Table F-2. 

The sizing analysis is in terms of individual actuator units of each 
type capable of producing equal magnitudes of peak torque and momentum. The 
actual control capabilities of a set of actuators depends not only on the 
individual unit characteristics but also on the number of units used and 
the mounting configuration {orthogonal , .skewed, pyramid, etc.). 

The reaction wheel is fundamentally a servo motor with a large inertia 
rotor. On command, motor torque accelerates or decelerates the rotor, 
changing its angular momentum. The reaction to this torque is transmitted 
through the unit mounting pads to the vehicle structure. The wheel motor 
can continue to apply torque to the vehicle until a maximum speed somewhat 
below the motor no-load- speed is reached. At this speed, the maximum amount 
of angular momentum is stored in the reaction wheel. Because the motor 
torque and rotor rate vectors are parallel in a reaction wheel, the electri- 
cal power input required to produce a given torque is relatively large . 

The electrical power (watts) required to produce a. torque T (N-m) with 
the rotor spinning at an angular rate w (rad/second) with a motor power ef- 
ficiency n can be expressed as P^ = T w/n. This relation restricts reaction 
wheels to low torque applications. 

The SG CMG operates with a fixed momentum magnitude with one. gimbal 
rotational degree of freedom and the output torque T * « x ^ Is thus con- 
strained to lie along the line normal to both the gimbal axis and the wheel . 
momentum vector F y/ . This is the SG CMG output torque axis, and since it 
lies normal to the gimbal axis there is no component of output torque about 
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Table F-2. Actuator Torque/Momentum Combinations for Sizing Analysis 


Actuator Characteristics J 

Operational Condition Restrictions 

Torque 

N~m 

Mbmentum 

N-m-sec 

Maximum 

RCS Limit Cycle Rate 

Minimum Altitude (km) for 
Orb iter Orientation 

(ft-lb) 

(ft-lb-sec) 

deg/sec 

Nose Forward 

Nose Down 

1.36 

678 

0.001 

180 

230 

0) 

(500) 

(0.1 x Nominal) 

13.6 

1360 

0.001 

180 

210 

(10) 

(1000) 

(0.1 x Nominal) 

13.6 

2710 

0.01 

. 180 

200 

00) 

(2000) 

(Nominal) 




the gimbal torquer axis. The SG CMG output torque is thus limited only by 
the radial load carrying capacity of the gimbal and momentum wheel bearings. 
The gimbal torquer must only accelerate the gimbal and since there are no • 
large reaction torques about the gimbal axis, relatively large gimbal rates 
are possible. The SG CMG is thus a torque multiplier, i.e., small torques 
applied about the gimbal axis to establish rate can produce large torques 
about the output axis. In broad, terms the single gimbal CMG is well suited 
to high torque applications, but has a constrained momentum storage capabil- 
ity. 


The momentum envelope of one SG CMG is the plane normal to the gimbal 
axis and the envelope of an actuator cluster is dependent on the mounting 
configuration of the individual actuators. While the basic cluster can be 
arranged to give a more or less regular momentum envelope, a single CMG 
failure strongly distorts the envelope and causes quite complex software 
problems in the steering laws. Another severe problem with a SG CMG cluster 
is the existence of singularity surfaces within the momentum envelope oc- 
curring when the torque output vectors of the individual actuators are co- 
planar, i.e., no control torque is possible about the axis normal to the 
plane. A total' analytic solution for the location of these surfaces and 
general software for singularity avoidance does not exist, and the normal 
strategy is local avoidance of the surfaces along with varying amounts of 
excess momentum storage capability over the basic requirements. 

The DG CMG has two rotational degrees of freedom, one about each of the 
two perpendicular gimbal axes. Motion of the two gimbals can produce an ef- 
fective gimbal rate w g anywhere in the plane normal to the wheel momentum 
H" w thus the output torque T = x H w can lie anywhere in this plane which 
can then be called the DG CMG output torque plane. However, the output 
torque always has a component along one of the gimbal axes, therefore, it is 
limited by torque capability of the gimbal torquers. Again in broad terms, 
the double gimbal CMG is more suited to applications with moderate torque 
requirements, but has a more flexible (i.e., optimum) momentum storage capa- 
bility. 
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The D6 CMG momentum envelope is spherical in the absence of gimbal 
stops and even with limited gimbal rotation the envelope of a single actua- 
tor is a portion of a sphere. The envelope of a DG CMG cluster of any num- 
ber of actuators is thus also spherical and the failure of a single actuator 
while decreasing the envelope size does not change its spherical character. 
Desaturation strategies can thus be devised to decrease or control the 
stored momentum with the guarantee that the momentum vector is always moving 
away from the envelope in both normal and failure modes of operation. DG 
CMG singularity occurs only when the individual actuator momentum vectors 
are col inear and this situation can be avoided through relatively simple 
logic involving redistribution of the individual momenta through the use of 
the excess degrees of freedom provided when more than a single actuator is 
in use. 

F.2.1 RW Sizing 

For given saturation momentum and peak output torque requirements, op- 
timization of a RW design involves basically the selection of the rotor sat- 
uration speed that minimizes a weighted combination of weight, size and 
power. The total weight is the sum of three components 

• The rotor and housing 

e The motor 

e The power supply and drive electronics 

The first two of these components are presented in graphical form in Refer- 

ence F-l, along with size and power relations. The weight of the power supply 
and drive electronics is largely related to the peak motor power requirements, 
but is subject to wide variations (based on cooling constraints) with re- 
spect to the required torque time history and duty cycle. A reasonable value 
for the power supply and drive electronics weight is 1.0 pounds/watt of peak 
power. 
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Figure F-l summarizes the weight/power tradeoff for H = 500 foot-pound- 
seconds, T * 1 foot-pounds. The straight solid line represents the peak 
power, while the upper solid curve represents the total RW + electronics 
weight with a 1.0 pound/watt power supply/drive electronics weight penalty. 
The minimum weight of 400 pounds occurs at a - 400 RPM , with a power input 
at peak torque p = 125 watts. For the required four units, the total weight 
is 1600 pounds (727 kg) and the total power is 500 watts. Both these fig- 
ures exceed the previously stated goals: maximum weight = 1276 pounds 
(580 kg,), maximum power = 250 watts. 

The tradeoff was reworked with an optimistic power supply/drive elec- 
tronics weight penalty of 0.1 pound/watt and the result is plotted as the 
lower solid curve in Figure F-l. The minimum weight per unit is reduced to 
less than 200 pounds for 1300 < a < 2400 RPM, but the power to run at these 
speeds ranges from 400 to 720 watts. The power can be reduced to .102 watts 
by reducing the wheel speed to 340 RPM and increasing the weight per unit to 
one-quarter of the total weight goal or 319 pounds (145 kg). The power for 
four units is 408 watts, exceeding the power goal. It is thus apparent that 
simultaneous satisfaction of the weight and power goals is not feasible with 
reaction wheels. The design coming closest to meeting the weight/power 
goals is indicated by the small circles in Figure F-l. The diameter of each 
of these RW's is 44 inches. If it is assumed that the volume occupied by 
the RW plus electronics corresponds to a cube with the RW diameter as a side, 

3 -5 O 

the volume is 85,184 inch per unit or 340,736 inch = 5.58 meter per 
cluster of four. This is almost twice the volume goal of 3 meter 3 . 

Figure F-2 shows the tradeoff curves for the H - 1000 foot-pound- 
seconds, T = 10 foot-pounds reaction wheel units. It is clear that regard- 
less of wheel speed, the size/weight/power goals will be far exceeded, and 
no satisfactory RW design exists. The H = 2000 foot-pound-seconds, T = 10 
foot-pounds units are even more futile and no curves are provided for these. 


F-8 



Weight (lb) or Power (watts) 



Figure F-l . Reaction Wheel Weight/Power Tradeoff 
H = 500 ft-lb-sec, T = 1 ft-lb 
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Reaction Wheel Weight/Power Tradeoff 
H = 1000 ft-lb-sec, T = 10 ft-lb 
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With the exception of the initial run-up periodj the CMG rotor operates 
at fixed speed and the input power to the spin motor only has to make up for 
losses. It is, therefore, practical to run the CMG rotor at a higher speed 
than with a RW, and all the CMG design data in Reference F-l is based on an 
optimum rotor speed of 11,700 RPM. The total weight and power is relatively 
insensitive to rotor speed over the range from 8,000 to 14,000 RPM, however. 

The SG CMG power requirement, after run-up, consists of the rotor run- 
ning power and the gimbal torque power. The rotor running power is 14 watts 
for the 500 foot-pound-seconds units, 20 watts for the 1000 foot-pound- 
seconds case, and 26 watts for the 2000 foot-pound-seconds .CMG's. The torquer 
power is in all cases negligible for the required output torque levels. The 
SG CMG total weight consists of the CMG rotor/housing/gimbal weight plus the 
power supply/drive electronics weight. The power/weight/ size values for the 
three momentum/ torque combinations are summarized in Table F-3. 

F.2.3 DG CMG Sizing 

The sizing analysis of the DG CMG is very similar to that of the SG CMG, 
with the basic difference of the addition of an extra gimbal and its associ- 
ated drive electronics. This increases the weight and size of each unit, 
compared to a SG CMG with the same torque and momentum capability. However, 
fewer units are required in the cluster. The larger reaction torques on the 
gimbals in the DG CMG increases the required gimbal drive torques so that 
the gimbal drive power is no longer negligible. TheDG CMG sizing results, 
derived from Reference F-l, are summarized in Table F-4. 



Table F-3. SG CMG Summary* 


Torque 

(ft-lb) 

Momentum 

(ft-lb-sec) 

Power (watts) 

Height (lb) 

Size 

Spin 

Motor 

Gimbal 

Drive 

Total 

CMG 

Electronics** 

Total 

Diameter 

(in) 

Vol ume 
(in 3 ) 

1 

500 

14 

< O.l 

14 

70 

14 

84 

24 

13,824 

10 

1000 

20 

< 0.1 

20 

95 

20 

115 

26 

17,576 

10 

2000 

26 

< 0.1 

26 

130 

26 

156 

30 

27,000 


*Per unit, four required 

**1 Ib/watt 
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Table F-4. DG CMG Summary* 


Torque 

(ft-lb) 

Momentum 

(ft-lb-sec) 

Power (watts) 

Weight (lb) 

Size 

Spin 

Motor 

Gimbal Drive 

Total 

CMG 

El ectronics** 

Total 

Diameter 

(in) 

Volume 

(in 3 ) 

1 

500 

14 

m 

15 

85 

15 

100 

27 

19,683 

10 

1000 

20 

2.5 x 2 = 5 

25 

140 

25 

165 

32 

32,768 

10 

2000 

26 

2.5 x 2 = 5 

31 

. 175 

31 

206 

38 

54,872 


*Per unit, three required 
**1 lb/watt 

























F.3 CONCLUSIONS 


The results of the momentum exchange actuator sizing analysis are sum- 
marized in Table F-5 using metric- units. Reaction wheels clearly fail to 
meet the requirements in all three size categories, and can, therefore, be 
eliminated from further consideration. Both single gimbal and double gim- 
bal CMG's, on the other hand, are consistent with the assumed torque/ 
momentum requirements and size/weight/power limitations. The size/weight/ 
power figures, on a total system basis (four single gimbal or three double 
gipibal CMG's) are sufficiently close that additional criteria must be con- 
sidered to make a final decision on the actuator type. 

For a given momentum, the cost of an individual double gimbal CMG is 
higher than a single gimbal CMG. However, because fewer double gimballed 
units are required, the total system costs are likely to be comparable. 
Likewise,’ it is technically feasible to build either single or double gim- 
balled CMG's in the momentum range required. 

On performance grounds, the double gimballed CMG has several advantages 

© Better match to the required nearly spherical momentum 
envelope, especially with one unit failed. 

e . Relatively simple steering law. 

e Less problem with singularities. 

The double gimballed CMG, therefore, appears to have the best combination 
of properties- for the EVAL mission. 
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Table F-5. Momentum Exchange Actuator Sizing Summary 


Capability per Actuator Unit 

Type 

Humber 

in 

Cluster 

Total 

Requirement for Cluster 

Torque 

(N-m) 

Momentum 

(N-m-sec) 

Weight (kg) 
Goal = 500 kg 

Power (watts) 
Goal = 250 watts 

Volume (m^) 
Goal = 3 m 3 



RW 

4 

580 

408* 

5.58* 

1.36 

678 

SG CMG 

4 

153 

56 

0.91 



DG CMG 

3 

136 

45 

0.97 



RW 

4 

> 964* 

•> 1200* 

> 7.25* 

13.6 

1360 

SG CMG 

4 

209 

80 

1.15 



DG CMG 

3 

225 

75 

1.61 



RW 

4 

* 

* 

k 

13.6 

2710 

SG CMG 

4 

284 

104 

1.77 



DG CMG 

3 

281 

' 93 

’ 2.70 


♦Does not meet goal 


RW = reaction wheel 

SG CMG = single gimbal control moment gyro 
DG CMG = double gimbal control moment gyro 
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